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ABSTRACT 


The goal of this thesis effort was to develop an 
interactive electromagnetic scattering algorithm for 
personal computer use in the analysis and computer aided 
design (CAD) of scattering structures. Extensive 
interactive computer aided design type features have been 
added to a previously developed finite element computer 
algorithm which considers scattering from inhomogeneous 
axisymmetric objects. These computer aided design features 
allow structures composed of multiple dielectric and/or 
magnetic materials to be graphically constructed or 
modified by the user by way of either a keyboard or mouse 
interface. Bistatic scattering patterns are graphically 
presented upon completion of the computation for a variety 


of user specified incident field orientations. 
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DN TNURODUCTTON 


The study of resonance region electromagnetic 
interaction with arbitrarily shaped and often inhomogeneous 
penetrable objects has increased dramatically in the past 
few years. This new concentration of effort is driven by 
the likelihood of expanded implementation of high frequency 
over-the-horizon radars, very high frequency satellite-borne 
radars and medium frequency surface wave radars. Such 
radars have the potential to observe "stealth" type 
platforms whose low radar cross section (RCS) design was 
limited to more conventional radar frequencies, above 500 
MHz. In addition, these new studies of electromagnetic 
interaction have proven to be of benefit in other areas 
such as propagation through rain, power absorption by 
biological tissues, focusing by microwave lenses, radiation 
from reentry vehicles and electromagnetic pulse (EMP) 
coupling into vehicles, installations and weapon systems. 
[Refs. 1 and 2] 

As frequencies become sufficiently low, structures 
having a low radar cross section may exhibit strong 
resonance characteristics. Design of such structures is 
complicated by these low frequency resonance 
Characteristics, due to the lack of accurate analytical 


methods. Even attempts at structure design utilizing 


"shaping" and "edge treatments" are ineffective at providing 
a highly reduced radar cross section at resonant 
frequencies. At higher frequencies, "asymptotic" 
RAUS: are available to rectify these problems (Ref. 
ΙΕ A more comprehensive approach is needed, one which 
allows the designer to rigorously evaluate the interaction 
of the scatterer with the illuminating field. At the same 
time, such a tool should be highly interactive with the 
user. 

The subject of this thesis is the development of a 
generalized, interactive, finite element scattering 
algorithm to be used in the analysis and computer aided 
design of scattering structures. In this, the first stage 
of effort, the programs developed for use on a personal 
computer calculate electromagnetic scattering Бу 
axisymmetric inhomogeneous structures. These structures are 
composed of user specified lossy or lossless penetrable 
(i.e., dielectric or magnetic) materials. The resulting 
programs provide significant improvement in computational 
speed and versatility over the current integral equation 
based computer codes being used. The physical structure can 
be quite general, other than the requirement for rotational 
symmetry. This general construction avoids the 
computational restrictions found in surface integral 
equation formulations, where, as frequency and geometric 


complexity increase, the time required to obtain a solution 


lncreases. The flnite element algorithms, based upon a 
volume variational formulation, are relatively immune to 
structural complexity driving the computation time. Another 
advantage of these alqorithms is their method of matrix 
manipulation. The finite element algorithms produce very 
"Sparse" matrix structures, whereas the integral equation 
formulations produce "full" matrices (Ref. 4]. Finally, 
when considering the construction of realistic material 
configurations, the computational time required utilizing 
the finite element algorithms can be substantially less than 
the time required by the use of integral equation 
formulations  [Ref. 5]. 

The development of the axisymmetric finite element 
program is based upon  Morgan's previous effort  [Ref. 1], 
which restricted structures to non-metallic configurations. 
Major programming modifications to improve computational 
efficiency and to incorporate computer aided design features 
such as interactive graphics are now integral parts of the 
finite element computer program. 

A brief overview of the finite element scattering 
algorithm, which constitutes the basis of this thesis, is 
presented in Chapter II. In addition to the appropriate 
background information, the concept of coupled azimuthal 
potentials, the finite element algorithm and the unimoment 


method are reviewed. 


Having established, in Chapter IL, the necessary theory 
utllized in the construction of the original version of the 
computer program EMSCAT, Chapter III presents methods to 
construct progran input data. Included in this chapter are 
discussions of computer programs written specifically for a 
given geometric shape, data file input and keyboard input. 
Basic flow charts provide an overview of the input-output 
interactions of the all computer programs. 

The first half of Chapter IV discusses the extension of 
the methods used to construct program input data to those 
methods that utilize graphics. Of special interest is the 
incorporation of computer aided design features. The second 
half of this chapter addresses the user interactive output 
graphics features provided. These features allow the user 
to obtain one of four different graphs, either on the 
computer screen or printer, of the EMCAD output data of 
cholce. 

In Chapter V, the validation of modifications and 
improvements to EMCAD is presented. Graphical comparisons 
are made between the output of the modified version of 
EMSCAT versus both theoretical calculations and actual 
experimental data. In addition, validation of new programs 
written specifically to be used in conjunction with EMCAD is 
presented. 

The brief conclusions in Chapter VI will summarize the 


work previously presented, as well as provide 


recommendations for future modifications to the interactive 
package of computer programs being used. 

The Appendices provide copies of the source code for 
all the E programs utilized in this interactive 
package. A sample listing of the output of EMCAD is 
attached. As an added feature, user manuals with hints, 
warnings and directions for all of the major computer 


programs are included. 


Ll. FINITE ELEMENT SCATTERING ALGORITHM 


A. BACKGROUND 

In recent years, great attention and effort has been 
directed toward the attempt to obtain the numerical 
electromagnetic scattering solutions for inhomogeneous 
material structures. Early efforts dealt primarily with 
improving the efficiency of existing computational 
algorithms. These algorithms utilized surface integral 
equation formulations, which required the inversion of large 
“full matrices. Accordingly, methods using these 
formulations were computationally slow while being highly 
restricted in the level of inhomogeneity that could be 
handled.  [Refs. 1 and 2] 

À major advance in matrix manipulation occurred with 
the development by Mei (Ref. 6] of the unimoment method. 
This new method decoupled the interior and exterior regions 
associated with radiation or scattering problems, resulting 
in uniformly banded system matrices for the interior region 
solution. The primary advantage of uniformly banded 
matrices is the capability to invert them directly using 
fast banded matrix inversion techniques, such as the Riccati 
transform. Mei's breakthrough opened the door for the 
development of applications such as the finite-difference 


solution of radiation by a biconical antenna, and the 


finite-element solution of scattering by two-dimensional 
inhomogeneous dielectric cylinders. (Refs. 1 and 2] 

Attempting to improve upon these applications, Morgan, 
Chang and Mei (Ref. 7] conducted research with the goal of 
extending the unimoment technique to problems involving 
three-dimensional, axisymmetric  inhomogeneous structures. 
During the initial period of this research, another major 
breakthrough occurred. This breakthrough was the discovery 
of the coupled azimuthal potential (CAP) formulation 
[Ref. 7]. The САР formulation became the basis for 
construction of several numerical algorithms. These 
algorithms incorporated the capability to perform scattering 
computations for inhomogenous axisymmetric metallic and 
penetrable structures  [Ref. 2]. 

Utilizing a triregional unimoment method, these 
numerical algorithms were incorporated into a versatile 
finite-element computer program. The triregional unimoment 
method was chosen because of its ability to minimize the 
bandwidth of the global system matrix, and its ability to 
minimize the number of nodes in the finite-element mesh. 
The finite-element computer program generates scattering 
solutions for an arbitrarily shaped inhomogeneous’ lossy 
penetrable body of revolution being impinged upon by 


multiple incident fields. [Ref. 2] 


Β. COUPLED AZIMUTHAL POTENTIALS 

The solution of electromagnetic field problems 
involving generalized axisymmetric penetrable structures 
using the classical TE and TM potential formulations are 
restricted to a handful of configurations, and have at most 
inhomogeneities in one-dimension. By utilizing the coupled 
azimuthal potential (CAP) formulation, solutions are valid 
in lossy inhomogeneous isotropic rotationally symmetric 
media, as shown in Figure 2-1, where the complex material 
parameters are invariant to the azimuthal coordinate. 
[Ref. 1] 

Utilizing normalized circular cylindrical coordinates, 
as also shown in Figure 2-1, the first step in developing 
the necessary equations is to decompose the total 
electromagnetic field into azimuthal modes through the use 


of an exponential Fourier series, 


0 
Fr a e (R 2) exp(jn$) aon 
π. Ὁ 
and 
= Ф амы. 
ngHC(R 2,8) s 2 h (R42) expCjnd) ο ο) 
m=-@ 


531УМ1049003 1VOISQONI 1۸3 ΥΠ ΒΙΟ Q03ZI ЛУИЧЯОМ 
МІ ИПІОЗИ JIYLAWWAS-ATIVIXY 1-с 7913 


ТЫ 


9 


Y 
LN 


| (е-)3 
πό νου 


ae, 


> 


where 





RET) = СКОР 02) (255 
with 
2n 
к, = | (29) 
0 
апа 
л, ІЛ ee ου 


Two uniformly continuous coupled azimuthal potential 


functions 
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can be used to generate the modal electromagnetic field 


pOcCcors, 
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while the multiplicative function is given by [Refs. 1, 2, 


and 7] 
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For a two-dimensional cut in the (R,Z) plane of Figure 
2-1, the coupled azimuthal potential functions satisfy a 
OE self-adjoint system of second-order partial 
differential equations. Utilizing either Euler-Lagrange or 
stationary-theorem approaches, it can be shown that there is 
an easily managed variational principle which can be applied 
in addition to the partial differential equation 
formulation. 

This variational formulation is based upon the 
existence of a stationary function which manifests itself as 


a surface integral over the planar cross section, S, 


F = LCR, 2, v ;V 
с | 
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where the Lagranqian L, which is unique to within an 
arbitrary constant multiplier and arbitrary independent 


additive function, is of the form (Refs. 1 and 2] 
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e. FINITE-ELEMENT ALGORITHM 

Having presented an overview of the coupled azimuthal 
potential formulation, it is now appropriate to explain its 
application to the electromagnetic scattering problem. 
Figure 2-2 illustrates the finite-element mesh geometry 
employed by the numerical algorithm. The (R,Z) planar 
configuration shown is named the semi-annular conformal 
(SAC) mesh [Ref. 1]. This mesh is comprised of numerous 
linear triangular elements, which conform to the surface of 
the axisymmetric scattering body. This surface conformation 
simplifies analysis of normal derivative discontinuities at 
the surface boundary. The semi-annular conformal mesh 


provides a distinct computational advantage to those 
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FIG. 2-2. SEMI ANNULAR CONECRMANE 
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INHOMOGENEOUS BODY OF REVOLUTION 
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scattering structures which can support κ. falrly large 
homogeneous spherical core [Ref. 2]. 

Figure 2-2 graphically demonstrates the regions of 
interest Mhen calculating the solution of the scattering 
problem. The homogeneous spherical core occupies Region I, 
the inhomogeneous cross sectional surface of the 
axisymmetric penetrable scattering body is in Region II, and 
an unbounded exterior region is in Region III. The surface 
of interest is shown with the capability of being both 
banded and layered. The layers are depicted as "films" that 
run the entire length of the surface, while bands are 
depicted as isolated "pockets". The possibility exists, as 
shown, that all the bands and all the layers have differing 
material parameters, thus establishing a structure that is 
truly inhomogeneous. Note that the surfaces of Region I and 
Region III bound the surface of interest located in Region 
II. These surfaces, r = a and r = b, can therefore be used 
to specify Dirichlet boundary conditions [Ref. 2]. The 
numerical solution for the electromagnetic fields in the 
finite-element domain (Region II) can be computed for each 
of a set of boundary conditions. By requiring continuity of 
the electromagnetic fields at the interfaces where г = гү 
and r- r5 , the unknown expansion coefficients for the core 
and scattered fields can be obtained (Ref. 1]. This is the 


essential concept behind the unimoment method. 
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In Region II of Figure 2-2, the semi-annular conformal 
mesh is composed of nodes, with each interior node 
Surrounded by six triangular regions, known as "elements". 


The azimuthal field components of the mesh are described by 
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where u,(R,Z) represents the area coordinate linear pyramid 
basis functions. Each of these basis functions has linear 
variation in the six triangular elements around each node, 
while maintaining a value of one at the central node and a 
value of zero at all bounding nodes. Each set ог е „(п) апа 
hq (n) represent one of the N complex nodal values found in 
the azimuthal field. To initiate the variational finite- 
element technique, the linear shape-function expansions in 
Equations 2.12 and 2.13 are substituted into the surface 
integral stationary function F described in Equation 2.10. 


Imposing the null derivative condition with respect to each 
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of the unknown coordinates representing the interior nodal 
fields, while constraining the coordinates of the boundary 
nodal fields, allows the computation of the stationary 
point Bene tunceion F. The resulting linear system for 


the nodal fields is described by 
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where n = each internal node. When these equations are 


combined with the semi-annular conformal finite element mesh 
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analysis, each equation, for a given interior node n, will 
relate eqQ(n) and hgQ(n) to 14 adjacent node azimuthal field 
values egq(i) and hqQ(i). Some of these adjacent field values 
may be EU on a boundary. The result of conducting 
these calculations for each value of n is a sparse global 
system matrix, containing complex elements and displaying a 
diagonal  tri-block  submatrix structure. Applying the 
Dirichlet boundary conditions to the semi-annular conformal 
finite-element mesh shown in Figure 2-2, the numerical 
solution for the electromagnetic fields in Region II is 
obtained by utilizing a two-sweep block-by-block matrix 
inversion algorithm related to the Riccati transform for 
difference equations. The diagonal tri-block submatrix 


structure is given by 


- (2.16) 
κο ο ον. 
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Ак, By, and Cy are sparse submatrices which constitute 
the tri-block diagonal structure of the global system 
matrix. Фр is generated from Equations 2.14 and 2.15 
through the use of the Dirichlet boundary conditions. 


[Ref. 2] 
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The mesh construction and matrix generation are 
conducted on a local basis in a sequential manner. The 
analysis starts at the positive Z-axis, and works clockwise 
in equal Eure Using a simple minimum distortion 
scheme at each step, the mesh is constructed to conform to 
multiple approximate piecewise smooth surfaces. The Riccati 
transform sparse matrix inversion technique conducts a 
"first sweep" and a "backsweep". On the "first sweep", the 


Values of Ry and Sy are calculated, where 


R = -( + “ 8 
К+] Уо C. боо 
^v ^v ^v ^v ^v 
z 
= E a а ) 
5% (δι + А, κι) (8. ALS)? Claes 
On the "backsweep", the values of the solution vector are 


obtained, where the solution vector is given by [Refs. 1 and 


2] 
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A principal advantage of this E method is that 
the Riccati transform technique provides a systematic 
procedure. This procedure divides a larger area into 
several erage areas defined by the tri-block matrix 
substructure, existing between any three adjacent radial 
rows of nodal values. Analysis using this type of smaller 
area procedure allows all operations including mesh 
construction, matrix element generation and matrix 


inversion to be performed at the same time. (Ref. 2] 


D UNIMOMENT METHOD 

The unimoment method is a computational procedure which 
allows numerical values to be calculated for the complex 
coefficients of the incident, scattered and spherical core 
electromagnetic field expansions when dealing with the semi- 
annular conformal finite-element mesh. The scattering 
solution is obtained by representing the electromagnetic 
fields of Region I and Region III, as shown in Figure 2-2, 
in terms of TE and TM spherical harmonic expansions, with 
unknown core and scattered field coefficients. Using 
superposition of boundary conditions and interior solutions 
for the azimuthal fields, the core expansion is established 
as a boundary condition along r = a while the scattered 
field expansion added to the incident field is established 
as a boundary condition along r = b. The two existing sets 
of boundary conditions are modal spherical harmonic core 
fields along r = a, with zero fields along r = b, and 
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spherical harmonic and incident fields at r = b, with null 
fields at r = a. The multiple boundary conditions will 
generate multiple solutions at the nodes where r = r) orr = 
rə, (i.e.,along the boundaries). The method used to obtain 
the expansion coefficients sets the weighted superposition 
of interior solutions to the appropriate original field 
expansions. Along the boundaries at r = rı and r = rs, this 
method enforces equality utilizing least squares analysis. 
[Refs. 1 and 2] 

This method of analysis is only valid if the scattered 
electromagnetic fields everywhere outside of the interior 
mesh region can be represented by field harmonics. These 
representations appear as infinite series, and must Бе 
Mathematically truncated to facilitate numerical 
calculations. Decomposing the plane-wave incident fields 


into azimuthal modes yields 
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where €, = 1 if m = 0, €m = 2 if m > 1; M is the truncation 
index, and the polarization convention (TE/TM) ,depicted in 
Figure 2-3, is used in the bracketed sinusoids. (Ref. 2] 
DECENT computations, the numerical scattering solution 
algorithm is run concurrently with the sequential Riccati 
transform finite-element solution in Region II. On the 
forward sweep, boundary conditions for both the incident and 
modal fields are generated at each node row being used in 
that particular sweep step. On the backsweep, boundary 
conditions along the гр = гру апа г = үу. contours are used to 
compute each nodal value. The incident field and spherical 
harmonic field for each node is formed, and used to generate 
the nodal residual fields. During the backsweep, inner 
product integrations are formed, and these inner product 
results are added to the field moment matrix, T, and the 
incident field driving vector, D corresponding to each 
incident field. The scattering solution is obtained by 
multiplying the inverse of the field-moment matrix with the 
driving vector to obtain the field expansion coefficients. 
The complex scattering amplitude ?(Ө,Ф) 15 defined in terms 


of the far-zone E-field, 
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where дах 18 the maximum scattering body dimension amd 


where the components of Ғ(Ө,%) аге given by 
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The scattering cross section is found to be [Refs. 1 апа 2) 


т | x: 2(n+1) ° (n+m)! 
σ------- Σε, Σ -------- 
. -m)! 
on | E! |? n=0 ΠΕ (n*1) (n-m 
n 0 
2 
ПИР ш 
Пуп My N 








Utilizing the theoretical principles and mathematical 
calculations just discussed, Morgan [Ref. 1] wrote the 
computer program EMSCAT to compute electromagnetic field 
scattering from layered bodies of revolution. Extensive 
interactive computer aided design features have been added 
to EMSCAT by Morgan and Connolly. These computer-aided 
design features, driving the change of name from EMSCAT to 
EMCAD, allow structures composed of multiple dielectric 
materials to be graphically constructed or modified by the 
user, and then evaluated by EMCAD for scattering. The 


remainder of this thesis is concerned with the utilization 
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and application of EMCAD, as well as with modifications and 


lmprovements of this transformation. 
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ἘΠΕ ΓΗ DATA CONSTRUCTION 


As has been previously explained, EMCAD computes 
electromagnetic field scattering from penetrable bodies of 
revolution using the coupled azimuthal potential formulation 
in conjunction with a variational finite-element technique 
and a tri-regional unimoment method. The numerical solution 
is performed using a two-sweep Riccati transform with 
circumferential marching. In order to provide relatively 
"fast", consistently accurate answers, yet be flexible 
enough to allow the user to vary any of a number of input 
parameters as conditions warrant, EMCAD requires sufficient 
input data to adequately depict the .penetrable body of 
revolution, as well as the desired incident and scattered 
field orientations of the user. This accuracy and 
flexibility is accomplished through the utilization of 
keyboard input and/or mouse input to form both material and 
structure input data files. Mouse input deals with the 
graphical interface, and as such will be discussed in 


Chapter ТУ. 


A. KEYBOARD INPUT 

Information obtained by ο input is fairly 
diversified in nature. This data ranges from input and 
output file names, to desired mesh parameters, from 
selection of phi angles, to selection of incidence angles. 
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To facilitate ease in understanding and = 5 the program, 
EMCAD allows a user to input the necessary information 
according to the level of expertise needed. For the 
"Novice" Ec the program gives brief explanations of 
input values, describes the required format, and gives 
typical values where applicable. For the other option, the 
"Expert" level assumes the user is familiar with the input 
parameters and formats, and simply prompts for the required 
inputs. A more detailed explanation of each facet of the 
keyboard input is contained in Appendix A, EMCAD USER'S 


MANUAL. 


B. DATA FILE INPUT 

Once the information obtained by keyboard input is 
established, the remaining data becomes that which describes 
the penetrable body of revolution being considered. This 
descriptive data can be separated into two distinct areas of 
Information: The first area of information is the data 
which conveys the material composition of each layer of the 
structure. Call this the "material file". The second area 
of information is the data which defines the structural 
shape of each layer of the penetrable body, referred to as 
the »"structur Filer". As shown in Figure 3-1, there are 
three methods of data construction which can be used to 
generate the "Material" file and the "Structur" file. The 
quotation marks indicate that the data files can be named at 
the user's discretion. Method 2 and Method 3 utilize 
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computer aided design features, and will therefore Бе 
discussed in Chapter IV, Interactive Graphics Features. 

Data construction by way of Method 1 is depicted in 
Figure 3-2 and is composed of special case programs which 
generate data files for specific geometrically shaped 
objects. For example, program SPHERE allows the user to 
define a layered penetrable sphere, offset from center if 
desired, with as many penetrable layers as can be handled by 
EMCAD at this time. The user enters material dielectric 
parameters, layer radii, and required resolution as 
prompted, and the program writes the necessary "Material" 
and "Structur" data into files іп the required  EMCAD 
format. CONE and CYLINDER were written specifically for 
EMCAD. As with SPHERE, they allow the user to establish the 
body of revolution as a multilayered structure. Source code 
for SPHERE, CONE, and CYLINDER can be found in Appendices J, 
K, and L respectively. Whenever a computer program 
generates data and data files, the capability exists for the 
data to be generated manually, or by any other means 
available, and then entered into the appropriate data files 
manually through the use of the keyboard. Keyboard entry is 
shown as a viable method of data construction in Method One, 


Two and Three. 
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Spheres, cones, and cylinders are simple geometric 
shapes, and are extremely useful for elementary 
applications; however, nature is complex, and not always 
made up of simple geometry and elementary shapes. To obtain 
realistic results, realistic penetrable bodies are needed. 
It would indeed be tedious to attempt to write a structure 
defining computer program for each new shape encountered. 
As such, the optimum solution indicates the need for the 


implementation of computer aided graphics. 
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ieee ERACTIVE GRAPHICS FEATURES 


Interactive graphics simply refers to the ability to 
take data at hand, and depict that data in some visual forn, 
where a graph is the most commonly utilized depiction. A 
special subset of interactive graphics, known as computer 
aided design, also allows the user to draw a picture on the 
screen, and convert that picture to data. 

EMCAD was written in the programming language of Ryan- 
McFarland FORTRAN (Appendix M). Microcompatibles GRAFMATIC 
(Appendix M), and later Jewell Technologies  GRAFPLUS 
(Appendix M) were chosen as the commercially available 
screen graphics and hard copy support for EMCAD output 
graphics. Because the output format that was used readily 
allowed conversion to an acceptable format, West Coast 
Consultants CURVE DIGITIZER (Appendix M) was chosen as the 


computer aided design package. 


A. COMPUTER AIDED DESIGN 

There are three means of constructing the data for 
penetrable bodies of revolution, as previously shown in 
Figure 3-1. Method 1 has already been discussed, and is the 
predecessor to Method 2. In Method 2, the "Material" data 
can be generated by one of the three programs written 
specifically for EMCAD, or by direct keyboard input into a 
data file if desired. The change from Method 1 occurs in 
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the construction of the structure defining data. As shown 
in Figure 4-1, the "Curveout" data is obtained by using 
CURVE DIGITIZER. In simple terms, this means that the user 
draws ene picture on the screen using a mouse with the 
computer aided design package. Once the user is satisfied 
with the single, or multilayered drawing on the screen, the 
data is saved in the format established by the computer 
aided design package. An example of the two-dimensional 
cross section of a three-dimensional body of revolution 
which was designed using CURVE DIGITIZER is shown in 
Figure 4-2. 

The data file that is saved is a valid description of 
the structure, but is not in a form which can be used by 
EMCAD. This format conversion is accomplished by a program 
written specifically for this purpose, EMCADIN. The 
computer aided design output data is written in X,Y format, 
with special codes interspersed throughout. EMCAD requires 
radius format and a normalizing factor of 27 must be taken 
into account. The output data is not necessarily evenly 
spaced, whereas EMCAD requires evenly spaced input data. It 
is conceivable that the data is of insufficient resolution 
for use by EMCAD. Primitive data should be of high 
resolution  (i.e., 361 points). The program  EMCADIN 
rectifies all of these problems, as well as computes the 


minimum and maximum radius values required. The output of 
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EMCADIN is the data file "Structur". A more detailed 
description can be found in Appendix D, EMCADIN USER'S 
MANUAL. 

Meena 3, as shown in Figure 4-3, is actually a union 
of Method 1 and Method 2. "Material" data is generated as 
described in both previous methods. Structural data is 
generated as in Method 1, but is then imported to the 
computer aided design package for modifications. The data 
to be modified must be changed by  CRVDIGIN, written 
especially for this purpose, into a format which can be 
accepted as input by CURVE DIGITIZER.  CRVDIGIN is explained 
in greater detail in Appendix F, CRVDIGIN USER'S MANUAL. 
Once the data has been read into the computer aided design 
package, the remainder of the development is essentially 


Method 2. 


B. INTERACTIVE OUTPUT GRAPHICS 

Having constructed the "Material" and the "Structur" 
input data files using one of the three methods discussed, 
it is now possible to conduct the scattering calculations. 
Figure 4-4 shows that the execution of EMCAD yields output 
data and output graphics. 

As also shown in Figure 4-4, EMCAD generates four data 
files which are used as input to the graphics program 
MCADGRAF. It is significant to note that the extensions of 
these filenames have special meaning. TEP refers to TE 
incidence (Figure 2-3),  F-Phi; TET refers to TE incidence, 
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F-Theta; TMP refers to TM incidence (Figure 2-3), F-Phi; TMT 
refers to TM incidence, F-Theta. 

MCADGRAF is a menu driven output graphics program. 
EE MCADGRAF USER'S MANUAL, gives examples of all 
the different input menus. The menu structure gives the 
user some general information concerning the use of the 
program, as well as user choice of aspect ratio, data file 
type selection, data file name input, incidence angle 
selection, and phi angle selection. Once these selections 
have been made, MCADGRAF produces four output graph types, 
giving the user the capability to make a printout of any or 
all of these. The first graph type is magnitude versus 
theta in a rectangular coordinate system, as shown in Figure 
4=5. The second graph type is phase versus theta in a 
rectangular coordinate system, as shown in Figure 4-6. The 
third graph type is magnitude versus theta in a polar 
coordinate system, shown in Figure 4-7. The final type is 
simply an axis rotated version of the third type, for a 


different perspective, as shown in Figure 4-8. 
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V. VALIDATION 


In зе to transform the existing program EMSCAT to 
the modified version entitled EMCAD, it was imperative to 
ensure that both the accuracy and the reliability of EMCAD 
were at least as good as that of EMSCAT. ТЕ was also 
essential to establish flawless compatibility between EMCAD 
and any program written to support EMCAD. 

The initial step used to validate the operation of 
EMCAD consisted of a comparison between the results obtained 
from EMCAD and the exact theoretical values of a Mie series 
far-field complex amplitude E-plane pattern for an offset 
lossless (€, = 2.6) penetrable dielectric sphere. Run 
parameters chosen included broadside incidence (incidence 


angle = 90 degrees), sphere radius equal to 0.8 (as a factor 


of wavelength), and z-axis offset = 0.16 (as a factor of 
wavelength). The series truncations at M=7, Nyl- and 
№м11=15 were chosen for the variable mesh scheme. these 


values were the same as those used by Morgan [Ref.1] in his 
initial validation of  EMSCAT. The results of this 
comparison are shown in Figure 5-1 for a phi angle of O 
degrees, and in Figure 5-2 for a phi angle of 90 degrees. 
Having demonstrated the close resemblance between the 
results obtained from EMCAD, and those predicted by theory, 


the next step used to validate EMCAD required a comparison 
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between the results obtained from EMCAD апа аслы measured 
data. Morgan (Ref. 1] validated the original version of 
EMSCAT with extensive experimental verification. Having the 
results of this experimental verification readily available 
presented the unique opportunity to validate the operation 
of EMCAD with the same experimental data used to validate 
EMSCAT. For each axisymmetric body constructed, parameters 
chosen included TE and TM incidence with incidence angles of 
0, 45, 90, 135 and 180 degrees. The series truncations at 
M=9, Nyt=15 and Nytt=20 were chosen for the variable mesh 
scheme. All experimental bodies were constructed of 
plexiglass (polymethylmethacrylate) with material parameters 
defined by €, = 2.59 - 10.017. (Refs. 1 and 2) 

The decision was made to demonstrate the versatility 
and reliability of the program by utilizing a number of 
Significant shapes which could be readily fabricated for 
experimental measurements. Since a sphere had already been 
used, and shown to work, a bicone was chosen as ап 
appropriate figure. À representation of the solid bicone 
utilized is depicted in Figure 5-3. The comparison for the 
solid bicone validation is shown in Figure 5-4. The hollow 
bicone of Figure 5-5 demonstrates increasing complexity. 
The sphere shown inside the bicone is a free space cavity in 
the object, although it could be made of any material. The 
comparison for the hollow bicone is presented in Figure 5-6. 


The next figure chosen was that of a hollow cylinder. Тһе 
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internal spherical cavity had the same ЕРИ as that 
found in the hollow bicone. The height of the cylinder was 
1.9 (in terms of wavelength) and the radius of the base of 
the SALONS was .945 (in terms of wavelength). The 
validation comparison for a phi angle of 0 degrees is 
depicted in Figure 5-7 and the validation comparison for a 
phi angle of 90 degrees is presented in Figure 5-8. 

To this point, all validation runs have correlated 
extremely well with the experimental data, and all 
experimental figures have consisted of combinations of 
Simple geometric shapes such as spheres, cones, апа 
cylinders. These are extremely useful for elementary 
applications, however, nature is complex, and not always 
made up of nice geometry or elementary shapes. To prove 
that EMCAD could handle more realistic bodies, the composite 
structure of Figure 5-9 was constructed. Since a program 
was not written to generate the "structur" data required for 
EMCAD input, the decision was made to utilize CURVE 
DIGITIZER (CAD package) to produce the required information. 
The representation of the composite body of Figure 5-9, 
required as input for EMCAD, is depicted in Figure 5-10. 
Note that this representation is composed of data defining a 
two dimensional axisymmetric body. The actual comparison 
for this structure is shown in Figure 5-11. As was the case 
before, the data generated by EMCAD correlates extremely 


well with the experimental data. 
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In addition to demonstrating the validity of the 
present version of EMCAD, the composite body of Figure 5-9 
helped demonstrate program compatibility. Since the output 
data of CURVE DIGITIZER was not in the format required for 
use by EMCAD, the computer program EMCADIN was developed and 
written to transform the output data of CURVE DIGITIZER into 
a format which could be used by EMCAD. Figures 5-12, 5-13 
апа 5-14 present figures which were generated by the 
programs SPHERE, CONE, and CYLINDER, transformed by CRVDIGIN 
into proper format, and imported into the computer aided 
design package CURVE DIGITIZER. The pictures shown are the 
graphical representations of the imported data as seen by 
CURVE DIGITIZER. This generation and transformation not 
only shows the validation of the programs written to 
generate specific geometric shapes, but also demonstrates 
the compatibility of these programs with CRVDIGIN, and the 


compatibility of CRVDIGIN with CURVE DIGITIZER. 
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This thesis has described a very useful numerical 
algorithm for computing resonance region electromagnetic 
interaction with arbitrarily shaped and even inhomogeneous 
axisymmetric penetrable objects. The background material 
associated with the development of this interactive, finite 
element scattering algorithm was presented in Chapter II. 
Included in this background information was the concept of 
coupled azimuthal potentials, the finite element algorithm 
and the unimoment method. Chapter III then presented 
procedures for constructing input data for the scattering 
program. Computer programs written specifically to generate 
a given geometry, and to allow data file input and keyboard 
input were examined in detail. The first half of Chapter IV 
discussed the extension of the methods used to construct 
program input data to those methods utilizing graphics. Of 
Special interest was the incorporation of computer aided 
design features. The second half of the chapter presented 
user interactive output graphics features. Validations of 
the modifications and improvements to EMCAD were presented 
in Chapter V. Extensive graphical comparisons were made 
between the output of the modified version of EMSCAT, and 
both theoretical calculations and actual measured 


experimental data. The impressive accuracy of the 
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comparisons validated the changes and modifications to 
EMSCAT зей іп developing EMCAD. Inter-program 
compatibility was also extensively demonstrated. The 
довава ааа provide detailed users' manuals for four of the 
major computer programs developed. 

This was only the first stage of a multistage project. 
Future developments of the axisymmetric code will proceed 
with major programming modifications to allow metallic 
surfaces to be embedded in the inhomogeneous penetrable 
material of the scattering object. Also, the near fields of 
the object will be made available for display. This is 
useful for identifying "hot spots" and field caustics. 

A second effort has already begun to develop a fully 
three-dimensional scattering object capability. The basis 
of this development was the discovery of a generalized 
coupled azimuthal potential formulation which does not 
require an axisymmetric material configuration [Ref. 8]. 
Such a computer algorithm will ultimately employ three- 
dimensional graphics in a powerful computer aided design 
package for user design of scattering objects as well as the 


output of scattering patterns and near field quantities. 
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APPENDIX À 


EMCAD USER'S MANUAL 


Welcome Ко  EMCAD. This program  computes 
electromagnetic field scattering from layered bodies of 
revolution using the coupled azimuthal potential formulation 
in conjunction with a variational finite element technique 
and a tri-regional unimoment method. The numerical solution 
is performed using a two-sweep Riccati transform with 


circumferential marching. 


A. BACKGROUND 

EMCAD is based upon an earlier code named EMSCAT which 
was written during the period 1974-1976 by Dr. M. A. Morgan 
while at the University of California at Berkeley. 
Extensive interactive computer aided design (CAD) type 
features have been added in this effort. These CAD features 
allow structures composed of multiple dielectric materials 
to be graphically constructed or modified by the user, and 
then evaluated by EMCAD for scattering.  EMCAD has been run 
in both Ryan-McFarland FORTRAN (Appendix M) and Microsoft 
FORTRAN (Appendix M), with the current version using 
Microsoft. Microcompatibles GRAFMATIC (Appendix M) and 
later Jewell Technologies GRAFPLUS (Appendix M) were chosen 
as the commercially available screen graphics and hard copy 
Support for EMCAD output graphics. Because its output 
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format allowed easy conversion to an acceptable format, West 
Coast Consultants CURVE DIGITIZER (Appendix M) was chosen as 
the computer aided design package. Any questions concerning 
colors, ΛΙ changes to the program, etc... should be 
referred to the respective user's manuals, Dr. M. A. Morgan 


or LT. E. M. Connolly (Appendix M). 


B. COMPUTER SETUP 
EMCAD was written for the following computer 
configuration: 
* IBM AT or compatible, including 80386 based 
compatibles 
* hard disk 
* EGA graphics card 
* EGA monitor 
* math coprocessor - not required, but 
highly recommended 
х memory ~ 640K RAM minimum 
EMCAD can be placed in any directory, at the user's 
discretion, however, if a RAM-drive is to be used, is highly 
recommended that during each run the user allocate a minimum 
of 4M of RAM for EMCAD and EMCAD input/output files.  EMCAD 
input is taken from a variety of sources. The data files 
from these sources are required to be in the same 


directory/sub-directory as EMCAD. 


66 


C. PROGRAM INPUT 

EMCAD requires both keyboard input and data file input. 
As shown in Figure A-1, two input data files, "MATERIAL".DAT 
and "STRUCTURE". DAT, are required as input to EMCAD. The 


quotation marks indicate that the data files can be named at 


the user's discretion. These two data files contain 
information which describes the penetrable body of 
revolution being considered. The "material file" contains 


the data which conveys the material composition of each 
layer of the structure. The "structure file" contains the 
data which defines the structural shape of each layer of the 
penetrable body. The three methods which are used to 
generate these data files are discussed in much greater 


detail in both Chapter II and Chapter III. 


D. PROGRAM EXECUTION 
EMCAD is a menu driven program. Once the correct 


directory has been selected, EMCAD can be invoked by typing 
EMCAD 


at the DOS prompt. This is shown in the top line of Figure 
A-2. The program will load itself, and the user will see 
the next two statements shown 42 A-2 requesting the 
names of the "material file" and the "structure file". 
Sample keyboard inputs are shown for these two requests. It 
must be pointed out that the extensions are included (.dat) 
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Û ‘\THESTS\MSFORTIEMCAD 


ENTER MATERIAL PARAMETER FILE (D:FILENAME. EXTENSION) : 
ndatl.dat 


ENTER STRUCTURE DATA FILE (D:FILENAME. EXTENSION) : 
sdat1.dat 


Offset Sphere h-8.8 D-8.16 
Erz(2.6,8.) Urz(1.,8.) 


PLEASE PRESS ANY KEY TO CONTINUE. 


ieee Ne UP SATA FILES 
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as required by EMCAD. The two lines beginning with "offset 
Sphere..." refer to two label lines which were previously 
inserted into the "material file" to allow the user to 
differentiate this set of data from every other set of data. 
By choosing labels carefully, the user can discern from the 
very first screen whether the proper set of data is being 
loaded. This becomes very important because as figures 
increase in complexity, run time for the program increases. 
It is not uncommon to have runs which take several hours on 
an AT class conputer. It would be very disturbing to run 
the program for seven or eight hours only to find out that 
the incorrect set of data had been used. Upon completion of 
this screen, the user depresses ENTER to see the following 


screens: 


* Welcome To EMCAD (Figure A-3) 


This screen presents some brief background information. 


* EMCAD Level of Expertise (Figure A-4) 


To facilitate ease in understanding and running the progran, 
EMCAD allows a user to input the necessary information 
according to the level of expertise needed. For the 
"Novice" level, the program gives brief explanations of 
input values, describes the required format, and gives 


typical values where applicable. The "Expert" level 
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ще 


assumes the user is familiar with the input parameters and 
formats, and simply prompts for the required inputs. The 
remaining menus will use the "Novice" level to cover all 


ranges of user. 


* Input Name for Output Data File 
(Figure A-5) 
This menu prompts the user for the filename, WITHOUT 
EXTENSION, that will be used to hold the output data of 
EMCAD. The program will append the extension .OUT to the 


user's filename. 


* Input Name for Output Graphics Data Files 
(Figure A-6) 
EMCAD produces four output data files for input to the 
graphing routines of MCADGRAF. А11 four data files have the 
same name, but each has a distinct extension related to the 
type of data the file contains. This menu prompts the user 
for the name of the output graphics data files, WITHOUT 
EXTENSION. This filename may be the same as that used for 


the output data. 


* Graphics Caption (Figure A-7) 


This menu affords the user the opportunity to input a 


graphics caption to help identify this set of data from all 
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THE OUTPUT DATA FILE IS THE FILE CONTAINING THE OUTPUT 
RESULTS OF ALL NUMERICAL CALCULATIONS CONDUCTED BY EMCAD. 
THE FORMAT FOR THIS INPUT IS FILENAME ONLY. NO EXTENSION 15 
REQUIRED OR DESIRED. EMCAD AUTOMATICALLY APPENDS AN 
EXTENSION OF .QUT 10 YOUR FILENAME. 


PLEASE ENTER THE FILENAME OF THE OUTPUT DATA FILE. 


FIG. A-9- BT NN NIS 
ОИТ DoT ЕЕ 
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THE OUTPUT GRAPHICS-DATA FILE IS THE FILE CONTAINING THE 
OUTPUT DATA FROM EMCAD TO BE USED AS INPUT TO GRAPHING 
ROUTINES. THE FORMAT FOR THIS INPUT IS FILENAME ONLY. МО 
EXTENSION 1S REQUIRED OR DESIRED. | EMCAD AUTOMATICALLY 
APPENDS AN EXTENSION OF TMT, . TMP, .TET, AND .TEP TO YOUR 
FILENAME AS IT PRODUCES FOUR OUTPUT FILES FOR GRAPHICS. 


„IMT --—> TM INCIDENCE, F-THETA 
‚ТИР --——) TM INCIDENCE, F-PH] 
„TET --—> TE INCIDENCE, F-THETA 
TET ---—) TE INCIDENCE, F-PH] 


PLEASE ENTER THE FILE NAME OF THE OUTPUT DATA FILE. 


ДВЕ А0: ТР ра а МАМЕ OR 
ш е κ [S DATA FILES 


ЛӘ 


THE GRAPHICS CAPTION IS A PERSONALIZED CAPTION ALLOWING THE 
USER TQ IDENTIFY THIS SET OF GRAPHS FROM ALL OTHER SETS. 
THE MAXIMUM LENGTH OF THIS CAPTION IS 64 CHARACTERS. 

NOTE: WHEN USED WITH THE GRAPHICS PACKAGE, THE PROGRAM IS 
ABLE 10 DIFFERENTIATE BETWEEN UPPER CASE AND LOWER CASE 
CHARACTERS. 


PLEASE ENTER ANY GRAPHICS CAPTION YOU DESIRE 


FIG. Απο ORA πο ο ο 
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other sets. The graphics caption appears on the first line 


of the header of each graph produced by MCADGRAF. 


* Mesh Density (Flgure A-8) 


As discussed in Chapter II, EMCAD utilizes a semi-annular 
conformal finite-element mesh in order to conduct scattering 
calculations. This menu prompts the user for the minimum 
and maximum mesh density parameters. At azimuthal modes M = 
О апа 1, Όπαχ is used as the density parameter. The density 
decreases linearly to Dnin at the final mode step. Notice 


that these values are REAL. 


* Incidence Angles (Figure A-9) 


EMCAD allows up to five incident fields to impinge on the 
object of interest. This menu asks for the number of 


incidence angles, and then the value for each in degrees. 


* Number of Scattering Field Theta Points 
(Figure A-10) 
EMCAD conducts scattering calculations by starting at Ө - 0, 
and progressing clockwise to @ = 180 in equal increments. 
This menu allows the user to determine the size of the 


increments, and indirectly the accuracy. 
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DMIN AND DMAX ARE PARAMETERS OF MESH DENSITY IN TERMS OF 
ELEMENTS INTERIOR LAMBDA. INPUT VALUES ARE EXPECTED TO ВЕ 
REAL, I.E. THE DECIMAL POINT MUST DE INCLUDED. 

TYPICAL VALUES ARE DMIN = 10. AND DMAX = 15. 

PLEASE INPUT DMIN 


16, 


PLEASE INPUT DMAX 


15, 


FIG. A-8. ME ο. 


ΤΕ 


THE NUMBER OF INCIDENT FIELD ANGLES [S THE TOTAL NUMBER OF 

INCIDENT FIELDS THAT IMPINGE ON THE OBJECT OF INTEREST. 

THIS PROGRAM ALLOWS A MAXIMUM OF FIVE INCIDENT FIELD 

ANGLES. WHEN ENTERING YOUR ANSWER PLEASE DO NOT INCLUDE À 

DECIMAL BECAUSE THE INPUT MUST BE IN INTEGER FORMAT, I.E. 
=3 , MÑ=1 

PLEASE INPUT THE NUMBER OF INCIDENT FIELD ANGLES. 


2 


ENTER INC FLD ANGLE (DEG) 1 


8. 


ENTER INC FLD ANGLE (DEG) FOR 7 


90. 


Ur Е (ea Nee ВЕЕ = 


E 


THE NUMBER OF SCATTERING FIELD THETA POINTS DETERMINES THE 
SPACING BETWEEN THETA POINTS DURING EMSCAT ITERATIONS AND 
CALCULATIONS. 

DELTA THETA = 188 / (NUMBER THETA POINTS - 1) 50... 
NUMBER THETA POINTS = 37 -—-> DELTA THETA = 5 DEGREES 
NUMBER THETA POINTS = 19 -—--> DELTA THETA = 10 DEGREES 
WHEN ENTERING YOUR ANSWER, PLEASE DO NOT INCLUDE A DECIMAL 
BECAUSE THE INPUT MUST BE IN INTEGER FORMAT. I.E. NT = 19 


PLEASE INPUT THE NUMBER OF SCATTERING FIELD THETA POINTS 


FIG. gg A IG. ΑΟ пин 
SCAT TERING F IE pR И И ee 


во 


* Phi Angles (Figure A-11) 


EMCAD allows up to eight phi angles which define the 
location from which the user views the scattering for all of 
the incident angles. This menu asks for the number of phi 


angles, and then the value for each in degrees. 


* Expansion Modes (Figure A-12) 


EMCAD partitions the incident and scattered fields into 
azimuthal modes, indexed by "M'. The solutions for these 
modes are accumulated during the "M-loop". If 0 degrees and 
180 degrees are the only incidence angles of interest, then 
the M = 1 mode is the only mode that is needed. The values 
of N1 and N2 represent the number of spherical harmonic 


modes in the spherical core and the external region. 


E. PROGRAM OUTPUT 

As previously mentioned, EMCAD produces text output data 
and output data for graphics. The graphics data is 
discussed in Appendix H, MCADGRAF User's Manual, and the 
text output data is demonstrated in Appendix B, Portion of 


an EMCAD Sample Output. 
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‘THE NUMBER OF PHI ANGLES IS THE TOTAL NUMBER OF PHI ANGLES. 

THIS PROGRAM ALLOWS A MAXIMUM OF EIGHT PHI ANGLES. WHEN 

ENTERING YOUR ANSWER PLEASE DO NOT INCLUDE A DECIMAL 

BECAUSE THE INPUT MUST BE IN INTEGER FORMAT. I.E. 
ПИТИ 

PLEASE INPUT THE NUMBER OF PHI ANGLES. 


2 


ENTER PHI FLD ANGLE (DEG) FOR 1 


98. 


ENTER PHI FLD ANGLE (DEG) FOR ὁ 


198. 


FIG. Απ μαι λος» 


за 


ENTER MSTOP (.LE. 13) 
ESTIMATED "MINIMUM" VALUE IS: 1 


1 


ENTERING NO. OF INTERNAL AND EXTERNAL EXPANSION MODES 
ESTIMATED (КІМІМ) “MINIMUM” INTERNAL NL IS: 12 
ENTER N1 (.GE. MSTOP) 


12 


ESTIMATED (KO*RMIN) “MINIMUM” INTERNAL N2 IS: 21 
ENTER NZ (.GE. MSTOP) 


21 


ΕΞ. Е τς Nol ON MODES 
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APPENDIX B 
PORTION OF AN EMCAD SAMPLE OUTPUT 


This Appendix demonstrates a portion of an EMCAD 
sample output for the specific case chosen where 
M indexes from 0 to 7 inclusive. Only the M = 1 
case is fully shown, as well as the accumulated 


solution for all eight modes. 


XXXXXXXXXXXXXXXXXXXX EMCAD OUTPUT DATA ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
COMPLEX Er(n) AND Ur(n) 
1 € 2.610Е+00 EL.000ES06) (С 1.000Е+00 -1.000E-06) 


INCIDENT FIELD ANGLES 


N ALPHA(N) 

1 0. DEG 

2 IU, DEG 

5 180. DEG 
SERIES SOLUTION PARAMETERS 
MSTART= 0 
MSTOP= 7 
N1 = 8 

оз 13 
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XXXXXXXXXXXXXx M - ] XX € € € € € € € € XK XK X 


MESH PARAMETERS 
RADII (WAVELENGTHS) 


RADMIN= -64000E+00 
RADMAX= .96000Е+00 
RADA= .55241Е+00 
RADB= .10676Е%01 


INTEGER CONSTANTS 


NDS= 15 
IMX= 59 
NR1= 1 
NR2= 1 


MESH DENSITY AND ELEMENT SIDE RATIOS 
DENSITY/LAMBDA-OUT = 24.233 
DENSITY/LAMBDA-IN = 15.000 

DSA/DR = СВ 

DS5ZDR = 1.565 
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ВІССАТІ MARCHING STE 
.60SETUL 
250 3E CUM 
2502 E: E01 
СОг Е ЕО 
.600Е+01 
.599E+01 
.597Е+01 
.595Е+01 
„ БЗ2ЕС 
.589Е%01 
‚ 586Е+01 
.5835ET01 
.579E+01 
ο σετ 
STILE DL 
.566E*01 
5625701 
.557Е401 
„ВБЗЕЕЦИ 
.547Е+01 
БЕСЕТІМ 
.537E+01 
TSSL EROT 
.5266+01 
25d UE EDI 
.515Ef0/0] 
.509Е+01 
.505E+01 
.G98E+01 
.G692E*01 
. 487 E01 
.482E+01 
.477Е+01 
.471Е%01 
.066Е%01 
.461Е+01 
.857E+01 
.452E*01 
.498Е+01 
.4%3Е+01 
.439Е+01 
„аз 5БЕБОШ 
za ScCETOT 
2 23S EON 
„425Е+01 
„422Е+01 
.419Е+01 
.416Е+01 
алас 
,412Е4ғ01 
„410Е+01 
.408Е+01 
.406Е+01 
.405Е+01 
.404Е+01 
.403Е+01 
.403Е+01 
.402Е+01 
.402Е+01 
.402E*01 


Гега бо во ОКО 00 4 Os Un 4S CA P) — t5‏ ی 


ЕКСЕ 
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ШЕ АШ ЭСАТТЕКЕШ ТЕК EXPANSION COEFFICIENTS 


AZIMUTHAL MODE NUMBER M- 
ANGLE ALPHA 


ο ο ΟΕ 


ОШАСБЕКЕР" СОЕГЕРТСТЕНТО: 


Oo MI ON tn FD A NF Z 


ΠΕΤΕΡ COEFFICIENTS: TE INCIDENCE 


COONAN DN HZ 


STAGGERED COEFFICIENTS: TM INCIDENCE 


U (NAN — O <O O9 — ON Un RNN Z= 


= m а а ә [-4 


STAGGERED COEFFICIENTS: 


Aan ANNE Z 


1 


0. 


DEGREES 


XXXXXXX INTERIOR COEFFICIENTS XXXxxxx 


SE DD ESO 
. 209۴00 
. 105E ۴00 
.207E+00 


.266Е%00 


TOE On 


ae Er 00 
πο ες] 
eHow es Ol 
239 0E T 0T 
.158E=01 


3296-07 


ΤΗ 
ACN) 


ACN) 


INCIDENCE 


.266E+00 
ιδ οἱ 
= I7 τοῦ 
που ει 
OUI 
"oo οἱ 
198E 01 
SESE ος 


.S00E 00 
.209E+00 
zT SEED 
S209 (00 
1995 E 302 
4906 E-UI 
.616E=01 
ПО ea On 


1. 


,318Еғ00 
"239E t00 
xLAPETDD 
21655-01 
ΕΕ 
РІП ЕД ШІ) 
pou oes 01 
625E 02 


.298E+00 
MOZ ETON 
.66 06-01 


-.167E+00 


ЭЕ КЕЕ 02 
τοι 
2S SEU 
.136 E201 


XXXXXXX SCATTERED COEFFICIENTS XXxxxxx 


ο, OL 
πι εκ οι 
SOS 010 
TSS ETOO 
.472E-01 
ШӨСЕСПІ 
ΞΘ ΕΠΣ 
.066Е-03 
.572Е-06 
. 192E 05 
. 15106 
606E 0, 
SIE ος 
σος ὃς 
-336 E=03 


.144Е-01 
.195E+00 
ПИБ E00 
T 31E+00 
.956E-01 


Е ТБЕЕОТ 


CON) 


CCN) 


-.164Е-01 
PISSETO0 
„157Е+00 
.131€+00 
.956E-01 
.616E-01 
КЕЕ ОИ 
12062 = 02 
ӨШСЕ - 015 
,497Е-04 
jog Ze=05 
IIE 05 
.207E-06 
ΕΠΙ; 
.776Е-08 


TE INCIDENCE 


οτε οι 
ИП БЕЗ 
ІШЕ) 
z LSOETOD 
“нана Ез ОЛ 
ος σοι 
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.104E+00 
955E 02 
Nd EST 
Seon. Ol 
35 oe = 0'1 
Ses 0 1 
πο ε 01 
Pogo mic 
2520 Е- 03 
.451Е-04 
'6356E=05 
1116-05 
255E 06 
SSE 07 
СЕЕ 07 


paces 00 
DESSERT) 
.S50E=01 
1696Е= 01 
απ ετοῦ 
soo Ε-δι 


Β(Ν) 


Β(Ν) 


129686000 
SLOZETUD 
.660E-01 
147 ЕЗ00 
991 E=02 
SJJ ESO] 
GSE UL 
νο. αι 


„ЗПЗЕЪО 


- ΕΙ 


DCN) 


.172E+00 
[ΕΠΙ 
SISSE 
.109E+00 
ΕΠΙ 
16206502 


. 22%Е+00 
252 ЕУ 00 


= oo Ua oll 


DON) 


ΕΙ 
.158E+00 
"6996701 
ΕΓ 
16 0Е= 02 
TOSE- 05 
Ες 
ΕΕ 
Ε-06 
l1 EUG 
woe ees Ol 
ὅσος- 06 


.104Е+00 
„ОББЕ=02 
.996E-01 
,928Е-01 
ο οευι 
.510E-01 


— — — — — — 
Un х оч IN) I со мо 09 344 
birt d d 04d» !' 


„17 1E OR 
cs SUE 
СІЗ Еа 5 
.897E-06 
26 9 ID 
LILIES 
-LOrE 06 
839507 
TSE Uo 


NEU 
«ΠΟΤΕ 
sof ce=04 
ICE U 
| 15d UG 
οὐ, πο, 
0E 0, 
9 Е 200 
59:56 E05 
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„135Е- 00 
КӨШЕСІ 
20 St 0S 
EES 
SIE ТӨ 
2422Е-06 
ΕΞ» 
2992 E207 
ο Ες 


кү гг ΠΠ 


.159 EUN 
: 293E ШЕ 
„эске 
.451Е-09 
GS EEUD 
.11ТЕ 8 
.e9 55S 
.2 2E SM 
«112Е ІШ 


PNT ER DOR) SCATTERED FIELD EXPANSION COEFFICIENTS 
MODE NUMBER M- 


AZIMUTHAL 


IHCIDEHCE ANGLE ALPHA 


Oo — ON n T" CN [Ñ (— = 


STAGGERED COEFFICIENTS: 


09 М O" On T" CN [SÑ i— Z 


Un +T." C4 ^9) i— О «О O° -.J ON Qn £" CN (S i— Z= 


hh а 


STAGGERED COEFFICIENTS: 


ON n +" CA P9 + Z 


το 
OIE 02 
.354E=01 
па6 0Е-05 
‘650E 02 
ПОСЕН 
1986502 
< 08 E703 


.6602E=-01 


.294E+00 


SSE 0I 


.445E~01 
603E O1 


LECE 01 
ο 1δ]ε 01 
-,57GE 02 


I ! 


N N N — CI m ee 


I I DNR ON DAO 
NOON) Orns O 
ENN NN OF 
rn rn m m m rm rm IT 
! 
оооооо о о 


ЕЕ 08 


SSE UE 


22605201 


poo cB 


ZOE 0 
oU PEST 
763 E202 


ACN) 


ACN) 


CEND 


CON) 


-- 


1 
207 


TM INCIDENCE 


ze bor 10 
.874Е-01 
sooo” | 
Poo ESUZ 
ΠΕΙ 
of SLES 
290E UE 
.106E=02 


TE INCIDENCE 


. 36GE-01 
953E Ul 
ОЧЕ ОТ 
.5%6Е-01 
231E 01 
ТӨЛЕШ 


-904E 2072 


. 1955202 


TM INCIDENCE 


ΚΕΕΠ 
:e0UESUN 
ЭТЕСЕ 
2606 E09 
26602 
.408E-02 


16ος 03 
Ig ES 


19E 05 
.169E-06 
520ЕС06 
Da SES 
219 СЕ ШИ 
.301ËE 08 
CCE US 


TE INCIDENCE 


291E 0 
267E 01 
.47 501 
.668E-01 
. 544E-01 
. 54E 02 
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DEGREES 


ШІЛ |І 


жжжжжжж INTERIOR COEFFICIENTS XXXXXXx 
ERAGOERED COEFEICIENTS: 


του 02 
. 246 E+00 
. 262E οἱ 
: 980E 01 
655E O1 
.206E-02 
61E UZ 
SLE UE 


.624Е%00 
οτε ο] 
56 2EËE=01 
ПЕ ат 
NITE UL 
SIGL ESUS 
.296Е-05 
206 Е-05 


XXXXXXX SCATTERED COEFFICIENTS ΧΧΧΧΧΧΧ 
STAGGERED COEFFICIENTS: 


r390 E01 
ЗОВЕ НО 
.600E-01 
moa. Ea Ol 
“902Е-01 
. 695E 02 
ISE 0L 
ПОЕТО 
.402Е-09 
2346Е-09 
.603E-06 
. 97E 07 
post 
56026-7095 
τους» 


EE UE 
.152E=01 
ЗеСТ 
Ἱυδα-υι 
ASTE OL 
zoo ABSUZ 


BCN) 


BCN) 


DCN) 


DCN) 


CUE UL 
152E F00 
626 Е-ОІ 
νυ ΕΠΙ 
G7 SE 01 
SEO 
«θεος 0 
. 2E OE 


267E 01 
SE Ul 


LOL EU 
-169E F00 
1696 Е=0 2 
EDO 01 
.368E-01 
SEE U 2 
POE ae 
USE 02 
ШӘЛГЕЗ 
„265Е- 04 
POE 
2596-06 
Та8ЗЕ-07 
552 bE*)9 
OE 209 


5%І92Е-05 
,294Е-01 
cr E OI 
619E 02 
SOLE 
ας 


— — — — — — 
Un 4^ CA P9 — O «o 09 54 


2h62E-92 
eo Oba 
;,5.5 Bi = 06 
2012205 
. S005٥ 
(20 9 EUM 
25E ος 
JR 5T EU 


-  З60Е-ШО 


ώρας UL 
.848Е-05 
-133E F05 
.186Е-04 
-1.ΠΕΕΞΗΡ 
КОБЕ Uy 
σεξ ος 
т 2062596 
3 1576502 


90 


„1635 Ема 
ма а е 
,187Е-04 
29206-05 
2550/6 
пее Е- Ше 
сдава ши 
SE 0o 
=, 153E-038 


.57 CES 
.lISESUE 
‚с6 Т ыш 
.e9 6 E05 
.G SP ESUD 
.700Е 
. 11 ШІ 
, 176E 20S 
-. 520 ERNE 


ACCUMULATED FORHARD AND BACK SCATTERED AMPLITUDES 
M= l 


STAGGERED AMPS: INCIDENT TM 


ALPHA 

CDEG) FWD BK 

Ü. -.9ШБЕШШІ ЕВИЯ Sloe 0 | E5SUEFTUI 
0 =.c9SE noe = cD eo Peuve tt) | molet OL 
S0. =. 7IOE OL SL ES ο Ετος .474E+01 


STAGGERED AMPS: INCIDENT TE 


ALPHA 
( DEO) FWD BK 
Ü. СВОЕ ОЯ - IN Eq ӘЗ ФЕН ΕΙ 

90. τε ΕΤ] Saco Se SU .196E=01 wee Ges 00 
180. =. L OER] τ. ZETO? ο πατε .97 4Е+01 
SCATTERING CROSS SECTION CONVERGENCE 

M= 1 

NINC = 1 

N NEW TERM PCT CHANGE 
l ες 100.0000 

2 СЕЗІ) 80.6270 

5 .7EF00 69... 592 

4 ТЕТО 58.7650 

5 Peer Ou 47.0760 

6 ЕТО 175.5939 

7 ЕТО 1. 7682 

δ πε πε 207-50 

9 796-04 “0017 
10 eE 05 .0000 
11 GE O ‚0000 
ше „ВЕТО .0000 
Що πες .0000 
19 ОЕ ‚0000 
15 OE . 0000 
NINC = 

N NEM TERM PCT CHANGE 
1 ΘΕῈ 100.0000 

2 „дЕЗОШ (oa 

Э 26Е- 0 18.8104 

4 ЭЕ 52.9934 

5 1600 17.4571 

6 .1E—-01 1.6261 

7 ЗЕ-ШЕ . 5695 

8 2E US 20501 

2 SE US 0012 
10 E 06 .0000 
11 .9Е-08 .0000 
12 «τες 10 .0000 
1» SES . 00060 
15 ESEZ . 0000 
15 5 ЗЕ-16 .0000 
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4 


NEW TERM 
.3E-0] 
. 3E*00 
JE DOT 
3E UL 
TE Ol 
ZE O0 
τς 01 
. 58-05 
CBE 
.GE- 06 
TE NS 
.4Е-09 
ЕО 
Е 
25Е-15 


= 
= 
< 
O 
и 


5 


-- 
ри 


Un .P* CA P9 ¬ © М Оо J O" Un P Ga P9 س‎ 


NINC = 6 

NEW TERM 
.6E=01 
.2E+00 
„(ЕЛО 
LEROI 
ΕΙ 
ΕΠΙ 
ТЕТОВО 
ЗЕЕ 0С 
.9E-04 
126505 
AE SS 
„2e 08 
ТЕЗІ 
δει 
δει 


Камара ко а а 
un fN ¿A N) — сә 4 O0 —J CN Un TN Ca N I 2” 


m 

== 

Ω 
M 


NEW TERM 


οι DNN — O .O O92 + G" Un P CAP) I ZZ 
N 
m 
! 
eo 
N 


н M) de ge — 


за 


СНАНСЕ 


‚0000 
92223 
. 452 
. 6074 


CHANGE 


. 0000 
.6270 
газ Е 
.7450 
.0760 
„5856 
2664 
07-30 
„Оа 
‚0000 
‚0000 
«0000 
«9000 
«0000 
‚0000 


CHANGE 


.0000 
5576 
‚6869 
„2985 
.8194 
20565 
‚4062 
пас 
‚ОП 
‚0000 
‚0000 
‚0000 
‚0000 
‚0000 
‚0000 


MODAL CRUSS SECTION EVOLUTION 
NORMALIZED W.R.T. WAVELENGTH SQUARED 


INCIDENCE ANGLE = υ DEGREES 
INCIDENT TM POLARIZATION INCIDENT TE POLARIZATION 

M SCT ABS TOT SCT ABS TOT 

0 . 000Е+00 . 000Е+00 .000E+00 .000E+00 .000E+00 .000E+00 
T 395 SE $0] 205 EU) soo Det 0) «ΡΟΣΕΤΘΙ 105E 01 ΕΠΙ 
2 ο Σετ SUSE >01 ΕΠΙ .558E+01 7:0 5E 0 Ὁ πο] 
3 ا ود‎ USE 0L ο ا‎ EOS GER ФАП UL 365510] 
ч КО5бЕ+О1 ОБЕ 01 ο 0. .D558E+01 7 UOSE a SOS SEU 
5 .558E*01 ο ο JD 5E TU ЭСЕП -7 05E301 ΘΕΊΩΙ 
ό Ο2δΕΤΟΙ OSE 0l Doo peo 255 SEIT «7096-01 πε ΕΠ] 
7 5810] ο ε-οι ЭСЕ РОТ 2998850] θε οι .505F101 
INCIDENCE ANGLE = 90. DEGREES 

INCIDENT TM POLARIZATION INCIDENT TE POLARIZATION 

M SCT ABS TOT SCT ABS τοῦ 

0 .247Е%00 .100Е%00 . 3567 E* 00 ΕΙ -1056+00 ЕД Е НО 
1 ΠΕΠ "9936-01 ΠΕΠ .834E+00 .105Er00 .960Е%00 
2 DOT ET ος οὐ СОВЕТОТ pai DL .110E+00 пат ЗЕ 
5 ТТЕРІ .106Е%00 Soccer 12855001 .106Е%00 ο ος ο] 
ú .440Е+01 .115E+00 .451E+01 .396Ef01 .118Е+00 .807Е%01 
5 свалено .1449Е+00 .536E+01 ӘПЕРШІ ΕΠ БӘЛЕНІ 
6 .525Ër01 „166Е#00 s56GcETU] ПОВЕО 2166000 „БЗБЕКОШ 
? ο. Ε.Ο .166E+00 ӘЗ ЕТІ ΕΠΙ sci oe 00 .5q5E+01 
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INCIDENCE ANGLE = ISG ευ κε» 


INCIDENT TM POLARIZATION INCIDENT TE POLARIZATION 


94 


Scy ABS TOT s ABS TOT 
.220Е-12 -.220Е-12 .000Е%00 .209E-12 .209E-12 . 000E+( 
546E+01 .255Е-01 .548E+01 546E+01 .255E-01 .548E+( 
.546Е+01 .255E-01 .548Е%01 .546Е+01 .255E-01 .548E*( 
546E+01 .255E-01 .548Е%01 .546E+01 .255E-01 ‚548Е+0 
546E+01 .255E-01 .548Е%01 546E+01 .255E-01 .548E«( 
.546E*01 .255E-01 .598Е+01 .566E*01 .255E-01 .568E+0 

| 
.546E*01 .255Е-01 .548Е%01 .546E+01 .255Е-01 548Е+0 
.546Е+01 .255Е-01 .548Е%01 546 Е+01 255 Е-01 548Е+0 









eS? ڪڪ‎ 





ΧΧΧΧΧ SCATTERED FIELD COMPLEX AMPLITUDES xxxxx 


INCIDENCE ANGLE - 0. DEGREES 
PHI- 90. DEGREES 


STAGGERED RAD FLDS: TM INCIDENCE 


TRETA F-THETA 

MAG DEG 

0. 651-06 .00 
SE "S 21E 06 .00 
10. 7506 . 00 
15. .608E-06 .00 
2U. .961Е-06 . 00 
25. Ко 2сЕ=06 B 
50. 1229606 200 
за: 2207 E-06 .00 
60.  225Е- 06 ‚00 
65. "259E=06 . 00 
SU. ESSE 06 .00 
955 па ТЕЗ Ше . 00 
60. ПОТ Е- 06 200 
65. SSE 06 .00 
70. ПИ БЕ -00 ‚00 
755 "16 ПЕ—06 .00 
80. οἱ σε πό 200 
δ... .111Е-06 20:0 
o BoE O07, ‚00 
95. E2995E-07 . 00 
100. Ш20Е-06 .00 
1057 "l 57E=06 .00 
τα 15E 06 . 00 
TL5? .1l21F 06 . 00 
120- EE SgIESUT .00 
ου E95 AES07 ‚00 
150 155 5Е-07 «00 
155 Poe = Om, ‚00 
140 П 29Е=06 ‚00 
145 (55Е-06 200 
150 » Wes SSK . 00 
155 1 7 ЗЕ -0е .00 
160 .186E-06 20:0 
165 7202Е-06 ‚00 
170 "Z 28E=06 . 00 
1755 .241Е-06 .00 
160. .268E-06 .00 
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MAG 


. 195+02 
„187Е+02 
ΙΕ; 
πμ τος 
POZE OU 
.999Е%01 
.116E+01 
.253E+01 
.886Е%01 
.603E+01 
.608E+01 
.515E+0] 
КӘЗЗЕТІШІ 
.159E+01 
PSSZ ETOO 
.187E+01 
"nog ep 
. 307 E* 01 
„21О0ОЕ+П1 
.186E*01 
δι Ες 
.290Е%00 
. 996 Е+ 00 
71:5 3 E350 
"I2l1ETÜ0I 
.716E-*00 
.189E+00 
.75%6E+00 
mie set OL 
ПЕРЕН ОТ 
.134Е+01 
.149Е%01 
m2 STEHT 
65101 
.460ξ101 
559ΕτΤΟ 
ЭССЕСІ 


F-PHI 


GN GN мо нс С) 
Хо Ун ч о 


‚81 


ΧΧΧΧΧ SCATTERED FIELD COMPLEX AMPLITUDES ΧΧΧΧΧ 


INCIDENCE ANGLE = 


PHI= 


DEGREES 


STAGGERED RAD FLDS: 


THETA 


£N TN C4 Са Гу) IN) 1 CO O КО КО Оо 09 — ON 


ә MÀ 3 M M ә 2 2 مر دمو‎ 


мл © Ln O ی‎ O Q O Ln O ہا‎ O Ú, O n O Q 


— — — — — — 
ی‎ J ON ON Un n 
Qn O On O Qn O 


180. 


MAG 


"ESSE SS 
ОЗЕН 
.168E+02 
E 
.105E+02 
τσ ο στο] 
„524Е+01 
.475Е%01 
.515E+01 
„567 Е+01 
ТБ ЛЕНО 
,697E*01 
„451Е+01 
SG TOES 
.596E*01 
.S66E*01 
σετ 
PP SSETOL 
„215 ЕТО! 
(227 ЕЖП 
.275E+*0] 
„314Е+01 
ПЭТ ЕТОЙ 
T EEDI 
.203E+01 
L2 AES 
ΙΙ ΕΠΙ 
215 ЕЕП! 
.295E*01 
„349Е+01 
237 $E E01 
.396Е+01 
„421Е+01 
.462E*01 
οἱ τοι 
ЭЭГЕ ЕО 
. 568E*01 


Pane TA 


DEGREES 


TE INCIDENCE 
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MAG 


foo LeE=06 
819606 
./ 28E=06 
2560E*U6 
.404Е-06 
215E 06 
1497 E=07/ 
ІШ ШЕР ШӘ 
па ЗЕ ЗОО 
„264Е-06 
12662 206 
25E 06 
2158 06 
Pon SESUA 
ШОТЕСІЛІ 
(912Е-07 
Eo SED 
2154Е-06 
ILSE 06 
.814Е-07 
ΕΠΙ 
ШЕТЕ- 07 
ИЗЗЕС- ШІ 
297 2Е-Ш7 
ӘСТЕ ӨЛ 
(8125-07 
.828Е-08 
529Е-07 
.565E-07 
ο Εμ 
.535F-07 
2653E-07 
"TOLE >06 
152E 06 
22016706 
1255 Е-06 
28-06 


F-PHI 


КО С) (0) С) С) (3) С) С) (3 С) (3 () С) (3) С) С) С) С) С) () (3) () (2? (С) (? С) С) С) С) () СС) С» С) <) С) С) С) С) (С) С) С) 


АРРЕНОТХ С 
ENCADESGURCE CODE 


PROGRAM EMCAD 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


RES ROGRAMSCOMPRUTES EM FIELD SCATTERING FROM LAYERED 

Ρο πες Ob SREVOLUTION USING THE COUPLED AZIMUTHAL POTENTIAL 
CCAP) FORMULATION IN CONJUNCTION WITH A VARIATIONAL FINITE- 
ELEMENT TECHNIQUE AND A TRI-REGIONAL UNIMOMENT METHOD. THE 
NUMERICAL SOLUTION IS PERFORMED USING A TWO-SWEEP RICCATI 
TRANSFORM WITH CIRCUMFERENTIAL MARCHING. UNIT ITS USED FOR 
SEQUENTIAL UNFORMATTED STORAGE OF GENERATED RICCATI ARRAYS TO 
DES Ες VERSION OF -CDC-7600 U.C. BERKELEY CODE "EMSCAT" 
WITH LAYERED MEDIA ADDITIONS. 


MODS BY M.A. MORGAN MAR 87-FEB 88 
EM: ΟΟΜΠΜΟ T MAR 62 FEB ὅδ. 


XXXXXXXXXXXXXXXxxxxx EMCAD PARAMETERS XXXXXXXXXXXXXXXXXXXXXXxxx 


xx NA = NO. OF INCIDENT ANGLES, DALPHACJ) (DEGREES) 

Xx DMIN,DMAX - MIN AND MAX MESH DENSITY CELEM/LAMBDA) 

Xx NR1,NR2 ~ RADIAL NODE SPACING BETWEEN RA-R1 & R2-RB 
xx LET SMO OF LAYERED SURFACES, _ NLAYtI IS FREE SPACE 
XX ERCNLAYO,URCNLAY) - COMPLEX ER AND UR OF PENETRABLE LAYERS 
xx MSTART,MSTOP - AZIMUTHAL MODE RANGE ON "M" 

xx MNE UPPER CIMITS ON INIERIOR S EXPBERIORSMODES "NH" 
xx ИЕНЕН = ЧО. OF FIELD PRINTOUT THETA & PHI POINTS 

XX NB,NBMX - NO. OF INTERIOR NODAL POTENTIALS/THETA STEP 
Xx NCMAX ~- NO. OF INTERIOR + EXTERIOR EXPANSION COEF'S 
xx МММАХ - NCMAX + 2XNA (NO. UNIMOMENT BC'S) 

xx МІНС,МІММХ - 2XNA (2 POLARIZATIONHS/IHC AHGLE)D 

xx DOSS O TOR 1I : ZTAATISSHOMOSPRNRSPCHÜUEPCOR CAUCHYI BC'S 

xx IMAS IBIG = NO. OF THETA STEPS IN MARCHING & DATA FILE 


XXXXXXXXXXXXXXXXXXXX ARRAY DIMENSIONS ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


Xx R,RINV - (NBMX,NBMX) RICCATT MATICES 

xx S,SINV - (ΝΒΜΧ,ΝΜΜΑΧ)  S-VECTOR ARRAYS 

xx FMOM,DUM - (NCMAX,HCMAX) COEFFICIENT MOMENT ARRAYS 

Xx ОМОМ, СОЕЕ - (МСМАХ, МТММХ ) DRIVING & COEF ARRAYS (SAME) 
xx ЕТ1 ТО ЕР2 - (НТ,МР, МА) SCOITERED E=FIELDS 

xx RAD = (2,5) RADIICTHETIA ) ΟΡ ΓΑΤΕΣ ΑΤ Т & I+] STEPS 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


CH O RACIERXI BELL,FF,LEVEL 

CHARACTERX8 GRFDAT,PRTDAT 

CHARACTERX12 GRFDTI,GRFDT2,GRFDTS,GRFDT6,PRTDTI 

CHARACTERX64 DUMMY,GRAFLAB,STRDAT,HDR,MATPAR 

ΤΕ А,ВСО,Р,Е1,Н1,Е2,Н2,ЕВ1,НВ1,ЕВ2,НВ2 

CELE. F50; F2ol,F12,F20;F11,F10,F1,F2,ET1, FP1,FT2,FP2 
COMPLEX R(58,58),S(58,80),RINV(58,58),SINV(58,80),ER(6),UR(6) 
COMPLEX FMOM(70,70),COEFC70,10) 

INTEGER BC,ELND, IDUMMY 

REAL RAD(C2,5) 
COMMOH/ONE/X(622,YC62) 
COMMON/TVWO/ELNDC60, 5), NDELC62, 35»), NDCTC 62) 
COMMOH/ZTHREE/EB1(3,802,HB1(C 5, 800, EB2( 5,800, HB2€ 5,80) 
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COMMON/FOURZAC58, 
COMMON? ΝΕ ο «ου 
ο ο ως. 
COMMON/SEVEN/SCAT , 
COMMONH/EIGHT/ZALPHAC5) 
СОИМОНСТЕПИЕЗ ЖЕЛГЕ 
EQUIVALENCE CRINV,SIN 
BELL=CHAR(7) 
FF-:-CHARCIZ) 

SETTING MAXIMUM MATRIX DIMENSIONS 

МСМАХ= 70 

МММАХ= 80 

NBMX=52 

NINMX=10 

. 168159; 

DTR=P17180. 

READING INC. FIELD, MESH, AND SOLUTION PARAMETERS 
WRITECX,1357) 

READCX,101) MATPAR 

ПРЕНнСе,ЕРБРЕМАЦФАЕО 


WRITECX,112) 
READCX,101) STRDAT 
OPENC10,FILE=STRDAT ) 
READC2,101) HDR 
WRITE(X,101) 
WMRITE(%x,101) HDR 
READCZ; δα НОЕ 
WRITECX,101) HDR 
REC x5. 25) 
READCX,101) DUMMY 
ИКТТЕСХ, ОПЛ ЕЕ 
NRITECX,10) 
READCX,101) DUMMY 
ΤΕΕ, ον FF 
WRITECX,12) 

READCX TOTO LEVEL 
HRITECX ЛОТТУ FE 


IF(LEVEL εἰ ΓΟ ΙΟ» 
IFCLEVEL E ο 1 Εμ 
WRITE Cerise 
READC¥,101) PRIDAT 
Е-Е 
WRITECX,13) 
READCX,101) PRTDAT 
ENDIF 
PRTDII = PRTDATZ2Z' QUT: 
OPENCS,FIUESRSETDIIDO 
LFCLEVEL.EQ СОРИ ТЕН 
HRITECX,114) 
READCX,101) GRFDAT 
ELSE 
WRITECX,194) 
READCX,101) GRFDAT 


ENDIF 

ОКЕШТІ - OREDATZAC DPI 
GREDI 2 = ОКОВА ΠΕ; 
GRFDT3 = СКЕВАШИ τει 
GRFDT4 = GRFDAT//'.TEP' 


OPENCG,FILE=GRFDT1) 
OPENC7, FILE-=GREBIZ) 
OPENC8, FILE=GRFDT 3) 
OPERC9 ,FILESORFDIS) 
ΓΓ(ΓΕνΕΙ ΕΠ; ἩΝΡΊΕΙ͂ 
HRITECX,LIS59 
READCX,101) GRAFLAB 
БЕСЕ | 
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22 


53 


44 


WVMRITE(x,15) 
READ(x,101) GRAFLAB 
ENDIF 
WHRITE(4,101) GRAFLAB 
WRITEC7,101) GRAFLAB 
ИКТТЕС8, 101) СКАРГАВ 
πατε θα) GRAFLAB 
ο ο ΕΕ οἳ. ο THEN 
HRITECX,116) 
READCX,X) DMIN 
КТ ТЕСЕН) 
READCX,X) DMAX 
ELSE 
WRITECX,16) 
READCX,X) DMIN 
ЧКІТЕСЖ, 17) 
κΕΑΡ(Χ,Χ) DMAX 


ENDIF 

wets, LOL) FF 
WRITEC3,100) 

WRITECS,X) 

ο τες 11) 
READ(2,102) ΝΙ ΑΥ 

DO 22 L=1,NLAY 
READ(2,103) ER(L),UR(L) 
WRITEC3,107) L,ERCL),URCL) 
CONTINUE 

БЕЗШОЛП» 150) IBIG 
READC10,103) RMIN,RMAX 


ШЕСОШЕЕС-ЕЦШС ТО"). ТНЕМ 
WRITECX,118) 
READCX,X) МА 

ELSE 
HRITE(x,18) 
READCX,X) NA 

ENDIF 

ΝΙΝΟΞΟΧΝΑ 

SAMX=0.0 

WRITEC3,109) 

DO 33 I=1,NA 

HRINREGX S23) [ 

READCX,X) DALPHACID 

ИКТТЕС 5,110) ГТ, ОА! РНА(Т) 

ALPHACI)=DTRXDALPHACT) 

SA=SINCALPHACT)) 

IFCSA.GE.SAMX) SAMX=SA 

CONTINUE 

ο ΕΕ ΕΕ. το. ο) THEN 

HRISPECX, 1195 
READCX,X) NT 

ELSE 
WRITECX,20) 
READCX,X) NT 

ENDIF 

ЕЕ ВЕЧЕ ЕО 21} THEN 
WRITECX,120) 
READ(CX,X) NP 

ELSE 
NRITECXxX, 21) 
READCX,X) NP 

ЕШТЕ 

DO 44 Ізі,МР 

WRITECX,24) I 

READCX,X) DPHICID 

PHICIO-DTRXDPHICI) 

MSTART=0 
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MSTOP=INT(RMAXXSAMX)+1 

HRITEC(X, 121) MSTOP 

READ(x,x) MSTOP 

IF(SAMX.LT Ι ετσι) MSTARTSI 

SETTING RADIAL STEPS BETWEEN RA ТО ЕТ ДБО ЕН 

МЕІзі 

HR271 

DEFINING OUTSIDE ο ατα WITH SLIGHT LOSS 

URCNLAY+])=Cl. 03-1 ВЕ со 

ERGHLAY +1 )=C1 ОТЕТ m 6) 

ESTIMATING NUMBER OF CORE AND EXTERIOR MODES 

NI=INTCRMINXCABSCCSQRT(2.0XER(1l1)XUR(1))))+1 

N2=INT(RMAXxCABS(CSQRT(2.0OxER(NLAY)%XUR(NLAY))))+1 

WRITECX,122) NI 

READ( X, x) N1 

WRITEC A, tess NZ 

READCX,X) N2 

WRITECX,X) BELL 

PAUSE 'Check that Nl*N2 is .LE. 35 .... Otherwise Abort' 

IFCCHL*N22.6T.35) STOP 

WRITECSII0S5) MSRARISMSTEOBSNTANZ 

WRITEC*,129) 

M1=MSTART+1 

M2=MSTOP+1 

ММ-М5ТОР-М5ТАЕТ-1 

IFCMM.LE.0) MM=1 

DO 77 MD=M1,M2 

OPENC1L,FILE='RS. DAT’, STATUS="UNKNOWN", SECESS="SEQUEI TAs 
FORM="UNFORMATTED! ) 

DD=DMAX-DMIN 

M=MD-1 

WRITEC3,101) FE 

МЕІТЕ(С5,204) М 

WRITECX,X) BELL 

WRITECX,1252) M 

BC=0 

IF (MEEQ 17 ες-] 

MF =0 

IF(M.EQ.0) ΜΕΞΙ 

COMPUTING MODAL EXPANSION TRUNCATION LIMITS 

НМ1<Н1 + 1-М-МЕ 

ІЕСММІ.ІТ.1) ММ1=1 

NM2=N2+1-M-MF 

IFCNM2.LT.1) NM271 

HCOEF=2%x(NM1+NM2 ) 

HMODES=NCOEF+NINC 

COMPUTING MESH PARAMETERS AND RADCI) FOR GIVEN MODE = M 

DENS=DMAX-DD¥(M+MFE-1.0)/MM 

WRITEC*,126) 


CALL MESH (CDENS,ER,UR,NR1,NR2,IMX,NDS,LMX,NB, DR, DT, 
] 


RA,RI,R2,RB,RMIN,RMAX,NLAY) 
CHECKING MAXIMUM MESH PARAMETERS 
IE CNS LE “NEMA De COS o> 
WRITEC3,104) 
STOP 
CONTINUE 
LOADING AND MARCHING ((R)) AND ((S)) 
МКТ ТЕС КОТО ДВЕ 
ИКТТЕС 5,201) 
WRITECX,127) ІМХ 
DO 66 I-1,IMX 
ІК-ІМХ-І 
WRITECX,128) I,IR 
ΗΕ ΤΕΣ, 


CALL LODER (I,M,BC,ER,UR,NMODES,NM1,NM2,NINC,NB,IMX,IBIG, 
1 


NDS,RA,RB,R1,R2,DR,DT,NR1,NR2,RAD,NLAY) 


100 


E 


66 


77 


10 


12 


IFCI- GT-11) GO TO 57 

ШЕТТЕСЪ З00ОН И ККАШСТ БО ЕТ, ЧЕАТ) 

CONTINUE 

τες S00) PX CRADC2 SD) EZISRUAY) 

HEITECX,150) 

CELL MARCH CIT IMX, BC МВ, ЧВМХ, НМОВЕЗ, НИМАХ, Е, 5, КТНУ) 

CONTINUE 

ERISSECX, 151) 

CALL SWEEP CM, IMX,BC,NML,NM2,NINC, NMODES, NMMAX,NCOEF, NCMAX, 

1 NRI,NR2,NB,HBMX,RI,R2,DT,ER,UR,R, S, SINV, FMOM, COEF, NINMX) 

WRITECX,I32) 

CALL MODECO CNMI,NM2,NINC,NCOEF,NCMAX,FMOM, COEF, M, NINMX) 

WRITEC(X,133) 

CALL RADFLD CM,NMI,NM2,NT,NP,NA,NCMAX,COEF,NINMX) 

WRITECX,134) 

CALL CRSSEC (М,ММ1,ММ2,МА,МСМАХ,СОЕҒ, МІММХ) 

REWIND 10 

READC 10,138) IBIG, IDUMMY 

READC 10,103) RMIN, RMAX 

CONTINUE 

ον ο. o TATUS: DELETE") 

CLOSE(Z2) 

CBOSECIO) 

BHRITECX,135) 

CALL DATOUT CMSTART,MSTOP,;, NT, NP, NA) 

WRITECX,X) BELL 

МЕТТЕСх, 156) 

ҒОРМАТ(// ///////////////7Х,!ххжххЖ MWHELCOME TO EMCAD Xxxxxx',77/7X, 
l1'THIS PROGRAM COMPUTES EM FIELD SCATTERING FROM PENETRABLE',77X, 
2 'BODIES OF REVOLUTION USING THE COUPLED AZIMUTHAL POTENTIAL',77xX, 
3SUCCAPO FORMULATION IN CONJUNCTION WITH A VARIATIONAL FINITE-',77X, 
G'ELEMENT TECHNIQUE AND A TRI-REGIONAL UNIMOMENT METHOD.  THE',77X, 
SD NUMERICAL SOLUTION 15 FERFORMED USING A TWO-SHEEP RICCATI',77xX, 
6'* TRANSFORM WITH CIRCUMFERENTIAL MARCHING.',77X, 

POUR DATED VERSION OF ORIGINAL U.C. BERKELEY CODE.',77X, 
8'MODIFICATIONS BY M.A. MORGAN MAR 1987 - ҒЕВ 1988!,/7Х, 

9! ΕΠ ΟΠ C JULDR1987 APR 1988 .'!,ZZ/ZZ7ZZZZ7; 
TTELEASE PRESS ANY KEY TO CONTINUE." ,7) 

ҒОРВМАТС// / / ////////////7Х, ΧΧΧΧΧΧΧ EMCAD INPUT XXXXXXX',//7X, 
ШРЕПСАр АГГОЫ5 А USER TO IHPUT THE NECESSARY INFORMATION',/77X, 

2 ACCORDING TO THE LEVEL OF HIS EXPERTISE.",777x, 

S UD UICE LEVEL — (GIVES@MORLEF EXFCANATIONS/DESCRIPTIONS!',;/7X; 
4! OF INPUT VALUES. DESCRIBES THE REQUIRED',/7X, 
5 FORMAT FOR THE INPUT VALUE. GIVES ",/7Х, 

6! TYPICAL UELUES MNERE APPLICABLE 777/7Х, 

VTE EXEERT LEVEL = ASSUMES THE USER IS FAMILIAR WITH THE',77X, 
o ΠΤ PARAMETERS AND FORMATS. SIMPEY',77X, 

9! PROMPTS FOR REQUIRED INPUTIS . !,ZZ7X; 

MS EXIT EMCAD!,; 777A, 

eP EA S ERS ELECT THE LEVEL OF EXPERTISE BY ENTERING 1, 2, ОК 5.! 
3⁄⁄7X) 


TS РОРМАТСУ/////////////////////7Х, 


ΠΕ OUTTUT DATA FILE IS THE FICE CONTAINING THE OUIPUT! ,Z7X, 

e RESULTS OF ALL NUMERICAL CALCULATIONS CONDUCTEP BY -EMCAD.',77X, 
SINE ΟΡΙΣΕ FOR THIS INPUT IS FILENAME ONLY. NO EXTENSION IS',/7X, 
G'REQUIRED OR DESIRED. EMCAD AUTOMATICALLY APPENDS AN ',77X, 
ΤΕΙ ΙΟ OF ου IO YOUR FILENMAME. T Xs 

ΠΕ ΕΝΙΕΙ ΤΗΕ FILENAME GF-THE- OUTPUÜT-DATA FILE.',777X) 


16% ЕО®ЁМАТС////////////////////////7Х, 


II TTHEFOUTPUT GRAPHICS DATA FILE IS THE FILE CONTAINING THE ',77X, 
e'QUTPUT DATA FROM EMCAD TO BE USED AS INPUT TO GRAPHING ',77X, 
SeSCULINES. THE FORMAT FOR THIS INPUT IS FILENAME ONLY.  WNO',77X, 
G'EXTENSION IS REQUIRED OR DESIRED. EMCAD AUTOMATICALLY ',/7/7X, 
S OP E IDSPAH E EXIENSIQH DE ITMI ОТЫРУ bet, AND „ПЕР TO YOUR',/77X, 
ο ΕΕ ΜΕ ΛΣ ТТ PRODUCES FOURTOUTPUT FILES FOR GRAPHICS.',777X, 
"EMI = ο) РЕНЕ ПРЕНЕСЕ ШЕЕ ТНЕТА АРХ, 
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8'.TMP ----- >TM INCIDENCE, F=PHI' 7 
ο: TET ----- > TE INCIDENCE, F-THETA TE s 
TE tie > ТЕ INCIDENCE, F=PHI' A7 7:5 
2'PLEASE ENTER THE FILE HAME OF THE OUTPUT ПАЛЛАС 

15 ΕΟΡΜΑΤέ ΑΟ ο τα... 
1' THE GRAPHICS CAPTION IS A PERSOHALIZED CAPTIQN ο ο σι D p 
2'USER TO IDENTIFY TAIS SET OF GRAPHS FROM ALDO IHER SES DD P 
S'THE MAXIMUM LENGTH OF THIS CAPTION IS 664 CHARACTERS. ',77X, 
+" НОТЕ: WHEN USED HITH THE GRAPHICS PACKAGE, THE PROGRAM IS',77X, 
5'ABLE TO DIFFERENTIATE BEIWHEEN UPPER CASE AND DOWNER CASE TETEE 
6 'CHARACTERS ID O A 
7 'PLEASE ENTER ANY GRAPHICS CAPTION YOUTDESIRE T777) 

16 FORMATCZZZA/7Z ТУ ТСО 
l'DMIN AND DMAX ARE PARAMETERS OF MESH DENSITY IN TERMS 0Е!,/7Х, 
A'ELEMENTS/INTERIOR LAMBDA. INPUT VAUUES ARE EXPECTED TO BETTS 
S'REAL, І.Е. ІНЕ DECIMAL POINT MUST BE ΙΝΕ ΓΕ)... 
4'TYPICAL VALUES ARE DMIN = 10. AND ОМАХ = 15.',//7Х, 
4'PLEASE INPUT DMIN',//7X) 

17 FOPMAT(ZZ7X,; 'PLEASE INPUT DMAX5 5 στ. 

18 ΕΟΡΒΜΑΤί,///  7Χ, 
1]1"'THE NUMBER OF INCIDENT FIELD ANGLES IS ΤΗΕ το μον... 
2' INCIDENT FIELDS THAT IMPINGE ON THE OBJECT OF INTEREST И 
S' THIS PROGRAM ALLOWS A MAXIMUM OF FIVE INCIDENT FIELD ',772X5 
G'AHGLES. WHEN ENTERING YOUR ANSWER PLEASE DO NOT INCLUDE A'*,/7xX, 
S'DECIMAL BECAUSE THE INPUT MUST BE IN INTEGER FORMAT, I.E 7 
S МА 53: , GNA = | 7XS 
6'PLEASE INPUT THE NUMBER OF INCIDENT FIELD ANGLES ee тх 

19 ΕΟΡΜΑΤί(, .,”χ, 
1'THE INCIDENT ANGLES DEFINE THE DIRECTION FROM HHICH THE !,Z7X, 
Q'INCIDENT FIELDSOORIGINATE 2X 


31 DEGREE  ===== > + 4 DIREQGIION', 77X; 
ΟΠ DEGREE =s > + Х DIRECTION АТ” 
S ISU DECREE Sas > = Z DIRECTION АСА 


20 ЕОЁМАТС///У/У/У/ УУГ /////////////7Х, 
1' THE NUMBER OF SCATTERING FIELD THETA POINTS DETERMINES THE ae 
2'SPACING BETWEEN THETA POINTS DURING EMCAD ITERATIONS AND',/7X, 
5S * CALCULATIONS: 27/70 
G'DELTA THETA = 180 7 CNUMBER THETA POINTS = 1) 3502.7 
5' NUMBER THETA POINTS = 37 ----- > DELTA THETA = 5 DEGREES 2400 
6 ' NUMBER THETA POINTS = 19 ----- > DELTA THETA = 10 DEGREES 7 
7'WHEN ENTERING YOUR ANSWER, PLEASE DO NOT INCLUDE A DECIMAL ',77X, 
8'BECAUSE THE INPUT MUST BE IN INTEGER ЕОЕМАТ. Г.Е. ИГІ: 19777770 
9'PLEASE INPUT THE NUMBER OF SCATTERING FIELD THETA POINTS" У 

21 ΕΟΠΜΑΤέ ΙΛ οτι. 
1 "ТНЕ NUMBER OF PHI ANGLES IS THE TOTAL NUMBER OF PHI ANGLES.: 177X; 
2' THIS PROGRAM ALLOWS A MAXIMUM OF EIGHT PHI ANGLES. WHEN ',77X, 
S'ENTERING YOUR ANSWER PLEASE DO NOT INCLUDE A DECIMAL ',77X, 
4" BECAUSE THE INPUT MUS) Be Tie tNnecer ГОЕМИИРТТТЕС Е 
Б! NP =< 3 5 NE - 7 x; 
6'*PLEASE INPUT THE NUMBER OF PHI ANGLES TAAT) 

à3 FÜRMATC/Z7X,'ENTER INC FLD ANGLE {ΠΕΡ БОК ТГ 

20 FORMATC/7X,*ENTER PHI FLD ANGLE CDEG) FOR 8", 13,.777%) 

25 FORMATC//7/777X, *PLEASE VE RES SoeANY КЕТ πα ΕΟΟ ΝΕ АА Т 

100 FORMAT CZ, ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ EMCAD OUTPUT DATA XXXXXXXXXXXXXXxx 
] XXXXXXXXXXX') 

101 FORMATCA) 

102 ҒОРМАТ(І5) 

102 ҒОКМАТ(48(Е14.6)) 

104 FORMAT(' xxx PROGRAM ABORTED BECAUSE NB > NBMX Xxx!) 

105 FORMATC/7X, 'SERIES SOLUTION PARAMETERS τε Ιστ ОБУ 
1'М5ТОР= ',15,/77X,'Nl= 1 Ακ ο... 

107 FORMAT(ZZ7X,15,2(2X,'(C tU ει ο ο ΤΕΙ ΙΙ 

109 FORMAT(ZZ7X, 'INCIDEIII FIELD 'ANGLES' C C O Oo Ox Veena Noe 

110 FORMAT(ZZ7X,13  7X7F5% 1 ΠΕΙ. 

111 FORMAT(/7X, COMPLEX Er) AND Un m 

112 FORMAT(//7X, "ENTER STRUCTURE DATAT FILES D FILENAME EXTENSTONDE и) 


bro FORMANGZ47x%, "ENTER OUTPUT DATA FILENAME HITHODESESTVENSION',Z7ZT7TX, 
TOC TRIMER ΕΙΤΕ ENTER ο στι ") 

114 FORMATC/4/7X, "ENTER GRAPHICS DATA FILENAME WITHOUT EXTENSION 2 

O FORHATC(ZZ7X, 'EHIER САРТТОП ОК LABEL E) 


EIOSEURMATCZZZXU'EHMTIER BHIM OCELEMEMTSZIHTERIOR LAMBDA): 52 

ШІП ЕОКИАТС//7Х) ТЕНТЕК ΜΑΝ CELEMENISZINTERIOR LAMBDA): n) 

ο ου οκ СТА ТЕППЕ DUNDEE DT ΡΙΕΙΡΕΝΤ FIELD ANGLES C.LE. 5) us) 

ΓΗ, ΕΤΕΙ NUMBER UF THETA POINTS FOR SCAT FIELD (.LE 37.) 
1',4X) 

meu FORMAIC/77X%, "ENTER NUMBER OF PHI ANGLES C.LE. 8): E 


ПВИЕОКМАТСИ7Х, "ENTER MSTOP (.LE. 13)',777X;, 
ESTIMATED MLN MUM VALUE І5:",6Х,14,//7Х) 

22 5 ΝΟ. ΟΕ INTERNAL AND EXTERNAL EXPANSION ΜΟΡΕΣ 
Πχ, 
САК ЛК АГЕШ κ "MINIMUM INTERNAL NI *',2X,164,777X; 
ШЕПТЕРІН ο. М5ТОР);: κ Z УА 

MS OURMATCZZ7X7UESTIMATED (KOXRMIN) U UMINIMUM")ƏW EXTERNAL N2 =!,2X,14, 
Ш 7, ТЕНТЕК М2 δε MSTOP): ЛУ А) 

124 FORMATC//7X, * INDEXING PROGRAM THROUGH VALUES OF M*) 

Meow rFORMAIC7A,°M-LOOP ....™M = ?,15) 

LDSOLCEORMATC/X, 'CALL MESH') 

I 7 FORMAT(?X, "ENTER I-LOOP; NO STE ES TO COMPLETE: I6) 

μα ο TIST 2х, 15,7, "МО 5ТЕР5 Т0 б0:',2Х,‚15) 

IO EP PORMAT(7X, 'CALL LODER!) 

ο ΟκΙα ο κ οι MARCH?) 

more ORMATC7X, "EX I-LOOP, SNEEP') 

E -NPORMATCZX, CALL MODECO') 

BS FORMATCZX,'CALL RADFLD') 

ШИ СВКМАТСТХ, "СА(ЦІГ СК55ЕС"”) 

ior ORMAIC7X,'EX M=-LOOP, ВАТ09Т”) 

156 FORMAT(//7X, 'XXXXXXXXXXXXXXX EMCAD COMPLETED XXXXxxxxxxxxxx"'!) 

137 FORMATC//7X, "ENTER MATERIAL PARAMETER FILE CD:FILENAME.EXTENSION) 

° 3 


: ) 
138 FORMAT(C2CIS5Ó) 
BASE UKRMAT(ZZ7X,' RICCATI MARCHING STEPS"') 
208 FORMAT(//7X,! XXXXXXXXXXXXxx M^ ',I5,!' Xxxxxxxxxxxxxx'//) 
ЕОМРГОКМАТСТ? ХА, І5,5Е12.5) 
99 5ТОР 
END 


The above constitutes only the main Program of EMCAD. The 
remaining subroutines would take approximately 26 more pages to 
print. The entire program code can be obtained by sending an 
IBM PC compatible formatted 5 144 " floppy disk То 


Prof. M. A. Morgan 

Code 62 MM 

Haval Postgraduate School 
Monterey, Ca. 93943 


APPENDIX D 


EMCADIN USER'S MANUAL 


Welcome to EMCADIN. This program is designed to accept 
properly formatted input data from the West Coast 
Consultants commercial shareware product CURVE DIGITIZER 
(Appendix M) and convert it to a form which can be used by 
the scattering calculation program EMCAD. CURVE DIGITIZER 
is a computer aided design (CAD) package which will allow 


user manipulation of the data. 


A. BACKGROUND 

EMCADIN was written in the programming language of Ryan- 
McFarland FORTRAN (Appendix M). Any questions concerning 
formats, changes to the program, etc... should be referred 
to the respective user's manuals, Dr. M. A. Morgan or LT. E. 


M. Connolly (Appendix M). 


B. COMPUTER SETUP 
EMCADIN was not written for any specific computer 
configuration, and as such, should run on most ІВМ 


compatible machines. 


C. PROGRAM INPUT AND OUTPUT 

EMCADIN obtains its input from a previously generated 
data file. The input data file for EMCADIN must be located 
in the same directory/subdirectory as EMCADIN. This data 
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file is produced through the utilization of the computer 
aided design package, CURVE DIGITIZER. As shown in Figure 
D-1, the output of CURVE DIGITIZER is the input to EMCADIN. 
For EMCADIN and CURVE DIGITIZER to be compatible, it was 
necessary to establish and adhere to the predefined data 
format of CURVE DIGITIZER. In order for EMCADIN to be 
utilized successfully, it is imperative that the user 
understand the need for the particular format being used, as 
well as the format itself. The transition from CURVE 
DIGITIZER to EMCADIN provides more of an opportunity for 
catastrophic error than at any other point in this entire 
package of interactive computer programs. 

Understanding the format of the input of the data for 
EMCADIN necessitates a thorough understanding of some of 
the idiosyncracies of CURVE DIGITIZER, referred to 
henceforth for brevity as CD. Computer aided graphic design 
deals in the basic business of handling points. Any graphic 
figure can be represented as some combination of points. 
Lines can be formed by connecting two points. Arcs can be 
represented by placing more points closer together and by 
approximating the curvature with sufficiently small line 
segments between points. The required input-output format 
for CD is utilized to insure that adequate information is 
present to direct CD how to connect the points. A similar 
process may be that of the simple children's game connect- 


the-dots. A series of dots exists which define a picture, 
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ШЕ connected properly. As more and more dots are used, 
more complex pictures with better resolution can Бе 
constructed. The fundamental assumption employed by CD is 
that ТИЯ should be connected to the next point listed 
in the data file unless some code is present to indicate 
otherwise. For simplicity and unrestricted accessibility, 
CD input and output data are represented in terms of XY 
rectangular coordinates. CD flags the break between two 
data points by what shall be referred to as a separation 
code. Since two coordinates, X and Y, are required to 
establish a data point, two separation codes, one for the X 
data and one for the Y data are required to separate two 
distinct data points. The separation code employed is of 
the form 999990. The last two digits, 90, may change to 
reflect scaling, color, or any of a number of other things. 
It was not possible to determine the specifics from the user 
manual. In the specific case of EMCADIN input data, 999990 
999990 codes are required to distinguish the structural 
data describing one layer from that describing another 
layer. f 

CD stores the data in the exact order in which the data 
points were drawn. There are other methods of generating 
the simple geometric shapes demonstrated in the next two 
examples, however, these examples were specifically chosen 
to illustrate potential pitfalls that could occur if extreme 


care is not taken. For the first example, the simple square 
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of Figure D-2 was constructed utilizing CD. The base of the 
square was drawn by defining point number 1, and then by 
drawing a line to point number 2. The top was likewise 
SCARE E by defining point number 4, and then by drawing 
a line to point number 3. The two sides were then formed by 
connecting point number 4 to point number 1 with a line, and 
by connecting point number 2 to point number 3 with a line. 
This is one of many methods of connecting the four data 
points that define this square. Figure D-3 lists the output 
data file which contains the CD defining data for this 
square. Notice there are significantly more than four data 
points describing this object which requires only four data 
points. Recall from earlier discussion that  EMCADIN 
utilizes the 999990 codes to distinguish structural data 
describing one layer from that describing another layer. 
Figure D-2 depicts only one layer, yet EMCADIN using the 
data file of Figure D-3 would interpret the data as four 
distinct layers. The 999990 codes of Figure D-3 are present 
due to the order in which the data points were entered. Had 
the square been generated by defining point number 3, 
connecting it to point number 2, then to 1, then to 4, 
there would have been no 999990 codes separating the data, 
and the data file of Figure D-3 would consist of only four 
pairs (XY coordinates) of data. 

The second and third major errors that might occur can 


also be demonstrated using Figure D-2. The sole 
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justification for utilizing CD was its ability to generate 
data to be used ultimately by the scattering program EMCAD. 
The main body of the thesis discusses the requirement of 
EMCAD τος axisymmetric bodies of revolution. Any three 
dimensional axisymmetric body of revolution can be defined 
by knowing the axis of revolution, and then by defining the 
body in planar terms. As an illustration, a given sphere is 
a three dimensional body of revolution. The sphere has 
axial symmetry and can be described by a semi-circle in any 
meridian plane intersecting the axis. Rotating this planar 
semicircle about the axis would produce the original 
sphere. EMCAD requires this meridian representation of the 
three dimensional body as its input. For simplicity's sake, 
the CD representation of the meridian plane will be referred 
to as the XY plane, with X being the horizontal axis and Y 
being the vertical axis, as well as the axis of revolution. 
This coordinate system is for discussion purposes only, and 
is not to be confused with the three dimensional coordinate 
system of EMCAD theory. Viewing the body of Figure D-2 as a 
two dimensional body, the Y axis is the required axis of 
revolution. Therefore all CD planar representations for 
EMCAD should be constructed with values on the X axis 
restricted to positive values. 

The third problem occurs when removal of data points is 
attempted in order to obtain a proper planar representation. 


EMCADIN performs a valuable service to the user of improving 


1: 


the resolution of input data. It suffices to say that due 
to the inherent mathematics involved, EMCAD requires that 
two data points be defined, both on the Y axis. The first 
data ου defines the beginning of the object, at X = O, 
and the last data point defines the end of the object, also 
at X = 0. In the example of the square, the beginning data 
point would be (0,2), and the ending data point would be 
(0,72). Without these beginning and ending data points for 
each layer of data, the EMCADIN interpolation will "bomb", 
and will produce erroneous results, if any results at all. 
For the second example, the simple circle of Figure D-4 
was chosen. CD allows the user to construct a circle by 
simply defining the coordinates of the center of the circle, 
and by specifying the desired radius. In every instance, CD 
generates the data points by starting at point 3 on Figure 
D-4, progressing to point 2, to point 5, to point 4, and 
back to point 3 in a counter clockwise pattern. A sample of 
the output data file for this circle is shown in Figure D-5. 
This data, for a chosen radius of 2, reflects the same 
counter clockwise pattern, which is directly opposite to the 
pattern required by EMCADIN input data. In EMCADIN, theta 
is defined in a clockwise pattern, starting on the positive 
y axis and progressing by equal step sizes clockwise to the 
negative y axis. By using this method of circle 
construction, the user must find someway to edit and 


rearrange the data before it can be used by CD. The other 
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problem that arises, particularly with circles, is that as 
the radius becomes smaller and smaller, fewer and fewer 
points are used to construct the circle, causing the data 
points εσας less and less accurate. 

For the case of multiple layered bodies,  EMCADIN 
expects that the data point values for each layer have been 
entered in the input data file starting with the layer 
closest to the center of the object, labeled as Layer 1 in 
Figure D-6, then advancing outward from the center until all 
layers have been accounted for. It has already been stated 
that data points are listed in the order in which they are 
entered, which means that not only might data points in a 
given layer be out of proper order, but data points in 
differing layers, or even the layers themselves might be 
interspersed. At the present time, there appears to be only 
two methods to preclude the abundance of possible errors. 


The first of these dictates the use of extreme care when 


entering the data points. The values must be entered in 
order, starting at Ө - 0 degrees and progressing clockwise 
to @ = 180 degrees for the layer closest to the center 


point, and then entered for the next closest layer, and then 
the next closest layer ... until all layers have been 
described. Even a small band defined by a different 
material parameter must be defined for all 6e from O degrees 
to 180 degrees. If this process is followed exactly, there 


is no need for the second method. The second method 
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involves the use of some text editor to rearrange the data 
as necessary, remove unnecessary 999990 codes, insure data 
for layers as well as layers are in the proper order, etc... 
The ident method would include extreme caution coupled 
with a data review utilizing a text editor. 

The EMCADIN program output data consists of the number 
of layers of the body being considered, the minimum and 
maximum values of radius, and the listing of radius values 
in a clockwise manner for 9 = 0 to 6 = 180 degrees for each 


layer. 


D. PROGRAM EXECUTION 
Once the correct disk or directory has been selected, 
EMCADIN can be invoked by typing 
EMCADIN 
at the DOS prompt. The program will load itself and the 


user will see the following screens: 
* Welcome to EMCADIN (Figure D-7) 
This screen presents some brief background information. 
* Input Filename (Figure D-8) 
This screen prompts the user for the name of the file 
containing the input data for EMCADIN. If the file is 


obtained from CD, the file must have an extension of „dat, 


117 


WELCOME TO EMCADIN 


THIS PROGRAM ACCEPTS PROPERLY FORMATTED INPUT DATA 
FROM CURVE DIGITIZER AND CONVERTS IT TO A FORM 
WHICH CAN BE USED BY EMCAD. 


PLEASE PRESS ENTER TO CONTINUE. 


FIG. D 2227015 17 
EMCADIN MENU 
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INPUT FILENAME 


INPUTFN IS THE INPUT FILENAME OF THE FILE CONTAINING 
THE DATA OF INTEREST WHICH WILL BE INTERPOLATED ON AND 
CONVERTED TO A FORM WHICH CAN SUBSEQUENTLY BE USED bY 
САЙ. 


PLEASE ENTER THE NAME OF THE INPUT FILE. THE EXTENSION 
OF THE FILENAME MUST BE .DAT, I.E. FILENAME. DAT 


ο ος. ο ο. ΕΕ ΕΝ Επ T 
ЕБТ n NA 
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and the extension must be entered by the user as part of the 
filename. The input data file must Бе іп the same 


directory/subdirectory as EMCADIN. 
* Output Filename (Figure D-9) 


This screen prompts the user for the name of the file that 
will hold the output data from EMCADIN. The extension of 
this file is left to user discretion, but must be included 


as part of the output filename. 

κ Resolution (Figure D-10) 
EMCADIN uses linear interpolation to provide the capability 
to generate output data with better resolution than the 


input data. This screen prompts the user for the desired 


output data resolution, in the units of degrees. 
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OUTPUT FILENAME 


OUTPUTFN IS THE OUTPUT FILENAME OF THE FILE CONTAINING 
THE DATA WHICH HAS BEEN INTERPOLATED AND CONVERTED TO À 
FORM WHICH WILL BE USED BY EMCAD. 


PLEASE ENTER THE NAME OF THE OUTPUT FILL, THE EXTENSION 
OF THE FILENAME MUST BE INCLUDED I.E. FILENAME DAT 


τ ee ene Ne ROME f 
EE ES TESTE NAME: 
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RESOLUTION 


DELTHE [S THE USER INPUT VALUE OF THE DESIRED THETA 
RESOLUTION IN DEGREES, 


PLEASE ENTER THE DESIRED DELTA THETA VALUE IN DEGREES. 


rib. Ш-ІП опери μα 
РОК OUT ου. s m 
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APPENDIX Е 
EMCADIN SOURCE CODE 


PROGRAM EMCADIN 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WELCOME TO EMCADIN WRITTEN BY LT E. M. CONNOLLY. 

EMCADIN IS A PROGRAM DESIGNED TO ACCEPT PROPERLY FORMATTED 

INPUT DATA FROM CURVE DIGITIZER AND TO CONVERT IT TO A 

FORM WHICH CAN BE USED BY THE PROGRAM EMCAD, WRITTEN BY 

DR. M. A. MORGAN. 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
VARIABLE DEFINITIONS 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


INPUTFN THE INPUT FILENAME OF THE FILE CONTAINING 
THE DATA OF INTEREST WHICH WILL BE INTERPOLATED 
ON AND CONVERTED TO A FORM WHICH CAN 
SUBSEQUENTLY BE USED BY EMCAD 

OUTPUTFN Paes OUTPUT FICENAME OF THE FILE CONTAINING 
THE DATA WHICH HAS BEEN INTERPOLATED AND 
CONVERTED TO A FORM WHICH WILL BE USED BY 
EMCAD 

I, J, К, INTEGER COUNTERS 

L, N 

RADMIN THE MINIMUM RADIUS VALUE 

RADMAX THE MAXIMUM RADIUS VALUE 

kI A CONSTANT 

DELTHE THE USER ΠΕ ο VALUE OF THE DESIRED THETA 
RESOLUTION IH ΡΕΟΡΕΕΣ 

DELTHETA THE USER ΙΡ VALUE OF THE DESIRED THETA 
RESOLUTION IN RADIANS 

THETACN) AN ARRAY HOLDING THE DESIRED VALUES OF THETA 
AS DETERMINED BY DELTHE 

DEGTORAD A CONVERSION FACTOR FROM DEGREES TO RADIANS 

NUMPTSIN IRE NUMBER OF FOINTS HELD IN A GIVEN INPUT 
DATA SET 

HUMSETS THE HUMBER OF INPUT DATA SETS 

DUMMY1 AN IHPUT VARIABLE ALLOHING SCREEN ADVANCING 
DURING MEHU SELECTION 

DUMMY 2 A CHARACTER STRING USED TO LOCATE THE END OF 
EACH INPUT DATA SET 

XARRAY (I) AN ARRAY CONTAINING THE VALUES OF X FOR ALL 
THE INPUT DATA SETS 

YARRAYCT) AN ARRAY CONTAINING THE VALUES OF Y FOR ALL 
THE INPUT DATA SETS 

ACI,J) AN ARRAY CONTAINING THE VALUES OF X FOR ALL 
THE J INPUT DATA SETS 

BCI,J) AN ARRAY CONTAINING THE VALUES OF Y FOR ALL 
THE J INPUT DATA SETS 

NUMPTOUT ΕΞ ΡΕ ος οὐ OUT DETERMINED BY THE 
DESIRED RESOLUTION 

ПНЕТАТ LEE DUE Or THEIR ΡΕΠΕΚΗΙΝΜΕΡ ВТ ТНЕ І-ТН ІНРУТ 
X AND Y VALUES ` 

THETAK THE VACUESO M HETAN DETERMINED BY THE K-TH INPUT 
X AND Y VALUES 

XFOUND THE DESIRED VALUE OF X CORRESPONDING TO A GIVEN 
VALUE OF THETACH) 

YFOUND THE DESIRED VALUE OF tT CORRESPONDING TO A GIVEN 


VALUE OF THETACN) 


E5 


OO 000000000 


ОООО ОООО 


0000000 
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RADIUS(N, J) AN ARRAY CONTAINING THE OUTPUT CALCULATED 
VALUES OF RADIUS FOR ALL THE INPUT DATA SETS 


M SLOPE OF THE LINE CONTAINING THE ITH AD 
DATA POINTS 
BB Y-INTERCEPT OF THE LINE CONTAINING THE ТЕЕ 


K-TH DATA POINTS 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧχχχχχχχχκχ 
VARIABLE DECLARATIONS 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
DOUBLE PRECISION XINPUT(2000), YINPUT(2000), RADIUS(2000,5) 
DOUBLE PRECISION DELTHE, DELTHETA, THETAC2000), THETAI, THETAK 
DOUBLE PRECISION A(C2000,5), BC2000,5) 

DOUBLE PRECISION DEGTORAD, INTERI, INTER2, RADMIN, RADMAX 

DOUBLE PRECISION PI, DUMMY, XFOUND, YFOUND 

DOUBLE PRECISION BB, M 

INTEGER І, Ј, К, L, N, NUMPTSIN, NUMPTOUT, NUMSETS 

CHARACTERX3 DUMMY1, DUMMY? 

CHARACTERX12 INPUTFN, OUTPUTFN 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INDTIAG VALUES 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
РІ - 35.16155 

J = 0 

Кезо 

L = 9 

NUMPTSIN = 0 

NUMSETS =NI 

DEGTORAD = РІ / 180. 

RADMIN = 10000 

RADMAX = -10000 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INPUT DATA OF INTEREST 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITECX,110) 

READ(CX,100) DUMMYI 

WRITE(X,120) 

READCX,100) INPUTFN 

WMRITE(xXx,130) 

READ(x,100) OUTPUTFN 

WRITECX,140) 

READCX,X) DELTHE 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxx 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
DETERMINE LOCATION OF END OF DATA 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
OPENC7, FILE * INPUTFND 

ОРЕНСО,ҒІГЕТ- ΠΤΙ 

READ (7,100) DUMMY2 

MNUMPTSIN = ΠΠ εν 

IFCDUMMYZ .NE. ΕΠΟ, ΠΗ 

CLO5FEF(7) 

OPEN(7, FILE = INPUTFN) 

М< | 

J = 1 

DO 20 I= 1,NUMPTSIN =- 1 

READ (7,ж) ХІМРОТСІ), ҮІМРОТСІ) 

IFCABSCXINPUTCIO — 999990) CI ον ПЫ! 

A(N,J) = ХІМРОТСІ) 

= FINPUTCI) 


BCN, J) 
NUMS Ets 597 
БЕБЕ 
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ОООО ОООО 


О 


20 


H = H + 1 
CONTINUE 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INTERMEDIATE CALCULATIONS 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ОИЕ ТООТ = ЕТСЕ: 7 DELTHE) + 1 

DELTHETA - DELTHE x DEGTORAD 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INTERPOLATION FOR THETA RESOLUTION 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
DO 30 J = 1, NUMSETS 

DO 40 N = 1, NUMPTOUT 

DETERMINE THETA WHERE DATA POINTS X AND Y ARE TO BE FOUND 
Γ ο = ( М - 1 ) x DELTHETA 

πο ο = 1, NUMPTSIN = 1 
NE І-! 

IE I FQ. 1595 K J) =S B(I,;3J) 

IFCABSCBCI JVD LT. 001) THEN 
WRETAL s PIZZ 

BEBE TFCBCIZJOS SLUT. 0) THEN 
Lam = PI - DATANCDABSCACI, JI7BCI,J))) 

[5 
ῬΑ = DATAN(ABS(A(1,J)ZB(I,J))) 

H 
DETERMINE IF K-TH DATA PAIR FORMS ANGLE GREATER Dn ЭП DEGREES 
BSSUHESSCENTER-OF-COORDINATE SYSTEM AT X = 0 Y 0 
IF(ABS(B(K,J)) .LT. .001) THEN 
THETAK * PI/2 

BESENIBCBCK,J) .LT. 0) THEN 
E = РІ - DATANCDABSCACK,J2/BCK, J22) 

(5 
ШАК = DATANCABSCACK, JO/BCK,J202) 

N 
CHECK TO SEE IF K-TH POINT IS THE DESIRED DATA POINT 
IFCDABSCTHETACN) - THETAK) .LT. .001) THEN 
XFOUND = ACK,J) 


YFOUND = BCK,J) 

RADIUSCN,J) - 2XPIXDSQRT(OXFOUNDXxx2 * YFOUNDX*xX2) 

IFCRADIUSCN,J) «ΙΤ. RADMIND RADMIN * RADIUSCN, J) 

Sues CN, ,J) .GT. RADMAX) RADMAX = RADIUS(ON, J) 
0 

CHECK TO SEE IF I-TH POINT IS THE DESIRED DATA POINT 

ELSE IF(DABS(THETA(N) = THETAI) .LT. .001) THEN 

XFO'ND 7 ACI,J) 

YFOUND = BCI,J) 

RADIUSCN,J) * 2XPIXDSQRT(OXFOUNDXX2 + YFOUNDXX2) 

IFCRADIUSCN,J) .LT. RADMIN) RADMIN = RADIUSCN, J) 

αρ ο ,J) .GT. RADMAX) RADMAX * RADIUSCN, J) 
40 


ες το SEE IF THE DESIRED DATA TUOINT 15 ΒΕΤΗΕΕΝ THE K-TH 
DATA POINT AND ΤΗΕ I-TH DATA POINT 
ELSE TECIRETACN) .GE. THETAK AND. 

ΠΕΙ; ΠΕ. THETAI) THEN 
CHECK FOR VERTICAL LINE 
IF(DABS(A(I,J)-A(K,J)) .LT. .0001) THEN 
XFOUND = ACI,J) 
YFOUND = XFOUND 7 DTANC THETACN) D 
RADIUSCN,J) *» 2XPIXDSQRT(OXFOUNDXx2 
IFCRADIUSCN,J) .LT. RADMIND RADMIN 
IFCRADIUSCN,J) .GT. RADMAX) RADMAX 
GOTO 40 


YFOUNDXx2) 
RADIUSCN, J) 
RADIUSCN, J) 


uou + 
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50 
50 


400 
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ооо о 
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120 


150 


CHECK FOR HORIZONTAL LINE 
Е 5Е ТЕСПАВЪСВСТ, 12 - ος ο τοσο alien 
YFOUND = BCI,J) 
XFOUND = YFOUND * DTANC THETACN) ) 
RADIUSCH,J) - 2XPIXDSQRTCXFOUHDXX2 + YFOUNDXx2) 
IFCRADIUSCH, J) .LT. RADMIN) RADMIN = RADIUSCH, J) 
νο .GT. RADMAX) RADMAX = RADIUS(N,J) 
GOTO 4 
S UE OF LINE CONTAINING I-TH AND K-TH DATA POINTS 

S 
М = ( В(1,2)-В‹К, 2) ) И Ас ;.,.. η 
Y-INTERCEPT OF LINE CONTAINING I-TH AND K-TH DATA POINTS 
BB = B(1,J) = (САСТЫ 
RECALL EQUATION OF ANY LINE IS Y = MX + B 
THE LINE CONTAINING THE IIH AND К=ПН DATA FORNTS INTERSECTS 
THE LINE DRAWN FROM THE ORIGIN THROUGH THE DESIRED DATA POINT 
AT THE DESIRED DATA POINT XFOUND AND YFOUND. THE EQUATION FOR 
THE LINE CONTAINING THE I-TH AND K-TH DATA POINTS MUST THEN 
ALSO WORK AT THE DESIRED DATA POINT. WRITING XFOUND AND 
YFOUND IN TERMS OF THEIR ANGLE THETACN), AND SUBSTITUTING 
YFOUND = BB 7 ( 1 - ( M wx DTAN( THETA(N) ) ) ) 
XFOUHD = YFOUND ж DTANC THETACN) ) 
RADIUSCN,J) * 2XPIXDSQRT(CXFOUNDXX2 + YFOUNDXX2) 
IFCRADIUSCN,J) .LT. RADMIN) RADMIN - RADIUSCN, J) 
IFCRADIUSCN,J) .GT. RADMAX) RADMAX = RADIUSCN, J) 
GOTO 40 
EHDIF 
ώς TO TRY A NEN I-TH AND K-TH DATA PAIR 

5 


DUMMY Ξ 0 
ENDIF 

COHTINUE 
CONTINUE 
CONTINUE 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
OUTPUT TO DATA FILE 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
НЕТТЕ(С, 103) NUMPTOUT, T NUMSETS + 1 

WRITEC8,102) RADMIN, RADMAX 

DO 400 N = 1,NUMPTOUT 

WRITEC8,101) N,CRADIUSON,J), J=1,NUMSETS) 

CONTINUE 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXx 
CLUSEC7) 

CLOSECS8) 

FORMAT (A) 

FORMATCI5,5€E149.6)) 

FORMATC2CE14.6)) 

ЕОКМАТ( 2(Т15)) 

ЕОВМАТС(////////7/7777// БА, ТИМЕГСОӘМЕ ТО ЕМСАР TS, 


1'THIS PROGRAM ACCEPTS PROPERLY FORMATTED INPUT DATA!,Z5X, 
2'FROM CURVE DIGITIZER AND CONVERTS ЕТ ІП FORM 7705902 
S'HHICH CAN BE USED By EMCAD 47 5x. 

G////5X, "PLEASE PRESS ENTER TOMGCONTINUES 7257) 


FORMAT(Z^^///44V4F4444444444445X, ' INPUT FILENAME',7777775X, 


L'TINPUTFN IS THE INPUT FILENAME OF THE FILE CONTAINING',Z5X, 
2'THE DATA OF INTEREST WHICH WILL BE INTERPOLATED ON AND'!,75xX, 
3'CONVERTED TO A FORM WHICH CAN SUBSEQUENTLY BE USED BY',/5X, 
G'EMCAD.',///5X, 

5'PLEASE ENTER THE NAME OF THE INPUT EILE THE EXTENSION А 
6'OF THE FILENAME MUST BE ΡΕ IT E TU EAE ΠΤ. 


FORMAT(Z⁄Z κ αι, OUTPUT FILENAME',Z/IIIZII5K, 


1 "ОУТРИТЕМ IS THE OUTPUT ΕΕ ΕΙΠΕ OE THE TILE CORI να. p 
2'THE DATA WHICH HAS BEEN INTERPOLATED АМО СОММЕКТЕО ТО А!,/5Х, 
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ος πο αι. BE USED BY EMCAD.'! ,ZZZ5X, 
ο ο PEI τι ΤΗΕ ΠΤ THE OUTPUT FILE. THE EXTENSION',^/5X, 
EE ECDHESPEPUBENAMESMUSTOBESTNCEUBED I.E. FILENMAME.DAT',77/5X) 

IOO EJ HAT(ZZ2ZZZZ/ZZ/ZZ/Z/7Z7ZZ/ZZZ75X, 'RESOLUTION!,ZZZZ/ZZZ5X, 
ον ιτ USER eae UT VALUE OF TRE DESIRED THETA',75x, 
МЕСПЕРРРОМФИПКПЕСРЕЕ 2 ZZZ5X; 
ο ο ο εΙτε τος DESIRED DELTA THETA VALUE IN DEGREES.',7//5X) 
EHD ` 
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APPENDIX F 


CRVDIGIN USER'S MANUAL 


Welcome to CRVDIGIN. This program is designed to 
accept properly formatted input data, and to convert it toa 
form which can be used by the West Coast Consultants 
commercial shareware product CURVE DIGITIZER (Appendix M). 
CURVE DIGITIZER is a computer aided design (CAD) package 
which will allow user manipulation of the data. Upon 
completion of using CURVE DIGITIZER, the user may employ the 
program EMCADIN to transform the data into a form which can 


be used by the scattering calculation program EMCAD. 


A. BACKGROUND 

EMSCAT was written by Dr. M. A. Morgan while at the 
University of California at Berkeley. Extensive interactive 
computer aided design type feature have been added to 
EMSCAT, which consider scattering from inhomogeneous 
axisymmetric objects. These computer aided design features, 
driving the change of program name from EMSCAT to EMCAD, 
allow structures composed of multiple dielectric materials 
to be graphically constructed or modified by the user, and 
then evaluated by EMCAD for the scattering. CRVDIGIN was 
written in the programming language of Ryan-McFarland 
FORTRAN (Appendix M). Any questions concerning formats, 
changes to the program, etc. should be referred to the 


128 


respective user's manuals, Dr. M. A. Morgan or LT. E. M. 


Connolly (Appendix M). 


B. COMPUTER SETUP 
CRVDIGIN was not written for any specific computer 
configuration, and as such, should run on most IBM 


compatible machines. 


C. PROGRAM INPUT AND OUTPUT 

CRVDIGIN obtains its input from a previously generated 
data file. The input data file for CRVDIGIN must be located 
in the same directory/subdirectory as CRVDIGIN. This data 
file may be produced by one of two methods. The first 
method consists of special case computer programs which 
generate data files describing the penetrable body of 
revolution for the specific geometrically shaped object 
being considered. A few examples of these special case 
programs are SPHERE, CONE and CYLINDER listed in Appendices 
J, K and L respectively. Each of these programs generates 
two descriptive data files which can be separated into 
distinct areas of information. The first area of 
information is the data which conveys the material 
composition of each layer of the structure. This data is of 
no concern at this particular time, however, it will be 
required before running  EMCAD. The second area of 
information, which is the data required as input for 


CRVDIGIN, is the data which defines the structural shape of 
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each layer of the penetrable body. This data is referred to 
as "STRUCTUR".DAT, and is shown in the flow diagram of 
Figure  F-1. The second method of input data file 
REE consists of user generation through direct 
keyboard input. 

As shown in Figure F-1, the output of CRVDIGIN, shown 
as "CURVEIN".DAT, is the input to the CAD package CURVE 
DIGITIZER. For CRVDIGIN and CURVE DIGITIZER to be able to 
accept input data from a variety of sources, it was 
necessary to establish and adhere to a predefined data 
format. The special case computer programs were written to 
incorporate this format. In order for the second method to 
be utilized successfully, as well as for other special case 
programs to be written, it is imperative that the user 
understand the need for the particular format being used, as 
well as the format itself. The input format for CRVDIGIN is 
relatively uncomplicated. An example of a portion of an 
input data file for CRVDIGIN is shown in Figure F-2. The 
first line of data contains two numbers in integer format. 
The first integer represents the number of data points used 
to describe any layer. The second integer lists the total 
number of layers + 1. The '+ 1' was used to correct some 
minor discrepancies during the design of the code. The data 
defining the physical structure of the boundary for each 
layer is written in radius versus theta format. The value 


of the angle theta associated with each radius is determined 
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ΓΗ 


361 


0.250300E-*01 


4 


0.382760Еғ02 


1 0.188496Е+02 0.251327Е+02 0.314159Е+02 
2 0.176975Е+02 0.241175Е+02 0.304952Е+92 
3 0.166793Е+02 0.231827Е+02 Я. 296292Е+92 
4 0.157731Е+02 0.223194Е+02 Ө. 288130Е+02 
5 0.149613Е+02 0.215196Е+02 Ө. 280428Е+02 
6 0.142301Е+02 0. 207767Е+02 0. 273146Е+02 
7 0.135680Еғ02 0.200848Е+02 0.266253Е+02 
8 0.129657Е+02 0.194390Е+02 Ө. 259719Е+02 
8 0.124155Е+02 0.188348Е+02 0. 253516Е+02 
10 0.119109Е+02 0.182684Е+02 0. 247621Е+02 
11 0.114466Е+02 0.177363Е+02 0.242013Е+02 
12 0.110180Е%02 0.172357Е+02 0. 236670Е+02 
13 0.106211Е+02 0.167638Е+02 0. 231575Е+02 
14 0.102525Е%02 Ө0.183182Е+02 0. 226712Е+02 
15 0.990939Е+81 0.158969Е+02 0. 222066Е+02 
16 0.958921Е+01 0.154980Е+02 0. 217622Е+02 
17 0.928975Е+01 0.151197Е+02 0.213369Е-02 
18 0.900910Е+01 0.147605Е+02 0. 209294Е+02 
19 0.874555Е+01 0.144191Е+02 0.205387Е+02 
20 0.849761Е+01 0.140941Е+02 Ө. 201838Е+02 
21 0.826396Е+01 0.137845Е+02 (0. 198040Е+02 
22 0.804341Е+01 0.134893Е+02 0.194582Е+02 
23: 0.783490E*01 0.132074Е+02 0.191256Е+02 
24 0.763750Е+01 Ө0.129380Е+02 0.188057Е+02 
25 0.745035Ез01 0.126803Е+02 0.184878Еғ02 
26 0.727270Е-01 0.124336Е+02 0.182009Е-02 
27  0.710385E«01 0.121972Е+02 0.179149Е+02 
28 0.684318Е-01 0.119706Е+02 0.176390Е+02 
29 0.679012Е+01 0.117531Е+02 0.173728Е+02 
30 0.664416Е+01 0.115442Е+02 90.171158Е-02 
FIG. F-2. PORTION GF SAMPLE 


INPUT DATA FOR CRVDIGIN. EXE 
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by starting at theta equal to 0 degrees and progressing to 
theta equals 180 degrees in step sizes determined by the 
number of points available. The program expects that the 
values EE the radius for each layer have been entered 
starting with the layer closest to the center of the 
object, then advancing outward from the center until all 
layers have been accounted for. The second line of data 
contains two real values. The first real value is the 
minimum radius value present, and the second value is the 
maximum radius value present. Following the description for 
layer entry above, this implies that the minimum radius 


value should come from the first layer, and the maximum 


radius value should come from the last layer. Figure F-3 
demonstrates the indexing of layers. The next A lines of 
data consists of B columns of data. The value of A is the 


number of points, and the value of B is the number of layers 
+ 1. The first column value is an index value starting at 1 
and incrementing by 1. The Second third, 2.2. «columns 
contain values of the radius associated with theta = (index 
- 1) * deltatheta for the first layer, the second layer 
For each successive row the index changes, the value of 
theta changes and the radii associated with that given value 
of theta are listed. 

The output of CRVDIGIN, used as an input to CURVE 
DIGITIZER, is just as simple in appearance as the input to 


CRVDIGIN, yet much more difficult in theory. The 
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explanation of this difficulty actually is embedded in the 
operation of CURVE DIGITIZER. Although these comments may 
apply to other CAD packages, they are specifically directed 


towards CURVE DIGITIZER, referred to henceforth for brevity 


as CD. Computer aided graphic design deals in the basic 
business of handling points. Any graphic figure can be 
represented as some combination of points. Lines can be 


formed by connecting two points. Arcs can be represented by 
placing more points closer together and by approximating the 
curvature with sufficiently small line segments between 
points. The required input-output format for CD is 
utilized to insure that adequate information is present to 
direct CD how to connect the points. A similar process may 
be that of the simple children's game connect-the-dots. A 
series of dots exists which define a picture, if connected 
properly. As more and more dots are used, more complex 
pictures with better resolution can be constructed. The 
fundamental assumption employed by CD is that each point 
should be connected to the next point listed in the data 
file unless some code is present to indicate otherwise. For 
simplicity and unrestricted accessibility, CD input and 
output data are represented in terms of XY rectangular 
coordinates. CD flags the break between two data points by 
what shall be referred to as a separation code. Since two 
coordinates, X and Y, are required to establish a data 


point, two separation codes, one for the X data and one for 


13 5 


the Y data are required to separate two distinct data 
points. The separation code employed is of the form 999990. 
The last two digits may change to reflect scaling, color, or 
any of SET of other things. It was not possible to 
determine the specifics from the user manual. In the 
Specific case of CRVDIGIN output data, 999990 2999990 codes 
are used to distinguish the structural data describing one 
layer from that describing another layer. Other than the 
insertion of separation codes at the appropriate locations, 
CRVDIGIN performs the simple mathematical task of converting 
the input data generated in radius-theta format to output 


data represented in XY format. 


D. PROGRAM EXECUTION 
Once the correct disk or directory has been selected, 
CRVDIGIN can be invoked by typing 
CRVDIGIN 
at the DOS prompt. The program will load itself and the 


user will see the following screens: 


* Welcome to CRVDIGIN (Figure F-4) 


This screen presents some brief background information. 
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WELCOME TO CRUDIGIN 


THIS PROGRAM ACCEPTS PROPERLY FORMATTED INPUT DATA 

AND CONVERTS IT 10 A FORM WHICH CAN BE USED BY THE 
COMMERCIAL SHAREWARE PRODUCT CURVE DIGITIZER. 

CURVE DIGITIZER I$ A CAD PACKAGE WHICH WILL ALLOW USER 
MANIPULATION OF THE DATA. UPON COMPLETION OF 

CURVE DIGITIZER, THE USER MAY UTILIZE THE PROGRAM 
EMCADIN 10 TRANSFORM THE DATA INTO A FORM WHICH CAN 

BE USED BY THE PROGRAM EMCAD WRITTEN BY DR. M. A. MORGAN. 


PLEASE PRESS ENTER TO CONTINUE. 


ЕНІ 5-2 webLUME Si 
CRVDIGIN MENU 
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* Input Fllename (Figure F-5) 


This screen prompts the user for the name of the file 
μη the input data for CRVDIGIN. The file may have 
any extension that the user desires, however, the extension 
must be entered by the user as part of the filename. The 


input data file must be in the same directory as CRVDIGIN. 
* Output Filename (Figure F-6) 


This screen prompts the user for the name of the file that 
Will hold the output data from CRVDIGIN. Because of the 
input requirements of CURVE DIGITIZER, the extension of this 
file must be .dat, and must be included as part of the 


output filename. 
* Resolution (Flgure F-7) 


CRVDIGIN uses linear interpolation to provide the capability 
to generate output data with better resolution than the 
input data. This screen prompts the user for the desired 


output data resolution, in the units of degrees. 


E. | PROGRAM OUTPUT 

A portion of an example of CRVDIGIN output is shown in 
Figure F-8. Notice the 999990 999990 codes separating the 
different layers of data. 
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INPUT FILENAME 


INPUTFN 1S THE INPUT FILENAME OF THE FILE CONTAINING 
THE DATA OF INTEREST WHICH WILL BE INTERPOLATED ON AND 
CONVERTED TO A FORM WHICH CAN SUBSEQUENTLY BE USED BY 
THE CURVE DIGITIZER PROGRAM. 


PLEASE ENTER THE NAME OF THE INPUT FILE. THE EXTENSION 
OF THE FILENAME MUST BE INCLUDED, I.E. FILENAME DAT 


cee oe ein ВТО КОМ | 
из i ο ENAME 


ШЕ с 


OUTPUT FILENAME 


OUTPUTFN IS THE OUTPUT FILENAME OF THE FILE CONTAINING 
THE DATA WHICH HAS BEEN INTERPOLATED AND CONVERTED 10 A 
FORM WHICH WILL BE USED BY THE CURVE DIGITIZER PROGRAM. 


PLEASE ENTER THE NAME OF THE OUTPUT FILE. THE EXTENSION 
OF THE FILENAME MUST BE .DAT, I.E. FILENAME. DAT 


FIG. TFG CRVBIGIN @RUMEaS 
FOR ОШТРО ЕЛШЕ НАШЕ 
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RESOLUTION 


DELTHE IS THE USER INPUT VALUE OF THE DESIRED THETA 
RESOLUTION IN DEGREES. 


PLEASE ENTER THE DESIRED DELTA THETA VALUE IN DEGREES, 


л ta ΕΠΊΠΤΕΙΝ PROMPT 
NS MES ΠΓΈΞΕΓΟΙ ΓΕΝ 


141 


4. 6939 99848 -18. 84999690 
4. 99180990 -18. 8495969390 
4.00664711 -18.84958080 
3.66403151 -18. 8499586990 
3.32372083 -18.8495883 
2.98551798 -18.84953120 
2.64817397 -18.84854830 
2.31447673 -18.84853880 
1.98121083 -18.84956740 
1.64915908 -18.84959410 
1.31813526 -18.84956740 
0.98790431 -18.84953120 
0.65828347 -18.84950830 
0.32306568 -18.84952740 
4.32930647Е-05 -18.84959980 
999990 998890 
0.00000000Е-01 25.13270000 
ϱ. 4045993559 23.17916870 
0.75102264 21.50648070 
1.85115640 20.05727580 
1.31384826 18.78891840 
1.54581928 17 .66880800 
1.75229287 16.67196666 
1.93734324 19. 7728407160 
2.104 29401 14. 8729597109 
2.25564241 14.2415719 
2.32 19. 9908286 
-IG. —F=S PD Fie nu R NS 


ӘЦТРГЦИТКПАТА KOT ШЕ 
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APPENDIX G 
CRVDIGIN SOURCE CODE 


PROGRAM CRVDIGIN 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


МЕШООМЕТІП СКУШІСЕН WRITTEN BY LT E. M. CONNOLLY. 

CRVDIGIN IS A PROGRAM DESIGHED TO ACCEPT PROPERLY FORMATTED 
INPUT DATA, AND CONVERT IT TO A FORM WHICH CAN BE USED BY 

THE COMMERCIAL SHAREWARE PRODUCT CURVE DIGITIZER. 

CURVE DIGITIZER IS A CAD PACKAGE WHICH WILL ALLOW USER 

πο ο «ΤΟ OF THE DATA. UPON COMPLETION OF CURVE DIGITIZER, 
THE USER MAY UTILIZE THE PROGRAM EMCADIN, ALSO WRITTEN BY 

LT E. M. CONNOLLY TO TRANSFORM THE DATA INTO A FORM WHICH CAN 
BE USED BY THE PROGRAM EMCAD WRITTEN BY DR. М. A. MORGAN. 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


VARIABLE DEFINITIONS 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


INPUTFN ITHE INPUT S FILENAME OF THE FILE CONTAINING 
THE DATA OF INTEREST WHICH WILL BE INTERPOLATED 
ON AND CONVERTED TO A FORM WHICH CAN 
SUBSEQUENT EY SBE USED OBEY SIRE CURVE DIGITIZER 
PROGRAM 

OUTPUTFN THE OUTPUT FILENAME OF THE FILE CONTAINING 
THE DATA WHICH HAS BEEN INTERPOLATED AND 
CONVERTED TTO A FORM WHHICHENILL BE USED BY 
THE CURVE DIGITIZER PROGRAM 


I, x: K, INTEGER COUNTERS 

|, 

NUMBRONS THE NUMBER OF RONS CONTAINED IN THE DATA 
MATRIX IN THE FILE INPUTFN 

NUMBCOLS THE NUMBER OF COLUMNS CONTAINED IN THE DATA 
MATRIX IN THE РШЕ ЕНЕ ОЛЕ М 

SEPRATN A FLAG USED BY CURSE DIGITIZER TO DISTINGUISH 
BETHEEN DATA POINTS OR SETS OF DATA 

RADARRAY A WORKING ARRAY USED TO HOLD ALL THE VALUES 
OF THE ARRAY SPECIFIED ABOVE BY NUMBROWS AND 
NUMBCOLS 

MINRAD THE MINIMUM RADIUS VALUE, NOT USED BY THIS 
PROGRAM BUT REQUIRED FOR EMCAD 

MAXRAD THE MAXIMUM RADIUS VALUE, NOT USED BY THIS 
PROGRAM BUT REQUIRED FOR EMCAD 

RADIUS A WORKING VARIABLE USED TO HOLD THE CURRENT 
VALUE OF THE RADIUS OBTAINED FROM RADARRAY 

ЕТ A CONSTANT 

DEL THE THE USER INPUT VALUE OF THE DESIRED THETA 
RESOLUTION IN DEGREES 

DELTHETA THE USER INF UT VALUE ОРАТ МЕ DESIRED THETA 
RESOLUTION IN RADIANS 

ΤΗΕΤΑΟΙ Ρ A WORKING VARIABLE USED TO HOLD THE CURRENT 
VALUE OF THETA OBTAINED FROM RADARRAY 

THETA AN ARRAY HOLDING THE DESIRED VALUES OF THETA 
AS DETERMINED BY DELIHE 

DEGTORAD A CONVERSION FACTOR FROM DEGREES TO RADIANS 

X THE VALUE OF X OBTAIHED BY CONVERTIHG THE POLAR 


COORDIHATES RADIUS,THETA TO X,Y COORDINATES 
THE X,Y COORDIHATE SYSTEM IS REQUIRED FOR THE 
DATA TO BE COMPATIBLE HITH THE CURVE DIGITIZER 
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ООО О ООООООО ООО О 


ООО О 


ОООО 


OOOO 


OOOO 


PROGRAM 

Y THE VALUE OF Y OBTAINED BY CONVERTING THE РОВА 
COORDINATES RADIUS, THETA TOX r COORDINATE 
THE X,Y COORDINATE SYSTEM IS REQUIRED ЕОК ΠΕ 
DATA TO BE COMPATIBLE WITH THE CURVE DIGITARE 
PROGRAM 

ΙΕ], ΤΕΕ: IHTERMEDIATE CALCULATION VARIABLES 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
VARIABLE DECLARATION 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CHARACTERx12 INPUTFN, OUTPUTFN, DUMMY 

INTEGER I, К СН 

INTEGER NUMBROMS, NUMBCOLS, NUMPTS, SEPRATN 

REAL RADARRAYC1000,10) 

REAL MINRAD, MAXRAD, RADIUS, PI, DELTHE, DELTHETA, THETAOLD 

REAL THETACI0005, DEGTORAD, X, Y 

REAL INTERI TINTER 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INITIAL VALUES 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
Γι SIGS) 

DEGTORAD - PI / 180. 

ЗЕРКАТМ - 989990 


X X XXXXXXXXXXDOOGGGOGGOGOCOXGXX X XX XXX X YX X XXX XX XX XX XX XXXx Xxxxxxxxxxxx 
INPUT DATA OF INTERE T 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITECX,110) 

READCX,100) DUMMY 

WRITECX,T29) 

READCX,100) INPUTFN 

МЕІТЕ(Ж%,150) 

READ(x,100) OUTPUTFN 

ИРІТЕСЖ, 140) 

READCX,x) DELTHE 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
READ DATA FILE 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
OPEH(7,FILE = INMEUDEN) 

OPENC8,FILE = OUTPUTFN) 

READC7,X) NUMBROWS, NUMBCOLS 

READC7,X) MINRAD, MAXRAD 

READC7,%) CCRADARRAYCI,J), J=1,NUMBCOLS), I = 1,NUMBROWS) 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INTERMEDIATE CALCULATIONS 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
NUMPTS = INTCISUL. ΕΙΠΕ, pl 

DELTAOLD = (180. / (МОМВКОН5 - 12) Xx DEGTORAD 

DELTHETA = DELTHE ж DEGTORAD 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INTERPOLATION AHD TRANSFORMATION 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
J = 1 

DO 10 K = 2, NUMBCOLS 

DO 20 N = 1, NUMPTS 

THETACN) = (N-1) x DELTHETA 

DO 30 I = 1, NUMBRONS 

THETAOLD = CRADARRAYCI,J) =~ 1) ® DELTAOLD 

IFCABSCTHETACN) - THETAOLD) .LT. .00001) THEN 

RADIUS = RADARRAYCI,K) 

X = RADIUS * SINC THETACN) } 
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ΩΩ 


Y = RADIUS жх COSC ТНЕТАСМ) ) 

HRITECS,; κο αχ 

COTO 20 

ЕМЕ πρ ΠΙΕΙ σε ο ος - DELTAOLD) .AND. 
l ТЕЛЕТА СНО T LCE OUNETAUDDOCTHEN 

кее 

ІМТЕКІ = RADARRAY(L,K) - RADARRAY(I,K) 

PTER? =  CTHETIACN N = THETAOLD) Z (=DELTAOLD) 
RADIUS = INTER1 xX INTER2 + RADARRAY(I,K) 

FE IO E JIOPMALIZINOG FACTOR OF 2xPI GENERATED BY SPHERE, CONE, 
ЕШ, FACTOR IS RESTORED IN EMCADIN 

PE RADIUS X SINC THETACND) ) x 1 Z ( 2 x PI ) 

ο E DIUS X COSC THETACND ο X* 1w «2 * PI. 
WRITECS,X) X,Y 


GOTO 20 

ELSE 

DUMMY BRANCH, TRY NEXT VALUE OF THETAOLD 
RADIUS = 0 

ENDIF 


30 CONTINUE 
20 CONTINUE 

WRITEC8,X) SEPRATN,SEPRATN 
10 CONTINUE 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 

100 ΕΟΚΜΑΤ{(Α) 

110 ΕΟΒΝΜΑΤ(2,,,7 77 "Χ, 'HELCOME TO CRVDIGIN',/7^75X, 
I'THIS PROGRAM ACCEPTS PROPERLY FORMATTED INPUT DATA',/⁄/5X, 
ANT CONVERTS IT TO A FORM WHICH CAN BE USED BY THE !,/5Х, 

Ss GOMMERGCIAL SHAREWARE PRODUCT CURVE DIGITIZER.', 75x, 

G*CURVE DIGITIZER IS A CAD PACKAGE WHICH WILL ALLOW USER', 75x, 
ο ο ο. τον OF THE DATA. UPON COMPLETION OF',75x, 

6'CURVE DIGITIZER, IHE USER MAY UTILIZE TRE PROGRAM ',75xX, 
7'EMCADIN TO TRANSFORM THE DATA INTO A FORM WHICH САМ!,/5Х, 
8'BE USED BY THE PROGRAM EMCAD WRITTEN BY DR. M. A. MORGAN.', 
775X; ee ASE PRESS. ENTER. TO CONTINUE P! yZZ5X) 

120 FORMATOSSSSSSSSSSSSSASASASA5SX,*INPUT FILENAME! ,7777775X, 
1'INPUTFN IS THE INPUT FILENAME OF THE FILE CONTAINING', 75x, 
ОЛЧЕ DATA OF INTEREST WHICH WILE BE INTERPOLATED ON AND',75xX, 
3' CONVERTED TO A FORM WHICH CAN SUBSEQUENTLY BE USED BY!,75xX, 
TIRE CURVE DIGITIZER PROGRAM.',7775X, 

PL IE SE ENTER THE NAME OF THE INPUT FILE. THE EXTEMSION',/5X, 
CUT THE FILENAME MUST BE INCLUDED, I-E. FILENAME.DAT',//5X) 

130 БОВМАТ(//////////////;////7/7/5Х, ООТРОТ ҒІПЕНАМЕ?,//////5Х, 
УОИ ОЕ Н 15 THE OUTPUT FILENAME OF THE FILE CONTAINING',/5X, 
2'THE DATA WHICH HAS BEEN INTERPOLATED AND CONVERTED TO A',^/5X,. 
Se FORM WHICH ATER BE USED BY THEWCURVE DIGITIZER PROGRAM.',7775x, 
ИРВЕА5Е ENTER THE NAME OF THE OUTPUT FILE. THE EXTENSION!,75x, 
ΜΓ THE FILENAME MUST BE .DAT, I.E. FILENAME. DAT',775xX) 

160 БОВМАТС(УИ/,////уу//////,,////5Х, "ВЕ5ОТУТІОМ!,//7////5Х, 
IGDEDBHESTSCTHE USER ТНРОТ VALUE OF ΤΗΕ DESIRED THETA'!' /⁄|kls,5Xx, 
EBERESOEUTTON IN DEGREES.',7775X, 

PI EAS SE ENTER THE DESIRED DELTA THETA VALUE IN DEGREES.',ZZ5X) 
CLOSE(7) 

CLOSEC8) 

END 
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APPENDIX H 


MCADGRAF USER'S MANUAL 


Welcome to MCADGRAF. This program was written to be 
used in conjunction with the FORTRAN program entitled EMCAD. 
MCADGRAF adds graphics features to support the analysis of 


data generated and formatted by EMCAD. 


A. BACKGROUND 

MCADGRAF was developed and written by LT. E. M. 
Connolly using the programming language of Ryan-McFarland 
FORTRAN (Appendix M). Microcompatibles GRAFMATIC (Appendix 
M), and later Jewell Technologies GRAFPLUS (Appendix M) were 
chosen as the commercially available screen graphics and 
hard copy support for EMCAD output graphics. Any questions 
concerning colors, formats, changes to the program, etc. 
should be referred to the respective user's manuals, Dr. M. 


A. Morgan or LT. E. M. Connolly (Appendix M). 


B. COMPUTER SETUP 
MCADGRAF was written for the following computer 
configuration: 
* IBM AT or compatible, including 80386 based 
compatibles 
* hard disk 


* EGA graphics card 
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ж EGA monitor 

* math coprocessor - not required 

* memory - 512K RAM minimum 
MCADGRAF E be placed in any directory, at the user's 
discretion, however, it is highly recommended that MCADGRAF 
and EMCAD be placed in the same directory.  MCADGRAF input 
is taken directly from EMCAD output. This output is 
required to be in the same directory/sub-directory as 
MCADGRAF, or on the same floppy disk if using a floppy 
system instead of a hard disk. The only other requirement 
for initial setup deals with printer initialization. In 
order to obtain a printer graphical output, GRAFPLUS must be 
installed on the computer, configured for the specific 
printer of use, prior to running MCADGRAF. An EPSON LX- 
800 series printer was used when writing this code. This 
printer, as with any other printer to be used, was 
Configured by inserting two lines in the AUTOEXEC.BAT file. 
The first of these lines was a path command to tell the 
computer where to find the appropriate files for GRAFPLUS. 
Assuming all the files for GRAFPLUS are contained in a 
directory on the C disk entitled Grafplus, an appropriate 


path entry would be 


path c:\/c:\Grafplus 


This setup will also work with other IBM graphics compatible 
printers such as the Okidata ML - 192/193. The second of 
these lines inserted actually established the initialization 
values fof the printer to be used. The second line must be 
inserted anywhere in the  AUTOEXEC.BAT after the path 


command. The printer initialization used for the LX-800 was 
GRAFPLUS = #115 


The equals sign (=) directs the printer to print the 
graphics horizontally. This means that the top of the image 
is at the top left corner of the page. For vertical 
printing, the object appears to be rotated with the top of 
the picture at the bottom left corner of the page. The 
pound sign (#) directs the printer to print the text screens 
as graphics screens. The I conveys the information that the 
screen mode 1s EGA - 640x350 resolution, with 16 colors 
available. The number 15 is the number assigned by GRAFPLUS 
to correspond with the EPSON  LX-800. Other options 
available are the capabilities to print in color and to 
reverse printing. All of these options are discussed in 
much greater detail on pages 3 through 8 of the GRAFPLUS 
user's manual. As a note of caution, some GRAFPLUS printer 
definitions are defined such that "graying" is used on a 
black and white print-out to represent different colors. On 


the Lx-800, the "graying" effect reduced printer clarity and 
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resolution dramatically. This problem was rectified by 
running the EDITGRAF.COM program contained in GRAFPLUS, 
answering all the questions without changes until the 
Eton. "Do you wish to print in Black & White only? 
(disable gray scale) " appeared on the screen. Changing the 
answer from no to yes removed the "graying" and improved the 
resolution. A much more detailed explanation of EDITGRAF 
can be found on pages 9 through 21 of the GRAFPLUS user's 
manual. Once these changes to the printer definitions and 
AUTOEXEC.BAT have been made, it 1s necessary to re-boot the 
computer to initialize these values. This re-boot can be 
accomplished by either turning the machine off, waiting 15 
E20 seconds, and turning it back on, or by hitting the 


Ctrl, Alt and Del keys all at the same time. 


C. PROGRAM INPUT 

As previously mentioned, MCADGRAF obtains its input from 
the output of EMCAD. This EMCAD output is four files in the 
form of "FILENAME".TET, "FILENAME".TEP, "FILENAME".TMT and 
"FILENAME".TMP. The quotation marks indicate that the 
FILENAME is a value that was previously chosen by the user 
while running EMCAD. The four extensions of FILENAME have 
significant meaning. The TET refers to TE incidence, 
F-Theta; the TEP refers to TE incidence, F-Phi; the TMT 
refers to TM incidence, F-Theta; the TMP refers to TM 
incidence, F-Phi. The values for these four files were 
taken directly from a much larger, more comprehensive EMCAD 
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output data file. The scattering calculation output and 


graphics relationship is depicted in Figure H-1. 


D. PROGRAM EXECUTION 
MCADGRAF is a menu driven output graphics program. Once 
the correct disk or directory has been selected, MCADGRAF 


can be invoked by typing 
MCADGRAF 


at the DOS prompt. “he program will load itself, and the 


user will see the following screens: 
* Welcome To MCADGRAF (Figure H-2) 
This screen presents some brief background information. 


ж EMCAD Output Graphics Presentation, Part 1 
(Figure H-3) 

This menu affords the user the opportunity to exit the 

program in the event that the printer initialization has not 


been accomplished. 


* EMCAD Output Graphics Presentation, Part 2 
(Figure H-4) 
This screen explains the means of obtaining an output graph 
on the printer. The user will be presented with four graphs 
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ТУузтнаучо ОМА 


ne 


женененене WELCOME TO. MCADGRAF seeeeeee 


THIS PROGRAM WAS WRITTEN TO BE USED IN CONJUNCTION WITH 

THE FORTRAN PROGRAM ENTITLED EMCAD. MCADGRAF ADDS 

GRAPHING FEATURES TQ SUPPORT THE ANALYSIS OF DATA 

GENERATED AND FORMATTED BY EMCAD. EMCAD WAS WRITTEN BY 

DR. M. A. MORGAN, WITH MODIFICATIONS BY HIMSELF AND 

LT E. M. CONNOLLY. MCADGRAF WAS DEVELOPED AND WRITTEN BY 
LT E. M. CONNOLLY TO HELP FULFILL THESIS REQUIREMENTS WHILE 
UNDER THE GUIDANCE AND THE SUPPORT OF DR. MORGAN,  MCADGRAF 
WAS WRITTEN AND COMPILED USING RYAN MCFARLAND FORTRAN, 
GRAFMATIC, AND GRAFPLUS. ANY QUESTIONS CONCERNING COLORS, 
FORMAT, CHANGES TO THE PROGRAM, ETC... SHOULD BE REFERRED 
TO THE RESPECTIVE USERS MANUALS, DR. M. A. MORGAN, OR 

LT E. M. CONNOLLY. 


PLEASE PRESS ENTER TO CONTINUE. 


FIG. H-2. ΝΕ ΕΕ ТЗ Е ВЕ 


= 2 


EMCAD OUTPUT GRAPHICS PRESENTATION 

NOTES CONCERNING GRAPHICAL HARDCOPY QUTPUT, PART 1 

This program allows the user to obtain printer hardcopies 
of the graphs being viewed. Hardcopy graphic outputs are 
obtained through the use of the program entitled GRAFPLUS. 
Proper utilization of MCADGRAF requires prior setup of 
GRAFPLUS. Directions for proper setup are contained 

in the GRAFPLUS users manual. 

l. CONTINUE WITH THE NEXT SCREEN 

b EXIT THE PROGRAM 


Indicate your selection by the choice of 1 or 2. 


ЕЕЕ O RE MCADTOUTPUT 
В E е SENAT ON PART 1 
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EMCAD OUTPUT GRAPHICS PRESENTATION 
NOTES CONCERNING GRAPHICAL HARDCOPY OUTPUT, PART Z 


The user makes the decision as to whether a printer 
graphical output is desired after viewing the graph. Once 
a graph is drawn on the screen, the user nay request a 
printer hardcopy by typing either a "P" or "p" and pressing 
the "ENTER" key. Selection of any other key vill sinply 
advance to the next graph. Spacing on 8 1⁄2 by 11 paper 


allows up to 2 graphs per page. 


PLEASE PRESS ENTER TO CONTINUE. 


FIG. HA ЕСЕ! 
GRAPHICS PRESIDEKNET SDN S и 
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mom each data run. As each graph is presented on the 
Screen, the user is given the opportunity to obtain a 
printer copy of the screen. The user decides to obtain a 
printer copy by depressing a P or a p, followed by ENTER. 
The reminder line telling the user what to do to produce a 
hardcopy is erased, and the graph is sent to the printer. 
Any other character input, followed by ENTER will simply 
advance to the next screen. Setting up the paper in the 
printer is the sole responsibility of the user for each and 
every graph. If one graph is to placed on every page, it is 
recommended that the paper be advanced a few inches so that 
the graph is printed closer to the center of the page. For 
the LX-800, it is possible to place 2 graphs on every page, 
if the paper is positioned properly in the printer. The 


best advice is to experiment. 


* Data File Name Input Menu (Flgure H-5) 


If the data file to be used was generated by EMCAD, the user 
need only choose option 1. If the file to be used is a 
self-generated file, it must be in the exact format of an 
EMCAD generated file, and it must have an extension of .TET, 
ШЕР, .TMT or .TMP . πε спе τίς 15 пос in the correct 
format, or if the file does not contain the proper 
extension, the user must exit the program and correct the 


problem prior to continuing. 
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EMCAD QUTPUT GRAPHICS PRESENTATION 
DATA FILE NAME INPUT MENU 
If utilizing a file or series of files created bu EMCAD, 
please select option 1. If utilizing any file other than 
one created by EMCAD, please ensure that the file is in 
proper format, and that the filename extension is one of 
the four extensions utilized by this program. These 
extensions are: „tet Дер .tnt and {мр 
If you are unsure as to whether your file is in the proper 
format with the proper extension, please exit the progran 
and verify before continuing. 

1. INPUT THE FILE OF USER CHOICE 

à. EXIT THE PROGRAM 


Indicate your selection bu the choice of 1 or 2. 


FIG. Πτα ШАЙ ЕЕ E 
INPUT MENU 
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* Data File Name Input (Figure H-6) 


This menu prompts the user for the filename of the input 
data file. Since the extension is limited to four cases, 
and has previously been selected in the Data File Type 
Selection Menu, no extension should be added to the filename 


to be typed in. 
* Data File Type Selection Menu (Figure H-7) 


As previously discussed, MCADGRAF requires that the four 
EMCAD output data files be located in the same directory. 
This menu gives the user the chance to choose which of the 
four files the data will come from. 
This menu also affords the user the opportunity to change 
the default setting of the aspect ratio. The menu for this 
change is shown in Figure H-8. Aspect ratio is a term that 
relates the number of available pixel columns on the screen 
to the number of available pixel rows on a screen. Since 
most computer screens are not square, this representation 
becomes very important. In simple terms, it deals with 
making a square look like a square, and not like a 
rectangle, or making a circle look like a circle, and not an 
oval. 

GRAFPLUS comes with a program to help determine the 


aspect ratio του any monitor. In the case of the 
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EMCAD OUTPUT GRAPHICS PRESENTATION 

DATA FILE NAME INPUT 

Please input the name of the file containing data that 
you wish to work with, Please input only the file nane 
1.6. do not include the file extension. The file 
extension will be selected in the DATA FILE TYPE 


SELECTION MENU. 


FIG. H-6. DATA FILE NAME ΠΝ 


ШЕН 


EMCAD OUTPUT GRAPHICS PRESENTATION 
DATA FILE TYPE SELECTION MENU 


Please select the type of data file containing the information of 
interest. EMCAD automatically created these four data file type s. 


1. TE INCIDENCE, Е-ТНЕТА 
2. ТЕ INCIDENCE, ЕН! 
3. IM INCIDENCE, — F-THETA 
4. IM INCIDENCE, — F-PHI 
9. CHANGE DATA FILE 

Ы. CHANGE ASPECT RATIO 

7. EXIT THE. PROGRAM 


Indicate your selection bu the choice of 1, 2, 3, 4, 5, 6 or ? 


ИСЕ 7500007 01 ЕСЕ ТҮРЕ 
ЕЕ Иа 


wegen ASPECT RATIO xxxxxxxxxx 


IN ORDER TO AVOID DISTORTION ON THE SCREEN, IT IS NECESSARY 
TO DEFINE THE NUMBER OF ROWS PER INCH AND THE NUMBER OF 
COLUMNS PER INCH. THESE DEFINITIONS CHANGE AS THE 
RESOLUTION OF THE SCREEN, AS WELL AS THE TYPE OF SCREEN 
CHANGE. THIS PROGRAM WAS WRITTEN IN MODE 16 (EGA) USING 

A GB-1 VIDEO BOARD AND A NEC MULTISYNC MONITOR. THE 
OPTIMUM ASPECT RATIO FOR THIS CONFIGURATION WAS .65 


PLEASE INPUT ASPECT RATIO. 


FIG. H-8, А5РЕШ E 
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development of this program, the aspect ratio given by 
GRAFPLUS delivered erroneous results. As such, this screen 
of MCADGRAF allows the user to input and change the aspect 
ratio of I For users unsure of their aspect ratio, 
start with an initial value of 1, and work up or down with 
increments of .1 or .2 until close to the desired results, 


then fine tune as much as needed. 
* Selection of Incidence Angle (Figure H-9) 


One of the input parameters for EMCAD was the number of 
incident fields, and their respective angles. MCADGRAF 
reads the input data file, and presents in this menu the 
number of incidence angles, and their respective angles in 
degrees. It is up to the user to select the incidence angle 


of interest. 

* Selection of Phi Angle (Figure H-10) 
As with incidence angles, phi angles were input parameters 
for EMCAD.  MCADGRAF presents in this menu the number of phi 


angles, and their respective angles in degrees. It is up to 


the user to select the phi angle of interest. 
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EMCAD OUTPUT GRAPHICS PRESENTATION 


SELECTION OF INCIDENCE ANGLE 

The data file you have selected contains data for 5 
incidence angles, Please select your choice from the 
following menu, Note, you may only work with one 
incidence angle at a time. 


INCIDENCE ANGLE = 4.08 

INCIDENCE ANGLE = 45,00 

INCIDENCE ANGLE = — 96.88 

INCIDENCE ANGLE = 135.08 

INCIDENCE ANGLE = 168,00 

RETURN TO THE DATA FILE TYPE SELECTION MENU 
EXIT THE PROGRAM 


—À СГ СП мь C1) TM ма 


INDICATE YOUR SELECTION BY THE CHOICE OF A NUMBER [.E. 1 2 3...) 


FIG. ΠΞ5. ΞΕΕΕΞΕΠΙ ες; 
INCIDENCE ANGLE 
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EMCAD OUTPUT GRAPHICS PRESENTATION 


SELECTION OF PHI ANGLE 

The data file you have selected contains data for 3 

phi angles. Please select your choice from the following 
menu. Note, you may only work with one phi angle at a tine. 


PHI ANGLE = 0.00 

PHI ANGLE = 90,00 

PHI ANGLE = 100.00 

RETURN TO THE DATA FILE TYPE SELECTION MENU 
EXIT THE PROGRAM 


C1 „Ре «АЗ P — 


INDICATE YOUR SELECTION BY THE CHOICE OF A NUMBER I.E.1 Z 3 ...) 


РИБНЕ TID. SELECTION OF 
PHI ANGLE 
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E. PROGRAM OUTPUT 
For each separate data set, MCADGRAF produces four 


different graphs. 


* Magnitude Versus Theta (Figure H-11) 


As can be seen in the first line of the title for this 
graph, as well as the other three graphs, a label specific 
to this data set appears. This label was used as an input 
parameter to EMCAD to assist the user in differentiating 
this set of data from every other set. In the third line, 
the values of the incidence angle and phi angle, as well as 
the type of incidence are listed. Either F-Phi or F-Theta 
appears in the fourth line followed by a short description 
of the type of graph. The line at the bottom referring to 
hardcopy is simply a reminder to the user of what is needed 
to request a hardcopy output. This line is erased before 
the printer graph is produced. This first graph depicts the 
scattering magnitude plotted versus the 180 degrees of 


Theta. 


* Phase Versus Theta (Figure H-12) 


The only change in this graph from the first graph is that 


phase instead of magnitude is being graphed versus Theta. 
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Offset Sphere Validation Run 
Scattered Field Conplex Anplitudes 
TE Incidence = 98. Degrees Phi 5 4. Degrees 


F-Phi Magnitude Versus Theta 

Magnitude 

с 2.01 

C 

à 

| 

Ë 

κ 18 

19 

x 

х 

1 

Πης------. E ee се ШЕГЕН 
8.8 45.0 90,9 135. H 188.8 


Hardcopy —-? Р ог р 


О КОШИ eee AGNI TE VERSUS 
ШЕРА. 


συ 


Conposite Dodu 
scattered Field Conplex Amplitudes 
IM Incidence = 135, Degrees Phi = 108. Degrees 





Phase F-Theta Phase Versus Theta 
108.8; 
98.8... 
5.0 98.8 (35.0 (вв aa 


2000р 


-ó88.0* | 
Hardcopu —? P or p 


εν ο 12 PHA ER ЕСКЕНЕ 
ШЕБЕР 
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ж Magnitude Versus Theta, Polar Plot 
(Flgure H-13) 
This graph represents the same data as the first graph, but 
from a Ae erent viewpoint. The first graph was drawn using 
rectangular coordinates. This graph utilizes polar 
coordinates. This graph type is the one predominately used 


in the industry to illustrate scattering. 


* Magnitude Versus Theta, Polar Plot 
Axis Rotated Minus Ninety Degrees 
(Flgure H-14) 

The only change in this graph from the third graph is the 
rotation of the axis. The limitations of GRAFMATIC made it 
difficult to generate the data, maintain spacing 
restrictions and still display all the data necessary to 
adequately describe the situation. Rotating the axis ninety 
degrees displays the data slightly differently from the 


third graph, giving the user a bit more useful information. 
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Offset Sphere Validation Run 

scattered Field Conplex Amplitudes 

TE Incidence - 98. Degrees Phi 5 4. Degrees 
Polar Plot, F-Phi Magnitude Versus Theta 





Ihetaz98 
18 1.0 ¿.B 


Thetaz188 Hardcopy ---> P or p 


FIG. H=-is. АШК НЕШЕ эзше 
THETA, EOL ЕЕЕ 
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Conposite Bodu 

Scattered Field Conplex ñnplitudes 

IM Incidence = 135. Degrees Phi = 188, Degrees 
Polar Plot, F-Theta Magnitude Versus Theta 

Axis Rotated Minus Ninety Degrees 
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Hardcopy ---> P or p 


US. (|2122 C ODE VERSUS 
ШЕП ШЕ ИШЕ 5 
ЕЕЕ ТЕЗА Νες DEGREES 


ne 


ООООООООООООООООООООООООООООООООООООООООООООООООООООООЭОООО 


APPENDIX I 
MCADGRAF SOURCES GODE 


PROGRAM MCADGRAF.FOR 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


THIS PROGRAM WAS WRITTEN TO BE USED IN CONJUNCTIUN WITH THE 
FORTRAN PROGRAM ENTITLED EMCAD. MCADGRAF ADDS GRAPHING FEATURES 
TO SUPPORT THE ANALYSIS OF DATA GENERATED AND FORMATTED BY EMCAD. 
EMCAD WAS WRITTEN BY DR. M. A. MORGAN, WITH SUBSEQUENT UPDATE BY 
HIMSELF AND LT E. M. CONNOLLY WHILE AT THE NAVAL POSTGRADUATE 
SCHOOL (MODS MAR 87 - МАК 88). MCADGRAF WAS DEVELOPED AND WRITTEN 
BY LT E. M. CONNOLLY TO HELP FULFILL THESIS REQUIREMENTS PI EE 
UNDER THE GUIDANCE AND THE SUPPORT OF DR. MORGAN. MCADGRAF WAS 
WRITTEN AND COMPILED USING RYAN MCFARLAND FORTRAN, GRAFMATIC, AND 
GRAFPLUS. ANY QUESTIONS CONCERNING COLORS, FORMAT, CHANGES TO THE 
PROGRAM, ETC... SHOULD BE REFERRED το τι КЕЗГЕ ЕК КИШЕН 
MANUALS, DR. M. A. MORGAN, OR LT E. M. CONNOLLY. 


A NOTE OF CAUTION, TO OBTAIN CORRECT PRINTOUTS FROM MCADGRAF, 
GRAFPLUS MUST BE PROPERLY CONFIGURED TO SUPPORT YOUR HARDWARE. 
QUESTIONS CONCERNING THE PROPER HARDWARE SETUP SHOULD BE REFERRED 
TO THE GRAFPLUS USER'S MANUAL. 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
X X X X X K K K K KK K KKK KK KKK K KKK K KXK KK KK KK K K KKK KKK KK XK XK XK KK KKK KKK KKK KX X X X X X X X 
PROGRAM STATEMENT DEFINITIONS 

ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


QSMODE (mode) 


Clears the screen and sets the mode of the CRT 
display. May be used in either the text or the 


graphics mode. The desired value of mode is 
chosen from the following table: 

MODE RESOLUTION (CxR) COLOR T/G 
0 40 Х 25 B/W T 

1 40 X 25 С T 

2 80 X 25 B/W T 

3 80 X 25 (2 Т 

4 320 X 200 С ο 

5 320 X 200 B/W G 

6 640 X 200 B/W G 

y 80 X 25 B/W TCM) 
13x 320 X 200 C G 
14x 660 X 200 (с б 
15x 660 X 350 50HZ MONO 
16x 660 X 350 4/16 COLORS 
X ----» Requires EGA card 


QPLOT © сос; jcol2, jrowl, jrow2, xmin, xmax, ymin, ymax, 
xorg, yorg, l1opt, yoverx, aspect) 
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5coll 
з) о 
jrowl 
2гом2 
xmin 
xmax 
ymin 


ymax 


xorg 
yorg 


іор% 
yoverx 


aspect 


QXAXIS ( xst, 


xst 


Sn 
xmajor 


minor 


label 


ndec 


QYAXIS ( vst, 


Plot initialization routine: sets plot window, 
corresponding x and y ranges, and scaling 
factors. 


left hand pixel column defining plot region 
right hand pixel column defining plot region 
lower pixel row defining plot region 

upper pixel row defining plot region 

value of xX variable assigned to correspond to 
column 3coll 

value of x variable assigned to correspond to 
column 3col2 

value of y variable assigned to correspond to 
row jrowl 

If iopt = 0, then ymax is the value of the y 
Variable assigned to row = jrow2. If iopt = l, 
then on input ymax is the assumed maximum value 
of the y variable required by the plot. 
vertical Cy) axis intersects the horizontal (x) 
axis at x = xorg 

horizontal (x) axis intersects the wertical (v) 
axis at y = yorg 

sets scaling method option, see ymax 

number of y unit changes for a unit change in 
value of x, ie.. the change in the y coordinate 
per inch divided by the change in the x 
coordinate per inch 

aspect ratio for the user's CRT, for NEC 
Multisync, aspect = 0.6 for mode < 16 


xfin, xmajor, minor, label, ndec) 


Draws the x axis along the horizontal direction, 
including tic marks and tic mark labels 


value of x at the start of the x axis; 

If a logarithmic scale is being used, then 

xst and xfin must be one of the following 

with xfin > xst : 

0.001, 0.01, 1.0, 10.0, 100.0, 1000.0, 10000.0 
value of x at the end of the x axis 

distance between the major (large) tic marks; 
If xmajor is positive/negative, then the tic 
marks will be placed above/below the x axis 
respectively. If xmajor is zero, then no tic 
marks are inserted 

integer flag for use of minor tic marks or 
logarithmic scale; If minor = 0, then no tic 
marks are inserted. If minor = +l, then minor 
tic marks are inserted halfway between the 
major tic marks, but are not labeled. If 
minor = -l, then a logarithmic set of tic marks 
are inserted. 

integer flag for auto labeling of tic marks; 
If label = 0, then there is no auto labeling. 
If label = 1, then the major tic marks are 
automatically supplied with a label. If label 
= =], then the origin is labeled as well. 
number of decimal places for auto labeling; 
The format of the label is Fl10.n with the 
decimal point lined up with the major tic marks. 


yfin, ymajor, minor, label, ndec ) 


171 


ООООООООООООООООООООООООООООООООООООООООООО ООО«ООООООООООООООООООООЮО 


yst 
yfin 
yma J or 
minor 
label 


ndec 


Creates the y axis for a plot. 


value of y at the start of the y axis 

value of y at the end of the y axis 
distance between major (large) tic marks on 
the y axis 

integer flag for use of minor tic marks or 
logarithmic scale, see QXAXIS 

integer flag for auto labeling of tic marks, 
see QXAXIS 

number of decimal places for auto labeling, 
see QXAXIS, 0 <= ndec <= 3 


QPREG ( ireg, ival) 


ireg 
ival 


QPTXT (С nchar, 


nchar 
char 


kolor 
icol 
ігом 


QPTXTB (С nchar, 


ОРТАТС (inchar, 


Changes the color table by modifying the 
Pallette Registers, EGA card only 


pallette register index (0-15) 
value to be placed in register (0-63) 


char, kolor, icol, irow ) 


Writes plot text. The character variable char 
of length nchar is written at ( column, row ) = 
( icol, irow ) in specified color. 


length of character variable 

character variable corresponding to desired 
text 

specifies color of the text 

starting text column of text 

text row containing text 


char, kolor ) 


Writes plot text at position B. A character 
variable char of length nchar is written at 
position ("B") just beyond the positive X axis 


char, kolor ) 


Writes plot text at position C. A character 
variable char of length nchar is written at 
position ("C") just above and centered on the 
Y axis 


QCMOV (С icol, irow ) 


QSETUP € ndots, 


ndots 


lkolor 
isymbl 


Moves the cursor to a specified column and 
row 


lkolor, isymbl, klrsym ) 


Sets up line type (solid or dashed) and symbols 
to be used by QTABLE when plotting data 


Number of pixels (dots) specifying length of the 


dash. If ndots = 0, then a solid line is drawn. 
The input ndots = 1 is not allowed. By varying 
ndots > = 2, you may create a line composed of 


variable length dashes 

Specifies color of the line 

Specifies symbol to be inserted at data points 
QTABL. If isymbl - -1, then a cross is drawn at 


ООООООООООООЭОЭОООООООООООООООООООООООООООООООООООООООООООООООООООЮО 


klrsym 


QTABL ( itype, 


itype 


npt 
xarray 
yarray 
selectl 


scalexy 
xymax 
fnamel 
fname2 


xtension 


grafprint 
dummy 
Ете 
trtle2 
title3 


title 
xarray 


уаггау 


хатах 
yamax 
gamax 


pi 
fname 


numsets 
numpts 


a,b,c 


їпсапот1 


паса РПІ ІТ :5yvmbl 2 =2, then no symbol 
1s drawn. Otherwise, the character with the 
ascii number "isymbl” is drawn 

Specifies color of the symbol 


NPT (аттау varray ) 


Plots"artabulated function 


Specifies point or line plot. When itype = 0, 
then the plot is a series of symbols at the 
tabulated points. When itype = 1, the symbols 
are joined by straight lines 

The number of points in the function table 

X values in the function table 

Y values in the function table 

Reenaracctes String üsed to hold the input 
values of incidence angle and phi angle 
selection 

A scale factor used in three of the four graphs 
The maximum scaled value out of all the scaled 
X values and all of the scaled Y values 

A character string used to hold the value of 
the input data file name without extension 

A character string used to hold the value of 
the input data file name with extension 

A character string used to hold the value of 
the input data filename extension, i.e. 

ОТЕТ, Τε το τής 

A character string used to allow the user to 
print out a graph from the screen 

Απ input variable used to advance the screen 
The first line of the title for each graph 

The second line of the title for each graph 
The third line of the title for each graph 

The fourth line of the title for each graph 
Data files contain numbers in THET, MAG, PHASE 
format. This polar coordinate format can be 
transformed to rectangular coordinate format. 
XARRAY is the array established to hold all of 
the transformed X values. 


Data files contain mumbers in THET, MAG, PHASE 


format. This polar coordinate format can be 
transformed to rectangular coordinate format. 
YARRAY is the array established to hold all of 
the transformed Y values. 


The maximum value of the X array (XARRAY) 

The maximum value of the Y array (YARRAY) 

The maximum value of the X array (XARRAY) 

and the Y array CYARRAY) 

A constant approximately equal to 22/7 

A character string which equals either F-Phi or 
F-Theta 

A number used for indexing purposes, equal to 
number of incidence angles times the number of 
phi angles 

The number of theta angles per given incidence 
angle and phi angle 

Dummy variables used for incrementation. 


a Index for incidence angles 
b Index for phi angles 
ς Index for theta angles 


A dummy variable used to read and hold the 
value of the incidence angle presently being 
utilized 
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incangl(a) A matrix containing the values of all the 
incidence angles 


1 A dummy integer used as a counter 

phil A dummy variable used to read and hold the 
value of the phi angle presently being ut: ized 

incangld A dummy variable containing an incidence a.gle 
value used for comparison purposes 

thetai A dummy variable used to read and hold the 
value of the theta angle presently being 
utilized 

magi The magnitude of the output for the existing 
incidence angle, phi angle, and theta angle 

phasei A dummy variable used to read and hold the 

value of the phase angle presently being 
utilized 

titleé A dummy character string used to hold an input 
title to be used in an output graph 

title/ A dummy character string used to hold an input 
title to be used in an output graph 

titles A dummy character string used to hold an input 
title to be used in an output graph 

таггау(с) A matrix containing the values of all the theta 


angles in degrees for a given incidence angle 
and phi angle 

rmarray(c) A matrix containing the values of all the 
magnitudes for a given incidence angle 
and phi angle 

parray(c) A matrix containing the values of all the 
phase angles for a given incidence angle 
and phi angle 


calcvar A variable used in intermediate calculatioins 

to determine a graphical scale factor (scalefac) 
scalefac A graphical scale factor 
dummyvar A dummy variable used as a limit in scale factor 


calculations for visual graphic asthetics 


theta(a,b,c) A matrix containing the values of all the theta 
angles for all the incidence angles and all the 
phi angles 

mag(a,b,c) A matrix containing the values of all the 
magnitudes for all the incidence angles and 
all the phi angles 

phase(a,b,c) A matrix containing the values of all the phase 
angles for all the incidence angles and all the 
phi angles 

rad(c) A matrix containing the values of all the theta 
angles in radians for a given incidence angle 
and phi angle 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
VARIABLE DECLARATIONS 

ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


INTEGER A, B, C, I, IREG, {VAR SCALE 

INTEGERX2 ICOL, ТОРТ, IRON, ISYMBL, JCOLI, JCOL2, JROMI 
INTEGERX2 JROW2, KLRSYM, KOLOR, LABEL, LKOLOR, MINOR 
INTEGERXx2 MODE, NCHAR, NDEC, SELECTI 

CHARACTERX2 SELECT, XTENSION 

CHARACTERX3 GRAFPRNT 

CHARACTERX8 DUMMY, FNAME], TITLE7 

CHARACTERX12 FNAME2 

CHARACTERX20 CHAR 

СНАКАСТЕКХ646 TITLEl, TITLE, TITLES, ΙΙ; eS 
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OOOO000000 


ОО 


60 


45 


44 


СЦАКАСТЕК Хъл ГТ ЕВ; ТІТ(Е6С5,19,1) 

REAL ASPECT, CALCVAR, DUMMYVAR, GAMAX, GMAX, INCAHGI 
СЕ ο ο ο. ΙΕ Σε РАТ РІ, МАСІ, HDOTS, NUMPTS 
REALT NUMSETS, RAD, THETAI 

REAL XAMAX, XFIN, XMAJOR, XMAX, XMIN, XORG, XST 

REAL YAMAX, YFIN, YMAJOR, YMAX, YMIN, YORG, YOVERX, YST 
REAL AARRAYC37), BARRAYC37), INCANGL(C5), ΡΗΙ(19) 
puse 19, 57), XPHASECS, 19,3720, THETAC5, 19,572 
COMMON/ONE/ RMARRAY(CS7), SCALEFAC, TARRAY(C3?7) 
COMMONH/ZTHO/  PARRAYC3?7) 

COMMON/THREE/ XARRAYCS7), YARRAYCS7), SCALEXY, XYMAX 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΚΧ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΚΧ 
READ DATA FILES, LOAD ARRAYS, CALCULATE MAXIMUM VALUES 

ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ESTABLISH_ INITIAL VALUES 

pi = 3.1416 

aspect = ‚65 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ESTABLISH INITIAL SCREEN MODE AND COLORS 

mode = 5 

CALL qsmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WELCOME TO MCADGRAF 

write(x*,103) 

read(*,100) dummy 

CALL aqsmode( mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
NOTES CONCERNING GRAPHICAL OUTPUT PART 1 

il te(x,101) 

read(*,100) select 

CALL asmode( mode) 

if(select .ne. 'l' .and. select .ne. '2') go to 60 

mu select „еа. '2') go to | 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
NOTES CONCERNING GRAPHICAL OUTPUT PART 2 

write(X,102) 

read(*,100) dummy 

CALL asmode(mode) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX XXX 
DATA FILE NAME INPUT MENU 

Writecx,91) 

read(*¥,100) select 

pitselect .ne. 'l'.and. select .ne. '2') go to 4 

1 (select δα. 12!) до to 1 

CALL asmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INPUT DATA FILE NAME 

write(xX,92) 

геаа(%,100) fnamel 

CALL asmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
DATA FILE TYPE SELECTION MENU 

write(x,90) 

read(*,100) select 

xtension = select 

CALL aqsmode( mode) 

lf(select .ne. 'l' .and. select .ne. '2' .and. select .ne. (45! 


l.and. select .ne. '64' .and. select .ne. '5' .and. select .ne. '6' 
cand. select .ne. '7') go to 44 


C select еа. '7') бо то ТІ 
Пп ее! ес ва. '6') go to 52 


OOO 


17 
16 


co 


26 


17 (зе! ест „ед, 5!) oo τοσα 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
OPEN DATA ΤΠΕ 
if€xtension θά. !1!) then 

fname2 = fnameb2 tet! 


else if(xtension .eg. '2') then 
fname2 * fnamel//'.tep' 

else if(xtension .eq. '3j3') then 
fname2 = fnamel//!.tmt! 

else 
fname2 = fnamel//'.tmp' 

endi f 


CALL aqsmode(mode) 
ореп(4, file=fname2) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
READ DATA FILE 
read(4,100) title8 
read(4,100) title7 
read(4,86) numsets 
рез и numpts 
а - 
incangld = -1000. 
do I$ i * l;int(numsets) 
c 
read(4,86) incangi 
if(incangi .ne. incangld) b 
if(incangi .ne. incangld) s 
if(i .eq. 1) а = 1 
incangl(a) * incangi 
incangld * incangi 
read(4,86) phii 
phil) = phil 
read(6,100) titleé(a,b,c) 
do 17 c * l;int(numpts) 
геаа(6,88) thetal, magi, phasei 
theta(a,b,c) = thetai 
mag(a,b,c) - magi 
phase(a,b,c) = phasel 
continue 
b = b + 1 
continue 
close(G) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SELECT DATA FILE SUBSET 
CIE. SELECT INCIDENCE ANGLE AND PHI ANGLE) 
write(X,93) a 
бо 25 1-13 

write(x,94) 1, 1псапа 1 (1) 
continue 
мі сес 95) arl на 
wrlte(x,96) 
read(X,X) selectl 
CALL asmode(mode) 
if(selectl .ne. l .and. selectl .ne. 2 .and. selectl .ne. 5 


"|! 
к 


l.and. selectl .ne. 4 .and. selectl пе. 5 .and. selectl .ne. 6 


.and. selectl .ne. 7) go to 5 
if(selectl .eq. atl) go to 44 
17 (зе! ес! еа атг) соса 
а = selecti 
write(x,97) b-1 
do 26 1=1,b-1 

write(X,98) i,phi(1i) 
continue 
wrl1te(x,99) b, b+1 
write(x,96) 


176 


ООО 


OO OO 


OO 


OO 


read(X,X) selectl 

CALL asmode(mode) 

CALL apregt(0,32) 

if(selectl .eq. b) go То 99 

1f(selectl .eq. b+1) go to 1 

b = selectl 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
LOAD ARRAYS 


с = 1 

хатах = 0.0 
уатах = 0.0 
gamax = 0.0 


title2 - title6(a,b,c) 
do 18 c * l,amt(numpts) 
tarray(c) * theta(a,b,c) 


rmarray(c) * mag(a,b,c) 

раггау(с) = рһаѕе(а, ЫЬ, с) 

rad = +аггау( с) X pi / 180 
xarray(c) »* rmarray(c) * sin Crad) 
yarray(c) * rmarray(c) X* cos(rad) 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CALCULATE MAXIMUM VALUES 

if (abs(xarray(c)) .gt. xamax) xamax * abs(xarray(c)) 

if (abs(yarray(c)) .gt. yamax) yamax * abs(yarray(c)) 

continue 

continue 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
БАРБӘӨПГАТЕ SCALE FACTORS 

gamax = abs(xamax) 

if (yamax .gt. хамах) датах * abs(yamax) 

if (abs(gamax - 0) .lt. .000001) датах = .000001 

calevar = alogl0(gamax) 

scalexy = int(calcvar) 

і? (датах .1t. 1.0) scalexy = int(calcvar) - 1l 

scalefac = int(scalexy) 

dummyvar * gamax/(1lO0XXscalexy) 

xymax 7 int(dummyvar + 1) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CLEAR SCREEN 

mode = 16 

CALL qsmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CALL rectmvt(title2,title7,title8,numpts,aspect) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CLEAR SCREEN 

mode = 16 

CALL asmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
Crymlmrectpyvc«titlec, title/,titles,numpts,aspect) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CLEAR SCREEN 

mode = 16 

CALL aqsmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CALL polarmvt(title2,ti1tle7,title8,numpts,aspect) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CLEAR SCREEN 

mode = 16 

CALL qsmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CALL polarrot(title2,title7,title8,numpts,aspect) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CLEAR SCREEN 

mode = 16 

CALL qsmode(mode) 
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XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxxxxxx 
goto 44 
ASPECT RATTO 
52 write(X,104) 
read(x,x) aspect 
mode = 3 
CALL asmode(mode) 
goto 44 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
1 continue 
CLEAR SCREEN 
mode = 5 
CALL qsmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
80 forma t(6X, F 5005 Э ELI οσο F7 " 
86 format(F5.0) 
88 Рогма! (Е5 0 Е SS E7 2) 
90 format(/5X,"EMCAD OUTPUT GRAPHICS PRESENTATION", 75X; 
1 "ПАТА FILE TYPE SELECTION MENOM 7 515 
s S s select the type of data file containing the information of 
2] A Ха 
сааш оез EMCAD automatically created these four data file type 
Sr. SAL OG 
БИ ΤΕ INCIDENCE; ΤΠ TIO, 
"RS TE INCIDENCE; Ғ-ЕНТТТ А 2 
Die TM INCIDENCE, F-THETA'; 773X; 
G. IM INCIDENCE; Επ το. 
57 CHANGE ПАТА ЕГРЕ ОХ, 
6. CHANGE ASPECT RATIO ZZ, 
Ж EXIT THE PROGRAM'!, //5х, 
: пие your Aue es by the choice of 1, 2, ὃ, 6, 5, 6 οι. 
/” 
91 format(775X, "ЕМСАрР OUTPUT GRAPHICS PRESENTATION T 733 
1'DATA FILE NAME INPUT MENU', 775X; 
2'If utilizing a file or series of files created Ьу ЕМСАЮр; 77555 
3'please select option 1. If utilizing any file other than',^/;5X, 
Gtane created by EMCAD, please ensure that the file is in',75x, 
5'proper format, and that the filename extension is one of',/5X, 
6'the four extensions utilized by this program.  These',/5X, 
7'extensions are: .ξοὲ . tep tmt and ешо 525; 
8'If you are unsure as to whether your file is in the proper',/5X, 
9'format with the proper extension, please exit the program',/5X, 
l'and verify before continuing. ',778X, 
2 “INPUT THE FILE OF USER CHDICE' VC P S 
Sue EXIT THE PROGRAM’ ,775x, 
á*Indicate your selection by the choice of 1 or 2.'//5X) 
92 Тогта%(//5Х, ТЕМСАр OUTPUT GRAPHICS PRESENTATION 5x3 
1'DATA FILE NAME INPUT SC 
2'Please input the name of the file containing data that',//5X, 
3'you wish to work with. Please input only the file пате!,//5Х, 
G'i.e. do not include the file extension. The file',/^/5X, 
5'extension will be selected in the DATA ЕЕ ПРЕС ллУх, 
6*SELECTION MENU 775) 
93 format(//5X, 'EMCAD OUTPUT GRAPHICS PRESENTATION',ZZ5X, 
l'SELEGTIONM OF INCIDENCE ANGLE "as 
2'The data file you have selected contains data for ',I2,/5X, 
j'incidence angles. Please select your choice from the',/5X, 
G' following menu. Note, you may only work with one'!,/5X, 
5'incidence angle at a time.',/) 
94 format(5Xx,;, Il, SX, T'INCIDENCE ANGLE = w! Fo O) 
95 Тогма σε SK, 
l'RETURN TO THE DATA FILE TYPE SELECTIONS ENUT S r D 
2'EXIT THE PROGRAM!) 
96 iuc d ao YOUR SELECTION BY THE CHOICE OF A NUMBER TE 
1 --- 0° 4775) 
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ο... δι ου "ЕПСАПр OUTPUT GRAPHICS PRESENTATION',775x, 
SELECTION ОЕ РНІ АМССЕ 97 5х, 
2'The data file you have selected contains сата for *,12,/5Х, 


3'phi angles. Please select your choice from the following',Z5X, 
G'menu. Note, you may only work with one phi angle at a time.',/) 
Ppuwrormateob5x 12,5X, FHL ANGLE = TE F6.2) 


o9crormat(5X,12,5X, 
MERE TURUA OTHE DATA FILE TYPE SELECTION MENU! !,Z5SX, 12; 5X, 
2'EXIT THE PROGRAM!) 

100 format(A) 

101 format(Z5X,'EMCAD OUTPUT GRAPHICS PRESENTATION',77/5X, 
1'NOTES CONCERNING GRAPHICAL HARDCOPY OUTPUT, PART 1',//5Х, 
2'This program allows the user to obtain printer hardcopies!,Z5X, 
3'of the graphs being viewed.  Hardcopy graphic outputs are',/5X, 
G'obtained through the use of the program entitled GRAFPLUS.!,75xX, 
5!'Proper utilization of MCADGRAF requires prior setup of',/5X, 
6'GRAFPLUS. Directions for proper setup are contained',/5X, 
7*in the GRAFPLUS users manual.!,ZZ5X, 
S CONTINUE WITH IHE HEXT SCREEN!, ZZ5X;, 
gU EXIT THE РКОСОКАИ",7/5Х, 
l'Indicate your selection by the choice о? 1 ог 2.'//5Х) 

Ши гогма Е СХЗХ, ЕМСАВО OUTPUT СКАРНТСЗ РЕЕЗЕНТАТТОН" ,//5Х, 
NOTES CONCERNING GRAPHICAL HARDCOPY OUTPUT, —PART 2',7^/5X, 
2'The user makes the decision as to whether a printer ',/5X, 
3'graphical output is desired after viewing the graph. Once',/5X, 
G'a graph is drawn on the screen, the user may request a ',/5X, 
5'printer hardcopy by typing either a "Р" ог "p" and pressing',/5X, 
6'the "ENTER" key. Selection of any other key will simply',^/5X, 
7'advance to the next graph. Spacing on 8 142 by 11 paper !,/5Х, 
8'allows up to 2 graphs per page.',///5X, 
PLEASE PRESS ENTER TO CONTINUE.',5X) 

105 format(^/5X,"'xxxxxxxxxx WELCOME TO MCADGRAF XXXxxxxxxx',/^7/5xX, 
l1'THIS PROGRAM WAS WRITTEN TO BE USED IN CONJUNCTION WITH ',75X, 
2'THE FORTRAN PROGRAM ENTITLED EMCAD.  MCADGRAF ADDS ',75X, 
S'ORAPHING FEATURES TO SUPPORT THE ANALYSIS OF DATA *',75X, 

G' GENERATED AND FORMATTED BY EMCAD. EMCAD WAS WRITTEN BY',/5X, 

ΠΠ. ἡ. MORGAN, WITH MODIFICATIONS BY HIMSELF AND *',75X, 
ΤΕ. M. CONNOLLY. MCADGRAF WAS DEVELOPED AND WRITTEN BY',75x, 
ОТЕ η CONNOLLY TO HELP FULFILL THESIS REQUIREMENTS WHILE!,75x, 
8'UNDER THE GUIDANCE AND THE SUPPORT OF DR. MORGAN. MCADGRAF!,75xX, 
9'WAS WRITTEN AND COMPILED USING RYAN MCFARLAND FURTRAN,!,75xX, 
ШІБКЕКАЕМАТТСЕ, AND бКАЕРІП5 ANY QUESTIONS CONCERNING COLORS,',75<X, 
ο ου ντ CHANCES TO THE PROGRAM, ETIC... SHOULD BE REFERRED ',75X, 
ШЕШІ ГНЕ RESPECTIVE USERS MANUALES; DR. M. A. MORGAN, ОК!,/5Х, 

ша Ее. М. CONNOLLY.',775xXx, 

PLEASE PRESS ENTER 10 СОМТІНЦЕ 5Х) 

104 format(/5X,'XXXXXXXXXX ASPECT RATIO XXXxXxxxxxxx',//5X, 
ITN ORDER TO οὗ DISTORTION ON THE SCREEN, IT IS NECESSARY',^75X, 
2'TO DEFINE THE NUMBER OF ROHS PER INCH AND THE NUMBER OF',/5X, 
S'COLUMNS PER INCH. THESE DEFINITIONS CHANGE AS THE!,75xX, 
CRE SOLUTION OF THE SCREEN, AS HELL AS THE TYPE OF SCREEN',75xX, 
5'CHANGE. THIS PROGRAM WAS WRITTEN IN MODE 16 CEGA) USING',75xX, 
6'A GB-1 VIDEO BOARD AND A NEC MULTISYNC MONITOR. THE!,Z5X, 
Z'OPTIMUM ASPECT RATIO FOR THIS CONFIGURATION MAS .65',77^75X, 
ο λος МРОТ ASPECT RATIO.',775X) 

105 format(F5.3) 

end 


SUBROUTINE RECTMVTCtitle2,title7,title&,numpts,aspect) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxxx 
MAGNITUDE VERSUS THETA 

REC ANGULAR COORDINATE GRAPH 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


O O С) 


O O 


ОО QOO 


ON 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧχΧΧχχΧ 
VARIABLE DECLARATIONS 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
INTEGER A; B; C) I F IREGI VARSA ALERS 

IMTEGERX2 ICOL, IOPI. IRO, I> MBC; СӨРЕЛЕР 

IHTEGERX2 JROL2; SLRSYM, КОГОК, CABEL ΠΠ ΠΠ 

IHTEGERx2 MODE, κα πες ЕЕС 

CHARACTERXZ SELECI pest an 

CHARACTERX3 GRAFPRNT 

CHARACTERXxX8 DUMMY, FNAME], TITLE? 

CHARACTERX12 FNAME2 

CHARACTERX20 CHAR 

CHARACTERX6G TLITEET; TITUE2, TITLES ree eee 

CHARACTERXx6G TITLES, TITEEOC5, T9715 

REAL ASPECT, CALCVAR, DUMMYVAR, GAMAX, GMAX, INCANGI 

REAL INCANGLD, PHASEI, PHII, PI, MAGI, NDOTS, NUMPTS 

REAL NUMSETS, RAD, THETAI 

REAL XAMAX, XFIN, XMAJOR, XMAX, XMIN, XORG, XST 

REAL YAMAX, YFIN, YMAJOR, YMAX, YMIN, YORG, YOVERX, YST 

REAL ААККАУ( 37), BARRAYC37), INCANGL(5), PHIC19) 

REAL MAG(5,19,37), PRASEC5, 19, 3/704) fHetn ιτ) 

COMMON/ONE/ RMARRAYC372, SCALEFAC, TARRAY(C37) 

COMMOH/THO/ PARRAYC37) 

СОММОМ/ТНКЕЕ/ XARRAYC37), YARRAYC3/7), SCALEXY, XYMAX 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SET DISPLAY MODE 

mode = 16 

CALL qsmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SET AXIS COLOR 

CALL apreg(7,6) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
DEFINE GRAPH REGION 


jcoll Ξ 080 
jcol2 = 520 
jrowl = 40 
Jrow2 = 320 
xmin = -0.5 
xmax = 180 
ymin = -0.5 
ymax = xymax 
xorg = 0.0 
vorg = 0.0 
lopt = 
yoverx * ./X(ymax-ymin)/(xmax-xmin) 


CALL qplot(jcoll,jcolZ2,jrowl,jrow2,xmin,xmax,ymin,ymax,xorg, 


lyorg,iopt,yoverx,aspect) 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SET UE LINE Τ.Ε 


ndots = 0 

lkolor = š 
isvmbl = -2 
klrsym = 4 


CALL qsetup(ndots,lkolor,1symbl,klrsym) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


DRAH X AXIS 
xst = 0.0 
xfin = 180 
хтајог = G5 
minor = 0 
label = -1 
ndec + | 


CALL qxaxis(xst,xfin,xmajor,minor, label,ndec) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
DRAN Y AXIS 

yst = 0.0 


180 


20 


yfin = xymax 
ymajor = 2 
ІіТСХИЛах ва I ог. хумах еа. 2 .ог. xymax .eq. 3 .or. xymax 


1.еа. 6 .ог. хутах . еа. 5) утајог = 1 
minor = 1 
label = -1 


аа 5 СУБЕ УИ Утазаг , птіпог; Табе1 , паес ) 
ΕΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


ΙΕ exis ONS SCREEN 


nchar = 20 
char = ‘Theta’ 
kolor = 14 


слава хтобпераг, спаго Ко!ог) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITE TEXT ON SCREEN 


nchar = 20 
char = 'Magnitude' 
kolor = 14 


CALL Gptxtc(nchar, char, kolor) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


KITE TEXT ON SCREEN 


nchar = 64 
kolor = 10 
icol = 15 
irow = 24 


CAL ὉὮαοηδο 11 -Тев , Ко! ог, 1со 1, 1гом) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxXXXX 


WRITE TEXT ON SCREEN 


nchar = 64 

titlel = ‘Scattered Field Complex Amplitudes' 
kolor = 10 

icol = 15 

irow = 23 


CALL qptxt(inchar,titlel,kolor,icol,irow) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


WRITE TEXT ON SCREEN 
4 


nchar = 

kolor = 10 
icol = 15 
irow = 22 


CALL aptxt(nchar,title2,kolor,icol,irow) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


RITE TEXT ON SCREEN 


nchar = 7 
kolor + 10 
Jcol = 15 
irow = 2] 


ccptxttnebar,title/,;kolor;icol,irow) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


HEITE TEXT ON SCREEN 


nchar = 40 

title3 = ! Magnitude Versus Theta' 
kolor = 10 

icol = 30 

"του - 21l 


ιτ аоръехъ(пспаг, 11 1 1е5, КоТог,1со 1, 1гом ) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ἔπ; VALUES 

do 30 i = l,numpts 

aarray(1) <- гтмаггау(т ) /( 10 ЖЖзса |! еху ) 

continue 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ORAPH DATA POINTS 

itype = 1 

npt = numpts 

CALL qtabl(itype,npt,tarrav,aarray) 
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ОООООО 


oi 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
НЕТТЕ ЗСАГЕ FACTOR SCREEII 


nchar = |! 

titleQq = "'S!#8 

Kolor = 14 

icol - 1 

irow = 18 

CALL aptxt(nchar, titles, kolor,i1col,inrow) 
titled = 'c! 

irow = 17 


CALL aptxt(nchar, title4á,kolor;,icol,irow) 
title“ = "а! 

irow = 16 

CALL aptxt(nchar,titleá,kolor,icol;irow) 
title4 = 11! 

irow = 15 

CALL aptxt (nchar, title4,kolor,icol,irow) 
title4 = 'e!' 

ігом 2 14 

CALL aptxt(nchar,title4á,kolor,icol;irow) 
title4 s 'X' 

irow 2 12 

CALL qptxt(nchar, title“4,kolor,icol,irow) 
nchar = 2 

titleá : '10' 

lrow = 

CALL aqptxtCnchar,titleG,kolor,Qicol»)irow) 
nchar 3 

titleG = 'x! 

irow = 9 

CALL aptxt(nchar,title4,kolor,icol;irow) 
titles = ах! 

1гом 3 

CALL «aptxt(nchar;,title4,kolor;icol, irow) 

CALL qemov(2,8) 

write(X,11) scalefac 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
HARD COPY OF GRAPH 7 


nchar = 20 

title3 = ‘Hardcopy ---> Р ог р! 
kolor = 14 

icol = 1 

згом = | 


CALL aqptxt(nchar, title3,kolor,icol,irow) 
CALL GemovC 25,1) 
read(*,100) grafprnt 


nchar = 40 

тез! : 
kolor = 10 

icol- oI 

irow = | 


CALL aptxtí(nchar,;title3,kolor:icol Iron) 

17 (ага ргл? „ед. "Р ог. дгаТргп ea ο ο δι 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
format( 13) 

format CA? 

return 

end 


SUBROUTINE RECTPVTCtitle2,title7,title&,numpts,aspect) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
PHASE ANGLE VERSUS THETA 

RECTANGULAR COORDINATE GRAPH 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


182 


O OOO 


O O 


O0 ОО 


ОО 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxxxxx 
VARIABLE DECLARATIONS 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
M OER A ο ου ΝΕΟ, Val, SCALEFAC 

K ΕΕ το ΟΠ, [τΟΝ, ο. JCOL l, JCOLZ, JROMI 
ШИЕВЕКХФФТЕЦИО ο τος KOLOR, LABEL, LKOLOR, MINOR 

ο Ένκκα MODE, NCHAR, HDEC, SELECT1 

ОНАКАСТЕКХО SELECT, XTEHSIOM 

CHARACTERX3 GRAFPRNT 

CHARACTERXx8 DUMMY, FNAMEl, TITLE? 

CHARACTERX12 FNAME2 

CHARACTERX20 CHAR 

SHOP O CIERX6G TITEEI, TITLEZ, TITLES, TITLE“, TITLES 

EHARACIHERXSSq TITLES, TITLE6C5,19,1) 

REAL ASPECT, CALCVAR, DUMMYVAR, GAMAX, GMAX, INCANGI 
КЕББИТЫСАНЫ: , PHASE], ΡΗΓΙ͂, Pl, MAGI, NDOTS, NUMPTS 

REAL NUMSETS, RAD, THETAI 

REAL XAMAX, XFIN, XMAJOR, XMAX, XMIN, XORG, XST 

REAL YAMAX, YFIN, YMAJOR, YMAX, YMIN, YORG, YOVERX, YST 

REAL AARRAY(37), BARRAY(37), INCANGL(5), PHI(C19) 

БЕДА ОК5 19.57), РНАЗЕ(5,19,35?7), ТНЕТА( 5, 19,57) 

COMMON/ONEZ RMARRAY(37), SCALEFAC, TARRAY( 37) 

COMMONS TWO PARRAY(C 37) 

COMMON THREEZ XARRAY(37), YARRAYC37 05 SCAL EXY ,""XYMAX 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxXXxx 
SET DISPLAY MODE 

mode = 16 

CALL asmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SET AXIS COLOR 

CALL qpreg(7,6) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxx 
DEFINE GRAPH REGION 


jcoll = 080 
gcol2 = 530 
jrowl = 30 

jrow2 = 320 
xmin = -0.5 
xmax = 180 

ymin = -180 
ymax = 180 

xorg = 0.0 

yorg = 0.0 

iopt = ] 


yoverx = 1.5 | | | | | 
CAL L qplot(3coll,3col2,jrowl,jrow2,xmin,xmax,ymin,ymax,xorg, 


lyorg,iopt,yoverx,aspect) 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SEV UF LINE TYPE 


ndots = 0 

lkolor = 3 
isymbl = -2 
klrsym = 4 


CALL asetup(ndots,lkolor,isymbl,klrsym) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


DRAM X AXIS 
xst = 0.0 
xfin = 180 
xmajor = 45 
minor = 
label = 1 
ndec = 1 


CALL axaxis(xst,xfin,xmajor,minor, label,ndec) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


DRAM Y AXIS 


183 


мес - =120 
Τη = (0) 
ymajor = 90 
minor = 1 
label =i 


CALL dqvaxis(vyst,yfi1no)vmasor,minor label OSC 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
WRITE TEXT ON SCREEN 


nchar = 20 
char * 'Theta' 
kolor = 14 


CALL ορι πσ ποσο ο α- kolor) 
ΧΧΝΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ2. «ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITE TEXT ON SCREEN 


nchar * 20 
char = 'Phase' 
kolor = 14 


CALL qptxtc(nchar, Char, Keon. 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITE TEXT ON SCREEN 


nchar = 64 
kolor = 10 
icol = 15 
1гом = 26 


CALL ар х+(песпаг, #1 Е Тед, Ко! ог, 1со! зге и 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


WRITE TEXT ON SCREEN 


nchar Ξ 66 

titlel с 'Scattered Field Complex Amplitudes' 
kolor = 10 

icol = 15 

1гом + 25 


CALL aptxtCnchar,titlel koloro1colo ow. 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


HRITE TEXT ON SCREEN 


nchar * 66 
kolor = 10 
τεσ] 
irow + 22 


CALL qptxttneobar,trtle2,kolor 69153057900) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITE TEXT ON SCREEN 


nchar = 7 
kolor * 10 


icol =s 

“Пом = αἱ 

CALL aptxti(nchar;title/7.kolor,:icol iron) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITE TEXT ON SCREEN 


nchar = 40 

title3 = ' Phase Versus Theta! 
kolor = 10 

1661 - 30 

гош = 2l 


CALL qptxt(nchar,titleš,kolor,l1col,lrow) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ΟΝΑΡΗ DATA POINTS 

itype = 1 

npt = numpts 

CALL atablCitype,npt,tarray,parray) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
HARD COPY OF GRAPH ? 


nchar - 20 

title3 = 'Hardcopy ---> P or p! 
kolor = 14 

1с01 = 40 
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QO0OO0 00400000 


ОО 


ОО ОО 


ОО 


10 


iro = l 

ЕЕ арс еСпспаг titlesS,kolor; 1col,irow) 
CALL астом (65,1) 

read(*,100) grafprnt 


nchar = 40 

tities = * я 
kolor = 10 

icol = 840 

1гом = l 


RGP í nchar:;:title3,kolór;icol;irow) 

OTAGO .ed. 'P' „ог. grafprnt .eq. 'p') CALL prtsc 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
format(A) 

return 

end 


SUBROUTINE POLARMVT(title2,title7,title8,numpts,aspect) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
MAGNITUDE VERSUS ΤΗΕΤΑ 

POLAR COORDINATE GRAPH 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
VARIABLE DECLARATIONS 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXkxXXXXxxxxxxxxxxxxx 
ШШПЕЗЕК А) 8 С, І, ІБЕС, ТІУДІ, 5САТЕРАС 

ШИПНЕБЕКЖЕТТЕПІ, ІОРТ, IROW, ISYMBL, JCOLI, /С012, JRON1 

INTEGERx2 JROW2, KLRSYM, KOLOR, LABEL, LKOLOR, MINOR 

INTEGERx2 MODE, NCHAR, МОЕС, 5Е! ЕСТ1 

CHARACTERX2 SELECT, XTEHSION 

CHARACTERX3 GRAFPRNT 

CHARACTERXx8 DUMMY, FNAME], TITLE7 

CHARACTERX12 FNAME2 

CHARACTERX20 CHAR 

ντ ετεκ τε πι ες TITLES, TITLE, TITLES 

СНАКАСТЕКж64 TITLES, ΤΙΤΙΕό(5 19,1) 

REAL ASPECT, CALCVAR, DUMMYVAR, GAMAX, ΟΜΑΧ, INCANGI 

СЕ INCANGLD, PHASELI, FPHIT, PL, MAGI, NDOIS, NUMPTS 

REAL NUMSETS, RAD, ΤΗΕΤΑΙ 

REAL XAMAX, XFIN, XMAJOR, XMAX, XMIN, XORG, XST 

REAL YAMAX, YFIN, YMAJOR, YMAX, YMIN, YORG, YOVERX, YST 

КЕП ДАККАҮСЗ5/), ВАКБКАҮС57), ІМСАМбІС(5), РНІ(19) 

ШЕДИИМАВС5, 199357)/ "РПАЗЕС5,19,57), ТПЕТА(5,19,37) 

COMMONZONEZ RMARRAY(37), SCALEFAC, TARRAY(37) 

COMMON/TMWO^/ PARRAYC37) 

COMMON/THREE/ XARRAYC372, YARRAYC37), SCALEXY, XYMAX 
XX3«XXXXXXXXXXXxxxxxxxxxxxxXxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SET DISPLAY MODE 

mode = 16 

CALL qsmode(mode) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SEIT AXIS COLOR 

CALL apreg(C7,6) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxXxxxxxxx 
s ας ο αρ πι VALUES 

do 10 1 = l,numpts 

aarray(1) = xarray(1)Z(10xxscalexy) 

barray(1) уаггау(1)/(10Х%%е<са1еху) 

continue 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
DEFINE GRAPH REGION 


1са 1 = 080 
2с012 = 400 
jrowl = 10 
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OM 


jrow2 = 280 


xmin < -хутах/2 

xmax 2 хутах ` 
ymin = ~xymax 

ymax = xymax 

xorg 5 7 0 

yorg = 0.0 

liopt - ] 


yoverx = По ва | 
CALL aplot(3coll,jcol2,jrowl,3row2,xmin,xmax,ymin,ymax,xorg, 


lyorg,iopt,yoverx,aspect) 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SET ὑΡ ΕΙΠΕ ΤΕ 


ndots = 0 

lkolor * 3 
isymbl = =2 
klrsym = 4 


CALL qsetup(ndots,lkolor,isymbl,klrsym) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


ΡΚΑΝ Χ ΑΧΙ5 

xst * 0.0 

xfin = xymax 

xmajor = 2 

if€xymax .le. 5) xmajor = 1 
minor = 1 

label = 1 

ndec = 1 


CALL axaxis(xst,xfin,xmajor,minor, label,ndec) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


DRAM Y AXIS 


yst * -xymax 
yfin = xymax 
ymajor = xmajor 
label = 0 


CALL ayaxis(yst,yfin,ymajor,minor,label,ndec) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


WRITE TEXT ON SCREEN 


nchar = 10 
char = 'Theta=90' 
kolor = 14 


CALL aptxtbl nchar, char, kolor) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


МЕТТЕ Per  ϱ SCREEN 


nchar = 10 
char = 'Theta=0° 
kolor = 14 


icol = 10 

irow = 

CALL aqptxt( nechar, char, kolor, icol; irow) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITE TEXT ON SCREEN 


nchar = 10 

char = 'Theta=180' 
kolor = 14 

icol = 10 

ігом = ] 


CALL qptxt( nchar, char, kolor, του τυ... 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


WRITE TEXT ON SCREEN 
64 


nchar 3 

kolor = 10 
icol = 15 
1γοι’ > 24 


CALL aptxt( nchar, *itle8, Ко! ог, icol “σι 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


WRITE TEXT ON SCREEN 
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nchar = 64 
titlel = 'Scattered Field Complex Amplitudes' 


kolor = 10 
Jeol = 15 
1που = 23 


ο ορ τί πσραος, ctitlel, kolor, icol, irow) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


ΤΕ ΤΕ ON SCREEN 


nchar = 64 
kolor = 10 
icol = 15 
irow = 22 


(C ΕΕ δν (πσῃβασς, t1tle2,kolor,1col,lrouw) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


ЧЕТТЕ τει ОН SCREEN 


nchar = 11 

title3Ë5 Polar Plot, * 
kolor = 10 

icol = 15 

irow * 2l 


CALL aqptxt(nchar, title3,kolor,icol,irow) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


WRITE y ON SCREEN 


nchar = 
kolor = T 
icol - 27 
ігом = 2] 


CALL aptxt(inchar, title7,kolor,icol,irow) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


WRITE TEXT ON SCREEN 


nchar = 40 

titles = ' Magnitude Versus Theta’ 
kolor = 10 

теа = 35 

irow = 2] 


CALL aqptxtCnchar, title3,kolor,icol,irow) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


WRITE SCALE FACTOR ON SCREEN 


nchar = 1 
title4 = (5! 
kolor = 14 
icol = 10 
irow = 17 


CALL qptxt(nchar,title4,kolor,icol,irow) 
ісе = "с! 

irow + 16 

СОЕ арахебпераг., 1 31 ед, Ко! ог, 1со! ,1гом ) 
title4 * "а! 

irow = 15 

CALL aptxt(nchar,titled,kolor,icol,irow) 
ша ед = 71! 

irow + 14 

CALL aptxt(nchar,titleá,kolor,icol,irow) 
кеа = 'е' 

irow = 13 

CALL aptxt(nchar,titleá,kolor,icol,irow) 
titled = 'X! 

lrow = ll 

CALL aqptxt(nchar,title4,kolor,icol,irow) 
псһаг = 2 

сеа = "I0! 

irow = 9 

CALL aptxt(nchar,title4á,kolor,icol,irow) 
nchar = 1 

titleG = "Хх! 

ігом = 
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OO O 


O, O) O O) O) O O, O) O) O O 


CALL qptxt(nchar,; їс еб; котос ТСО Iron) 

title, : (Хх! 

irow = 7 

CALL qptxt(nchar,tí1itle%,kolory”icol,l1row) 

CALL “qemnov Cs, 7) 

write(*,12) scalefac 

See discussion of QCMOV for problem requiring this fix 
ἘΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
GRAPH DATA POINTS 

itype = 1 

npt = numpts 

CALL qtablCitype,npt,aarray,barray) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
HARD COPY OF GRAPH ? 


nchar = 20 

titles = "Нагасору ο. 
kolor = 14 

icol- = 50 

irow = | 


CALL aptxt(nchar, title3,kolor,icol,irow) 
CALL acmov (75,1) 
read(x,100) grafprnt 


nchar = 26 

titles =! : 
kolor = 10 

icol = 50 

irow = ] 


CALL qptxt(nchar, title3š;,kolor,icol, irow) 

1f(grafprnt .eq. 'P' .or. ογδίεργος - μι r t= 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
format(9x,13) 

format(A) 

return 

end 


SUBROUTINE POLARROT(Ctitle2,title7,title8,numpts;,aspect) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
MAGNITUDE VERSUS THETA 

POLAR COORDINATE GRAPH 

AXIS ROTATED MINUS NINETY DEGREES 

X X X XX XXX X X Xx Xxx xxx JE JE JE JE JE JE JE JE JE JE JE XE X x 1 XX x JE JE JE JE JE JE JE JE JE E JE JE J€ JE ЭЭЭЭ 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
VARIABLE DECLARATIONS 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INTEGER A, B, С, 1, ТКЕС, IVAL, SCALEFAC 

INTEGERX2 ICOL, IOPT, IROM, ISYMBL, JCOLl, JCOL2, JROMI 
INTEGERX2 JROW2, KLRSYM, KOLOR, LABEL, LKOLOR, MINOR 
INTEGERx2 MODE, NCHAR, NDEC, SELECTI 

CHARACTERX2 SELECT, XTENSION 

CHARACTERX3 GRAFPRNT 

CHARACTERx8 DUMMY, FNAME], TITLE7 

CHARACTERX12 FNAME?2 

CHARACTERX20 CHAR 

CHARACTERX64 TITLEI, ТІШРЕ2, ГІНЕ, TOMEG ИВЕ 
CHARACTERX64 TITLE8, TITLE6(5,19,1) 

REAL ASPECT, CALCVAR, DUMMYVAR, GAMAX, GMAX, INCANGI 
REAL INCANGLD, PHASEI, PHII, PI, MAGI, NDOTS, NUMPTS 
REAL NUMSETS, RAD, THETAI 

REAL XAMAX, XFIN, XMAJOR, XMAX, XMIN, XORG, XST 

REAL YAMAX, YFIN, YMAJOR, YMAX, YMIN, YORG, YOVERX, YST 
КЕА| ААККАУ( 57), ВАККАУ( 57), INCANGL(52, PHI(19) 

REAL МАО(5,19,57), РНАЗЕС5,19 22790 ТИЕП ОЗ 
COMMON/ONEY RMARRAY( 37), SCALEFAC, TARRAY(C37) 
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О О 


ОО ОО 


О О 


ΟΟΜΜΟΝΜ/ΤΗΟ/ PARRAYC37) 

COMO THREEZ XARRAYC37), YARRAY(37), SCALEXY, ΧΥΜΑΧ 
ΧΧΧΧΧΑΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SEIT DISPLAY MODE 

mode = 16 

CALL asmode(mode) 
ΧΚΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΚΧΧΧΧΧΧΧΧΧΧΧΧ 
SET AXIS COLOR 

CALL apreg(7,6) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
DEFINE GRAPH REGION 


jcoll = 080 
3со12 = 550 
Jrowl = 10 
Jrow2 = 520 
xmln = -хутах 
хтах = хутах 
ymin 2 -хутах/4 
ymax = xymax 
xorg = 0. 
yorg = 0.0 
iopt = 1 


yoverx = 0.5 | | | f f 
CALL qplot(jcoll,jcol2,3Jjrowl,Jrow2,xmln,)xmax,vymln,ymax,>xorg, 


lvorg,iopt,yoverx,aspect) 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SEV UF LINE TYPE 


ndots = 0 
lkolor 5 2 
isymbl = -2 
klrsym = 4 


CALL asetup(ndots,lkolor,isymbl,klrsym) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


DRAN X AXIS 

xst = -xymax 

xfin = xymax 

xmajor = 2 

if(xymax .le. 5) xmajor = 1 
minor = ] 

label = 0 

паес = 1 


CALL axaxis(xst,xfin,xmajor,minor,label,ndec) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


DRAM Y AXIS 


yst = 0.0 

yfin = xymax 
ymajor = xmajor 
label = 1 


CALL ayaxis(yst,yfin,ymajor,minor,label,ndec) 
XXXXXXXXXxXxxxXXXXXXXXXXXXXXXXXXXxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE TEXT ON SCREEN 


nchar = 10 

char = !"ТНеХа-180" 
kolor = 14 

icol = 65 

ігом : 2 


CALL aptxt€nchar,char,kolor,icol,irow) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
WRITE TEXT ON SCREEN 


nchar = 10 
char = 'Theta=90'! 
kolor = 14 


CALL aptxtc( nchar, char, kolor) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITE TEXT ON SCREEN 

nchar = 10 
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20 


char = 'Theta =: и" 


kolor = 14 
icol = 6 
irow = 2 


CALL aptxt(nehar, Char,.kolor,2¢col, row 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


WRITE TEXT ON SCREEN 


nchar < 69 
kolor + 10 
1с01 15 


irow * 26 

CALL aqptxt€ mehar, titles, kolor: σσ. r ou 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITE TEXT ON SCREEN 


nchar = 66 
titlel 7» 'Scattered Field Complex Amplitudes' 
kolor = 10 


σοι 5115 

1гом + 25 

CALL qptxt( nchar, titlel, Ко! ог, ἀσοὶ το. 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
NRITE TEXT ON SCREEN 


nchar * 66 
kolor = 10 
icol = 15 


irow = 22 

CALL aptxtCnchar, title2,kolor,icol,irow) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITE TEXT ON SCREEN 


nchar = ІІ 

titles = ροιςς Riot," 
kolor = 10 

icol = 15 

το - εἰ 


CALL aptxt(nchar,title3,kolor,icol,irow) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITE TEXT ON SCREEN 

nchar = 7 


kolor = 10 
icol =: 27 
irow = 2] 


CALL qptxt(nchar,title7,kolor,icol,irow) 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


WRITE TEXT QN SCREEN 


nchar = 60 

title3 = ' Magnitude Versus Theta' 
kolor * 10 

icol = 35 

irow - 2l 


CALL ар хЕ( пспаг, 1 ЕТ е5, Ко! ог, пса Iron) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


WRITE TEXT ON SCREEN 


nchar < 64 

title5 = 'Axis Rotated Minus Ninety Degrees! 
kolor = 10 

1600" 315 

1гом = 20 


CALL qptxt(nchar,title5,kolor,l1col,l1row) 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
SCALE TAND YF VALUES 

do 20 i = l,numpts 


aarray(i) 2 -уаггау(1)/(10ХХе<са!1еху) 
Баггау(1) * xarray(i)/(10xxscalexy) 
continue 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
GRAPH DATA POINTS 
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ООО 


itype = 1 

npt = numpts 

CALL qtabl(itvpe,npt,aarrav,barray) 

XXX X XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxXXXXXXXXXXxxxXxx 
RITE SCALE FACTOR OH SCREEN 


nchar = 1 
title = 15! 
kolór = 14 
icol 5 


ігом 16 

CALL aqptxt(inchar,titleG,kolor,icol,irow) 
IIIS, = 5c! 

irow = 15 

CALL aptxt(nchar,titleó,kolor,icol,;irow) 
titleG = Τα! 

ігоы 5 16 

CALL aqptxt(nchar, title¢,kolor,icol,irow) 
ticle@ = "1! 

irow = 15 

CALL aptxt(nchar,titleó4,kolor,icol,;irow) 
«#1 #1е4 = “е! 

ігом - 12 

CALL aqptxt(nchar,title4,kolor,icol,irow) 
titles = 'X' 

irow = 10 

CALL aptxt(nchar, title4,kolor,icol;,irow) 
nchar 2 

title = 110! 

1гом = 8 

ο αρ πίποπανγ, title&6,kolor,;icol,irow) 
nchar = 1 

title = 'x' 

irow = 

CALL qptxt(nchar,title4,kolor,icol,irow) 
titleG = 'x! 

irow = 6 

CALL aptxt(nchar, title4,kolor,icol,irow) 
Call QCMOV(4,6) 

write(X,13) scalefac 

See discussion of QCMOV for problem requiring this fix 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
HARD COPY OF GRAPH ? 


nchar * 20 

trtle5 = ‘Hardcopy ---> P or p' 
Kolor = 14 

icol = 30 

irow = 1] 


CALL aptxt(nchar,title3,kolor,icol,irow) 
CALL астоу(55, |) 
read(*,100) grafprnt 


nchar = 40 

ticles = ' : 
Kolor = 10 

icol = 30 

irow = 1 


CALL aptxt(nchar,title$,kolor,icol,irow) 

if orafprnt θα. "Р" ог. агарга .eq. δ!) CALL prtsc 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
format( 4x, I 3) 

format( A) 

return 

end 
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OO00000 


11 


22 
100 


APPENDIX J 
SPHERE SOURCES ODE 


PROGRAM SPHERE 


GENERATING THE RADII & Er, Ur AT ROWS I & I*1 

FOR A LAYERED DIELECTRIC SPHERE HAVING 

Z-AXIS OFFSET/LAMBDA = D. Written By M. Morgan 2414788 
Modified By E. Connolly 473/788 


REAL RADC5),ROC5) 

COMPLEX ER,UR 

CHARACTERX64 FNAME,HDRI,HDR2 

PI=3.1415927 

WRITECX, 110) 

READ(X, 100) FNAME 

OPEN(1,FILE=FNAME,STATUS='UNKNONN'!,ACCESS='SEQUENTIAL!, 
FORM='FORMATTED') 

WRITECX, 111) 

READCX,100) FNAME 

OPENCZ, FILE=FNAME, STATUS= "UNKNOWN! ,ACCESS="SEQUENTIAL', 
FORM=*FORMATTED®) 

WRITEC(X, 112) 


READ(X, 100) HDRI 
ARITECX; 1153) 

READ(X; 100) HDR2 
FIRIIEC1,105)”HDRI 
WRITEC1,105) HDR2 
МЕІТЕСЖ, 116) 

READ OX) NL 
NRITEC1,101) NL 

DO 11 L*1,NL 

WRITEC(X,X) 'Layer No:',L 
WRITECX,x) 

WRITECX,X) 'Enter Radius/Lambda: ! 


READ(x,%x) RO(L) 

WRITEC(X,X) 'Enter Real Er of This Layer: у 
ΚΕΑΡ(Χ,Χ) ΕΚΚ 

WRITECX,X) 'Enter Imag Er of This Layer: ! 
READCX,X) ERI 

ERZCMPLXCERR, ERI) 

URSCMPLXC1.0,-1.E-6) 

WNRIIEC1,102) ER,UR 

CONTINUE 

WRITECX,X) 'Enter Z-Axis Offset/Lambda: T 
КЕАрСХ ХУ) D 

HWRITECX,X) 'Enter IBIG: ' 

READ(x,x) IBIG 

WRITEC2,104) IBIG,NL+1 

RMIN-:2.XPIXCROC1)D -D) 

RMAX-z2.XPIXCROCNLO-*D) 

МЕТТЕС2,102) ЕМІН,ЕМАХ 

DT=PI/CIBIG-1.0) 

DO 33 I=1,IBIG 

T=(I-1.0)xDT 

DS2=(DxSIN(T))%xx2 

СЕСЕ СТ 

DO 22 L=1,NL 

RADCL222.0XPIXCDXC* SQRTCROCLOXXx2-DS2)) 
WRIPEC2, 103) I; CRADCLO;U S1 NI 

FORMAT (A) 
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FORMATCI 

FORMAT (4 ) 

ҒОКМАТ(1 4.6)) 
РОКМАТС 2 
paga C я, 
КОЕ ТОО TNO FILES ARE REQUIRED TO HOLD PROGRAM OUPUT',777X, 

TT TUE E ENTE THE НАМЕ ОЕ THE FILE TO BE USED TO ',777X, 

ЕЙИШЕШ ІНЕ PROGRAM MATERIAL PARAMETER DATA. PLEASE ',777X, 

EN ο ο THE EXTENSION. -I. E. ГІГЕМАМЕ.ПАТ!,//7Х) 

ο Το EASE ΕΝΜΤΕΕ ΤΗΕ NAME OF THE FILE TO BE USED ',777X, 

' πο ο ο ΤΠῈ PROGRAM OUTPUT DATA. PLEASE INCLUDE THE',777X, 

ΤΕΕ ΤΟΝ. Т. C ESNEEETLENAME:DAT',777X) 

ОКИЛ ТОО THE HEADERS ALLON THE USER TO IDENTIFY THIS',ZZ7X,; 

ISET OF DATA FROM ALL. OTHER SETS. '; 

TT REESE ENTER HEADER #1 (64 CHARACTERS MAX)',777X) 
FORMAT(//7X, 'PLEASE ENTER HEADER #2 (6ά CHARACTERS MAX)',.,22Z7X) 
ο οτι ος Ετος THE NUMBER OF LAYERS (.1Е. 5)",//7Х) 

CONTINUE 
EBOSECI) 
CLOSEC2) 
STOP 

END 


— (n — n 
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00000000000 


OOOO 


ООО О 


APPENDIX K 
CONES SOURCE CODE 


PROGRAM CONE 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


GENERATING THE RADII & Er, Ur AT RONS I & I+1 
FOR А LAYERED DIELECTRIC CONE 
Written By E. Connolly 47/19/88 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
VARIABLE DECLARATIONS 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
REAL RÁDC1000,5)5, PI, DELTHETA, DEGTORAD 

REAL THETA, H(525, X(5), Y(5) 

REAL ALPHA1(5), ALPHAC5), B, M, RL,GAMMA 

INTEGER 1510, 

COMPLEX ER,UR 

СНАКАСТЕКх64 ЕМАМЕ, НОК1, НОК2 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INITIAL VALUES 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
КМІМ = 10000 

RMAX = -10000 

Fl-5. 1815927 

DEGTORAD = PI/180. 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 

INPUT NECESSARY DATA 

ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 

ИКТТЕ( ж, 110) 

READCX,100) FNAME 

OPENC], FILE=FNAME, STATUS=*UNKNOWN',ACCESS="SEQUENTIAL', 
FORM='FORMATTED' ) 

VRITECX,111) 

READ(X,100) FNAME 

OPENC2,FILE=FNAME, STATUS="UNKNOWN!,ACCESS="SEQUENTIAL', 
FORM=' FORMATTED!) 

WRITECX, 112) 

READCX,100) HDR1 

HRITE(X,113) 

READ(CX,100) HDR2 

WRITEC1,105) ΗΡΕΙ 

WRITEC1,105) HDR2 

WRITEC 


WRITECX, 
WRITECX, X 
WRITECX,10) 

READ(X,X) ALPHAICL) 

ALPHACL) = CALPHAICL)I72.) * DEGTORAD 
WRITECX,*x) 

WRI:ECX,15) 

READCX,X) HCL) 


Layer No:', 1 


чм» 


VMRITE(x,x) " Enter Real Er of This Layer: '' 
ΕΕΑΡ(Χ,Χ) ΕΕΚ 
WNRITE(Cx,x) ' Enter Imag Er of This Layer: '! 


ОООО 


OOOOOON O N0 


ОООО 


11 


КЕАШХ,Ж) ЕКТ 
ER-CMPLXCERR, ERI) 
URSOMPPEXCI-U 1. E-5) 
ОИЕ СИ 102) ЕК, ШК 
CONTINUE 

ШМЕГТЕСГЖ, Ж) ‘Enter IBIG: ! 
READ(x,x) IBIG 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
CALCULATIONS 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΚΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
DELTHETA * PIZ(IBIG-1) 

DO 22 ІЗІ,ІВІб 

THETA = (І-1) ж ОЕ! ТНЕТА 

DO 33 L=1,NL 
YCL) = HCEL)Z2 

XCL) = CYCLI)¥TANCALPHACL)) 
CHECK FOR HORIZONTAL LINE 
GAMMA = ATANCXCLIZYCL)) 

IFCCPI-GAMMA-THETA) .LT. .0001) THEN 
YFOUND = -YCL) 
XFOUND = YFOUND * TANC THETA ) 

e = 2XPIXSQRTC(XFOUNDXX2 + YFOUNDXX2) 

S 

SLOPE OF LINE CONTAINING CORNERS OF CONE 
M = ( -HCL)/XCL) ) 
ΕΕΤΤ OF LINE CONTAINING CORNERS OF CONE 

= YCL 

RECALL EQUATION OF ANY LINE IS Y = MX + B 
THE LINE CONTAINING THE I-TH AND K-TH DATA POINTS INTERSECTS 
THE LINE DRAWN FROM THE ORIGIN THROUGH THE DESIRED DATA POINT 
AT THE DESIRED DATA POINT XFOUND AND YFOUND. THE EQUATION FOR 

THE LINE CONTAINING THE I-TH AND K-TH DATA POINTS MUST THEN 
ALSO WORK AT THE DESIRED DATA POINT. WRITING XFOUND AND 
YFOUND IN TERMS OF THEIR ANGLE THETACN), AND SUBSTITUTING 
YFOUND = B Z ( 1 - ( M x TANC THETA ) ) ) 
XFOUND = YFOUND * TANC THETA ) 

ee = 2XPIXSQRTCOXFOUNDXX2 * YFOUNDXx2) 

N 

IFCRADCI,L) .GT. RMAX) RMAX = RADCI,L) 
IFCRADCI,L) „ЕТ. RMIN) RMIN = RADCI,L) 


33 CONTINUE 
22 CONTINUE 


44 


10 


15 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
OUTPUT 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
ИКТТЕ( 2,104) IBIG,NL+1 

WRITEC2,102) RMIN, RMAX 

DO 44 I = 1,181С 

ШЕӘШПЕС2, 105) I,CRADCI,L),L=1,NL) 

CONTINUE 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
FORMAT(/8X,'Please enter the apex angle, ALPHAl, in degrees, of ' 


1/78X,'the offset cone. Hote, apex angle is the angle from one side' 
2/8X,'of the cone to the other side of the cone. ) 


FORMATC/8X, 'Please enter the height, H, of the '/ 
8X,' cone. 2) 

FORMATCA) 
FORMATCI5 
FORMAT (4¢ 
FORMATCI5, 
FORMATC2CI1 
FORMAT." и 
FORMAT(C/77 
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1'DATA. PLEASE ENTER THE NAME GF THE FILE TO Ρε ει TIME 
2'HOLD THE PROGRAM MATERIAL PARAMETER DATA. PLEASE ',777X, 
S'INCLUDE THE EXTENSION. Т. Е FIL ENANE ΙΡ 
lll FORMAT(//7X,'PLEASE ENTER THE МАМЕ ОЕ THE ΕΙ ΕΤ ΡΕ ΟΕ) P АЫ 
l'TO HOLD THE PROGRAM OUTPUT UATE PLEASE TE E ΤΕΥ, 
2'EXTENSION. I. E. ΕΙΕΕΙΕΙ  .. 
112 FORMAT(ZZ7X, "THE HEADERS ALLON THE 15; τ ΤΕΙ; ОИУ 
1'SET GF DATA FROM SLL OTTER SET 
2//7X,"*PLEASE ENTER HEADER #1 (64 СНАКАСТЕК5 МАХ)",//7Х) 
113 FORMAT(//7X,'PLEASE ENTER HEADER #2 (64 CHARACTERS MAX)',777X) 
116 FORMAT(ZZ7X; 'ENTER THE NUMBER OF ΙΑΥΕ ο μΕ ο Т 
CLOSE D) 
ο εςς) 
STOF 
END 
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APPENDIX L 
CYLINDER SOURCE CODE 


PROGRAM CYLINDER 
ΚΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 


GENERATING THE RADII & Er, Ur AT RONS I & I+1 
FOR A LAYERED DIELECTRIC CYLINDER 
Written By E. Connolly 4/19/88 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
VARIABLE DECLARATIONS 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
REAL RADC1000,5), PI, DELTHETA, DEGTORAD 

REAL THETA, H(5), R(5), Y(5) 

REAL RL, GAMMA 

INTEGER IBIG, I 

COMPLEX ER,UR 

CHARACTERX64 FNAME,HDR1,HDR2 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
INITIAL VALUES 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
RMIN - 10000 

RMAX = -10000 

PI=3.1415927 

DEGTORAD - PI/180. 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 

INPUT NECESSARY DATA 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxxxx 

WRITE(X,110) 

КЕАО(ж,100) FNAME 

OPEN(1,FILE=FNAMWME,STATUS='"UNKNOHWN!,ACCESS='*SEQUENTIAL!', 
FORM='FORMATTED') 

WRITECX,111) 

READC¥,100) FNAME 

OPENC2,FILE=FNAME,STATUS="UNKNOWN*®, ACCESS="SEQUENTIAL', 
FORM=* FORMATTED!) 

WRITECX,112) 

ΚΕΑΡ(Χ,100) HDR1 

WRITE(X,113) 

READC*,100) HDR2 

WRITEC1,105) HDR1 

WRITEC1,105) HDR2 

ИКТТЕ(х, 114) 

READ(x,x) NL 

WRITEC1,101) NL 

DO 11 L721,NHl 

WRITECX,X) 'Layer No: *',L 

WRITE(CX,X) 

HRITECX,15) 

READCX,x) HCL) 

ИЕІТЕСЖ, 16) 

READCX,X) R(LD 

WRITE(X,X) 'Enter Real Er of This Layer: : 

READC(X,X) ERR 

WRITECX,X) ‘Enter Imag Er of This Layer: ! 

КЕАр(х,х) ERI 

ER-CMPLXCERR, ERI) 
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URSCMPLX(1.0,-1.E-6) 
ИКТТЕС1, 102)" ЕК: ШЕ 

11 CONTIHUE 
WRITECX,X) 'Enter IBIG: ! 
READ(x,%x) IBIG 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧχΧΧΧΧΧ 
CALCULATIONS 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
DELTHETA = PIZ(IBIG-1) 

DO 22 I=1,IBIG 
THETA - (1-1) * DELTHETA 
DO 35 L*1,NL 
YCL) = HCL)/2 
CHECK FOR HORIZONTAL LINE 
GAMMA = ATAN(R(L)ZY(L)) 
IF((PI-GAMMA-THETA) .LT. .0001) THEN 
YFOUND = -Y(L) 
XFOUND = YFOUND x TANC THETA ) 
RADCI, L) = 2XPIXSQRT(XFOUNDX*X2 * YFOUNDXx2) 
CHECK FOR HORIZONTAL LINE 
ELSE IF((GAMMA-THETA) „СТ. 0) THEN 
YFOUND = YCL) 
XFOUND = YFOUND * TANC THETA ) 
бл - 2XPIXSQRTCOXFOUNDXx2 * YFOUNDXxx2) 
LS 
CHECK FOR VERTICAL LINE 
XFOUND = RCL) 
YFOUND = XFOUND 7 TANC THETA ) 
RADCI,L) = 2ZXPIXSQRTCXFQUNDxx2 * YFOUNDXx2) 
ENDIF 
IF(CRADCI,LO .GT. RMAX) RMAX * RADCI,L) 
IFCRADC(I,L) .LT. RMIN) RMIN = RADCI,L) 
33 CONTINUE 
22 CONTINUE 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
OUTPUT 
ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΚΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
WRITEC2,104) IBIG,NL*1 
WRITEC2,102) RMIN, RMAX 
DO 44 I = 1,181С 
WRITEC2,103) I,CRADCI,L),L=1,NL) 

44 CONTINUE 


ΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧΧ 
15 FORMAT(/8X, 'Please enter the height, H, of the *^/ 

18X,' cylinder (from base to base). 25 
16 FORMAT (/8X, "Please enter the radius R, of the '/ 

18X,' cylinder. e 
100 FORMAT (A) 


101 ҒОКМАТ(15) 

102 ҒОКМАТ(9(Е16.6)) 
105 ҒОКМАТ(15,5(Е14.6)) 
106 ҒОКМАТ(2(15)) 


105 ҒОКМАТ(" ',A) 

110 FORMAT(//7X,'TMO FILES ARE REQUIRED TO HOLD PROGRAM OUPUT ',777X, 
1'DATA. PLEASE ENTER THE NAME OF THE FILE TO BE USED TO ',ZZ7X, 
e'HOLD THE PROGRAM MATERIAL PARAMETER DATA. PLEASE ',//7X, 
S'INCLUDE THE EXTENSIDBN. I. E. FILENAME DAT 777) 

111 FORMAT(ZZ7X, "PLEASE ENTER THE NAME OF THE FILE TO BE USED '!',ZZ7X, 
1'TO HOLD THE PROGRAM OUTPUT DATA. PLEASE INCLUDE ТНЕ!,//7Х, 
¿'EXTEHSION. I. Е. FILENAME. DAT ZIL) 

112 FORMAT(C//7X,'THE HEADERS ALLOW THE USER TO IDENTIFY THIS', 7/7/73 
1'SET OF DATA FROM ALL OTHER SETS. 
2//7X,*PLEASE ENTER HEADER #1 (64 CHARACTERS МАХ)",//7Х) 


198 


113 FORMAT(//7X,'PLEASE ENTER HEADER 32 (64 CHARACTERS MAX) !,777X) 
114 FORMAT(//7X,'ENTER THE NUMBER OF LAYERS (.LE. 5)',ZZ7X) 
CLOSE(1) 
CLOSE(2) 
STOP 
END 
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SOFTWARE SOURCES 


Ryan=McFarland FORTRAN 
609 Deep Valley Drive 
Rolling Hills Estates, CA 90274 
(213) 541 - 4828 


Microcompatibles, Inc. GRAFMATIC 
301 Prelude Drive 
Silver Spring, MD 20901 
(301) 593 - 5151 


Jewell Technologies, Inc. GRAFPLUS 
Seattle, WA 
(206) 937 - 1081 


West Coast Consultants CURVE DIGITIZER 
4202 Genesee Avenue, Suite 309 
San Diego, CA 92117 
(619) 565 - 1266 


Microsoft FORTRAN 
16011 NE 36th Way 
Box 97017 
Redmond, WA 98073 - 9717 


Prof. M. A. Morgan 
Code 62Mw 
Naval Postgraduate School 
Monterey, CA 93943 


LI. E. M: Connolly 


c/o Jd. J. Connolly 
208 Kent Drive 
Greenville, NC 27858 
(919) 752 - 2085 


OR 
С/ОЛС 5 лел 
228 Windsor Road 


Greenville, NC 27858 
(919) 756 - 4918 
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